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Resumen 



En esta Tesis presento un estudio espectroscópico detallado de una muestra de Regiones 
Hii Gigantes en galaxias observables con telescopios del Hemisferio Sur. 

Las Regiones Hll Gigantes proveen el vínculo entre regiones de formación estelar de menor 
escala y los violentos procesos de formación estelar que ocurren en las galaxias con brotes 
estelares intensos. Las Regiones Hii Gigantes (GHiiRs, de sus siglas en inglés) son objetos 
extendidos y luminosos, observados en los discos de galaxias espirales, irregulares y galaxias 
enanas compactas azules (BCDs, de su sigla en inglés). Las GHiiRs se forman debido a la 
presencia de un gran número de estrellas jóvenes y masivas cuyo flujo ultravioleta ioniza el gas 
que las rodea, y por lo tanto estas regiones indican una presencia de formación estelar violenta 
reciente o activa. El principal objetivo de este trabajo de Tesis es confirmar la presencia de 
Regiones Hii Gigantes en galaxias del Universo Local y realizar un estudio comparativo de 
GHllRs con distintas metalicidades y estados evolutivos. 

La identificación de Regiones Hii Gigantes fue posible realizarla mediante espectroscopia 
de alta resolución, la cual me ha permitido distinguir la presencia de distintas componentes 
cinemáticas en el gas ionizado, con diferentes velocidades y verificar la naturaleza supersónica 
de dichas componentes. 

Definir la naturaleza y condiciones físicas de los brotes de formación estelar masiva permite 
analizar las abundancias químicas que producen el enriquecimiento químico del Universo 
Local. Para poder caracterizar los procesos físicos que gobiernan la formación estelar masiva 
es necesario estudiar un variado rango de entornos y propiedades. Gracias al amplio rango 
espectral y la buena calidad de los datos he llevado adelante el estudio de las propiedades 
físicas del gas ionizado. El estudio de los parámetros físicos tales como densidad, temperatura, 
composición química, etc., comprende el análisis de las intensidades relativas de las líneas de 
emisión del gas. Este tipo de análisis resulta esencial para poder mejorar el estudio de brotes 
de formación estelar detectados a distancias cosmológicas. 

En esta Tesis se presenta el estudio realizado en seis regiones brillantes de las galaxias 
espirales, NGC 6070 y NGC 7479, y en cinco brotes de formación estelar de la galaxia enana 
compacta azul Haro 15. 

En el Capítulo 1 presento una breve introducción de las propiedades generales de las Re- 
giones Hii Gigantes y del conocimiento que se tiene sobre este tipo de objetos. En el Capítulo 
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2 propongo un método para proceder a un análisis detallado de la cinemática en aquellas 
regiones que presentan más de una componente revelada cinemáticamente en los espectros 
de alta resolución. En ese capítulo presento el análisis cinemático de las seis regiones estu- 
diadas en las dos galaxias espirales y se analiza cómo la presencia de múltiples componentes 
cinemáticas influye en la ubicación final en el plano L-cj. En el Capítulo 3 presento un análisis 
cinemático utilizando una metodología similar que en el capítulo anterior, pero ahora enfo- 
cado en los brotes de formación estelar violentos de la galaxia BCD Haro 15. Del análisis 
cinemático detallado en los brotes estudiados con más de una componente, se deriva el estu- 
dio de las propiedades físicas y abundancias químicas presentado en el Capítulo 4. En este 
capítulo, sumado a los datos espectroscópicos de alta resolución, presento datos de baja reso- 
lución espectral observados en la misma galaxia. Aplicando el método propuesto por el Grupo 
de Astrofísica del Departamento de Física Teórica de la Universidad Autónoma de Madrid, 
España (UAM), para la determinación de los parámetros físicos y abundancias, a partir de un 
programa generado por dicho grupo, realicé un cuidadoso y detallado estudio sobre las pro- 
piedades físicas, y abundancias químicas iónicas y totales del gas ionizado discriminado para 
cada componente cinemática y para el flujo global de las líneas de emisión. En el Capítulo 5 se 
presentan las conclusiones finales de la Tesis, haciendo un conciso resumen sobre los aportes 
más significativos, y se describen las líneas de trabajo a futuro que se pueden plantear como 
continuación de lo aquí expuesto. 
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Capítulo 1 

Generalidades 



Las Regiones Hii Gigantes (GHiiRs, sigla en inglés de Giant Hii Regions) son objetos 
extendidos y luminosos, con luminosidades en la emisión de Ha de alrededor de lO'"^ ergs~^, 
observados en los discos de galaxias espirales, irregulares y galaxias enanas compactas azules 
(BCDs, sigla en inglés de Blue Compact Dwarfs). Las GHiiRs se originan debido a la presen- 
cia de un gran número de estrellas jóvenes y masivas cuyo flujo ultravioleta ioniza el gas que 
las rodea y, por lo tanto, estas regiones indican una presencia de formación estelar violenta 
reciente o activa (Chu & Kennicutt, 1994). Esto último se ha podido verificar para el caso de 
la GHiiR visible más cercana, 30 Doradus en la Nube Mayor de Magallanes (ver Figura 1.1), 
con espectroscopia del núcleo del cúmulo estelar, conocido como R136, realizada con el Teles- 
copio Espacial Hubble (Massey &¿ Hunter, 1998). Tienen tamaños típicos de un centenar de 
parsecs de diámetro. Existen menciones en la literatura de tamaños de un orden de magnitud 
mayor, pero en general se refieren a regiones múltiples, las cuales pueden resolverse en varias 
componentes (Bosch et al., 2002). Las GHiiRs podrían ser clasificadas, de alguna manera, 
en una escala intermedia entre regiones de formación estelar de menor escala, como Orion en 
nuestra galaxia (ver Figura 1.2), y regiones con violentos procesos de formación estelar que 
ocurren en las galaxias con brotes estelares intensos (llamadas starburst en inglés) (Searle &¿ 
Sargent, 1972). 

Estas regiones albergan cientos de estrellas masivas que, desde el punto de vista químico 
y morfológico, son parcialmente responsables de los cambios en la estructura de las galaxias 
a través de una constante emisión de masa a alta velocidad (vientos) o durante uno de los 
procesos energéticos más impresionantes del Universo como son las explosiones de supernovas. 
La gran cantidad de fotones ultravioleta emitida por las estrellas OB jóvenes y calientes junto 
con los de las estrellas menos masivas, con menor nivel de emisión pero muy numerosas, tiene 
como consecuencia que en muchas zonas de estas regiones el gas, principalmente hidrógeno, 
está altamente ionizado, con temperaturas típicas del orden de 10000 K. 

Junto a las zonas de gas caliente ionizado, las GH iiRs contienen gas en otros dos estados. 
Por un lado, se encuentra en un estado molecular frío y denso. En este estado el gas puede 
ser re-generado en zonas más densas que pueden formarse por la compresión en esas regiones 
dinámicamente complejas. Es en tales nubes moleculares frías que la actual formación estelar 
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Figura 1.1: Imagen Ha de la Nebulosa Tarántula 30 Doradus, o NGC 2070. Imagen tomada 
con un telescopio refractor TEC-140 a f/7, en Hakos, Namibia, por Johannes Schedler - 
PantherObservatory. 

se está produciendo. Por otro lado, grandes volúmenes de gas caliente de baja densidad son 
producidos por la acción combinada de vientos y explosiones de SN que ocurren en la región. 
Un entorno tan complejo da lugar a más de un evento de formación estelar o a una nueva 
generación de estrellas. Uno de los principales propósitos para estudiar estas regiones radica en 
ampliar el conocimiento sobre la formación estelar dado que desde muchos aspectos astrofísicos 
es esencial para comprender desde la evolución de las galaxias, a través del enriquecimiento 
químico de las estrellas masivas, hasta la formación de los planetas, con el fin de entender 
las condiciones iniciales de los discos proto-planetarios. Es muy posible que la mayoría de 
las estrellas que se formaron alguna vez, y la mayoría de las que existen en la actualidad, lo 
hayan hecho en condiciones más similares a las de las Regiones Hii Gigantes o "starbursts" , 
que a las condiciones más apacibles en las que se forman hoy en días las estrellas en nuestra 
vecindad galáctica. 

En los últimos años ha habido grandes progresos para entender la formación estelar. Los 
modelos teóricos han tenido éxito en la reproducción de las principales características de las 
estrellas formadas en las nebulosas, tales como la función inicial de masa (IMF, sigla en inglés 
de Initial Mass Function). Los principales procesos físicos que pueden cambiar la masa final 
de las estrellas y que pueden establecer eficiencia en la formación estelar de las nebulosas son: 
la fragmentación durante el colapso, la acreción de masa, los vientos estelares y la radiación 
ultravioleta (Klessen, 2001; McKee & Tan, 2003; Tillcy &¿ Pudritz, 2005). Sin embargo, es 
muy complicado saber hasta qué punto estos procesos son capaces de modificar la función 




Figura 1.2: Imagen Ha de la Nebulosa (M42) de Orion. Imagen tomada con un telescopio 
refractor Takahashi Sky 90 a f/4.5, por Don McCrady 



inicial de masa. En la actualidad, hay dos modelos principales que compiten para explicar el 
origen de la IMF: i) el modelo que plantea que la IMF estelar se origina en la distribución 
de masas de los núcleos densos de las nubes (CMF, sigla en inglés de core mass function) 
(Alves et al., 2007; Simpson et al., 2008). Bajo esta hipótesis, sólo se necesita explicar por 
qué la turbulencia en las nubes produciría una distribución de masas de los núcleos igual a 
la distribución de masas de las estrellas. Se han construido teorías para determinar la CMF 
a partir de la turbulencia en el medio interestelar (Padoan & Nordlund, 2009; Hennebelle & 
Chabrier, 2008), y se utilizan simulaciones hidrodinámicas para comprobar si se verifican las 
predicciones de las teorías (Smith et al.,2009aa, 2009bb). ii) El otro modelo alternativo del 
origen de la IMF es el de la "acreción competitiva" de Bate &: Bonnell (2003). En este modelo, 
las estrellas adquieren su masa compitiendo por el material disponible para ser acretado. En 
este caso la turbulencia es poco relevante para las masas finales de las estrellas, siendo sólo 
el mecanismo que originalmente fragmenta a la nube. Este modelo ha sido criticado por 
requerir, aparentemente, condiciones demasiado "restrictivas" para las estrellas en formación, 
y que posiblemente sólo puedan darse en la formación de cúmulos estelares muy poblados. 
Como se puede ver, hoy en día las discusiones siguen abiertas en este tema. 

Por otro lado, el advenimiento de la generación de los grandes telescopios ha permitido 
estudiar y analizar las regiones de formación estelar no solo del Universo cercano, sino en 
galaxias distantes, como por ejemplo, las galaxias Hii que constituyen un potente indicador 
de distancia a escalas cosmológicas. Ya han sido identificadas numerosas galaxias Hii en los 
seguimientos realizados con el Sloan Digital Sky Survey (SDSS) (que en castellano significa 
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Figura 1.3: Dibujo esquemático que simula la geometría de las capas de mezcla turbulenta 
mostrando el gas frió y caliente y la capa fotoionizada lindera a la capa de mezcla. El gas 
caliente se mueve a una velocidad transversal Vf relativa a la capa de mezcla (Slavin ct al., 
1993). La figura pertenece al mismo trabajo. 



"Relevamiento Digital del Cielo Sloan"). 

Las galaxias Hii son galaxias con formación estelar violenta que producen una señal es- 
pectral fuerte e inconfundible la cual puede detectarse hasta z~3 (Lowcnthal et al., 1997). 
Melnick, Terlevich & Terlevich (2000) propusieron el uso de las galaxias Hii como una sonda 
cosmológica que resulta útil para medir la densidad de masa del Universo (ÍÍa/)- Sin embargo 
este método que es potencialmente muy poderoso, se basa en la calibración del punto cero de 
la regresión observada en el plano fundamental para Regiones Hii Gigantes Extragalácticas 
(su sigla en inglés es GEHiiRs). 



1.1. Características de las GHllRs 



Smith &¿ Weedman (1970) encontraron que los anchos de los perfiles de las líneas de emi- 
sión de las GHllRs corresponden a velocidades supersónicas en el gas. Dicha fuente cinemática 
observada permanece aún sin una clara explicación, pero permite identificar y distinguir una 
región Hii gigante de una mera aglomeración de regiones Hii clásicas por medio de espectros- 
copia de alta resolución. 

Melnick (1977) encontró una correlación estrecha entre la luminosidad de las GHiiRs y 
la dispersión de velocidades observada en el perfil de las líneas de emisión. Diversas inves- 
tigaciones siguieron a este descubrimiento y confirmaron la existencia de una correlación. 



1.1. características de las ghurs 




log(velQCitj' dispergían) 



Figura 1.4: Regresión en el plano Luminosidad (L) vs. Dispersión de Velocidades (a) de Bosch 
et al. (2002). La figura pertenece al mismo trabajo. 



aunque exhibían desacuerdos en las pendientes encontradas para dicha regresión (Tcrlcvich 
& Melnick, 1981; Roy et al., 1986; Bosch et al., 2002). Una gran variedad de mecanismos son 
los responsables del ensanchamiento que presentan los perfiles de las líneas de emisión en el 
espectro, lo que ha originado varías teorías paralelas para explicar dicho fenómeno, invocando 
dinámica de sistemas virializados (Tenorio- Tagle et al., 1993), superposición de múltiples bur- 
bujas de gas en expansión (Chu & Kennicutt, 1994), o turbulencia del mismo gas interestelar 
(Medina Tanco et al., 1997). Bosch y colaboradores obtuvieron fotometría CCD de banda 
angosta en líneas de emisión para un grupo de GH ilRs, analizando cómo incidía la incerteza 
de los fiujos publicados sobre las discrepancias de los valores hallados para la regresión en 
el plano Luminosidad (L) vs. Dispersión de Velocidades (a). Se halló una doble lectura, ya 
que los nuevos datos fotométricos no mejoraban la dispersión de los puntos en la regresión. 
Sin embargo, los autores realizaron un análisis más detallado de las regiones, considerando 
información sobre el estado evolutivo de algunas GHiiRs. Malumuth et al. (1996) y Drissen 
et al. (1993) analizaron la componente estelar de NGC 595, una de las GHiiRs de la galaxia 
espiral M33, y encontraron signos de evolución. Con espectroscopia Fabry-Perot de alta re- 
solución espacial y espectral de algunas regiones de M33 y MlOl, Muñoz Tuñón et al. (1995; 
1996) resolvieron el perfil supersónico global observado en varias componentes discretas, algu- 
nas con anchos subsónicos. De acuerdo al modelo de evolución propuesto por Muñóz-Tuñón 
et al. (1996), cuando las burbujas de gas sopladas por estrellas masivas dentro de las regiones 
Hii gigantes alcanzan el borde de la región, alteran el perfil global observado de la línea de 
recombinación del gas. Al descartar las regiones con evidentes signos de evolución, Bosch et 
al. (2002) obtuvieron una regresión muy estrecha en el plano L vs. a, aunque basada en una 
muestra muy reducida (ver Figura 1.4). Si se desea estudiar la relación entre luminosidad y 
dispersión de velocidades en una muestra de GHiiRs que tenga sentido estadístico es necesario 
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Figura 1.5: Dibujo esquemático que representa las ideas desarrolladas en el trabajo de West- 
moquette et al., (2009) sobre el estado del ISM en la galaxia M82. El dibujo esquemático 
pertenece a dicho trabajo y la composición es de autoría de Linda Smith. 



descubrir nuevas Regiones Hii Gigantes. 

Varios estudios han tratado el hecho de que un perfil simple Gaussiano puede no ser sufi- 
cientemente realista a la hora de ajustar las líneas de emisión observadas. Y en varios estudios 
se ha favorecido la presencia de una componente ancha omnipresente que podría estar expli- 
cando las alas del perfil integral (ver Figura 1.6) (Muñoz-Tuñón et al., 1996; Melnick et al., 
1999; Díaz et al., 1987; Terlevich et al., 1996; Méndez & Esteban, 1997; Homeier & Gallagher, 
1999; Hágele et al., 2007, 2009, 2010). Westmoquette et al., (2007a; 2007b) concluyen que la 
componente angosta representa el medio interestelar (ISM, sigla en inglés de Interestellar Mé- 
dium) ionizado y perturbado, que surge a través de la convolución de los efectos de agitación 
del gas debido al brote de formación estelar y los movimientos gravitacionales del material 
virializado. Y por otro lado, la componente ancha resulta del campo de velocidad altamente 
turbulento asociado con la interacción de la fase caliente del ISM con nudos de gas más frío 
(ver Figura 1.3), lo que genera capas de mezcla turbulenta (ver Figura 1.5). Sin embargo, 
otros estudios han apoyado la existencia de perfiles tipo "shell" con dos alas, azul y roja, a 
ambos lados de la componente principal (Chu &: Kcnnicutt, 1994; Rclaño & Bcckman, 2005; 
Rozas et al., 2006). 

Para verificar la naturaleza gigante de estas regiones Hii es necesario realizar espectros- 
copia de alta resolución. La espectroscopia échelle proporciona un medio para buscar diferentes 
componentes cinemáticas en las líneas de emisión del gas y analizar estos efectos, ya que al- 
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1.1. características de las ghurs 




4875 
Wavelength (Á) 



Figura 1.6: Ajuste de dos componentes en el perfil de emisión de íí/3 en una región Hii 
circunnuclear de la galaxia NGC 3310 (Hágclc ct al., 2010). La figura pertenece al mismo 
trabajo. 



canza la resolución espectral necesaria tanto para resolver los anchos de los perfiles de las 
líneas de emisión como para determinar las dispersiones de velocidades intrínsecas del gas 
ionizado y verificar si estas son supersónicas (>13 kms~^). 

La principal fuente de información sobre la metalicidad en todo tipo de objetos, desde 
los más cercanos hasta los más distantes, lo constituyen las regiones Hll, desde las menos 
brillantes en el disco de la Galaxia, pasando por las GEHllRs hasta las galaxias Hll. Sus 
espectros presentan líneas de emisión brillantes las cuales son visibles en todas las regiones 
donde se hayan producido episodios de formación estelar recientes. El análisis de dichos es- 
pectros nebulares constituye el método más efectivo, a grandes distancias, para determinar 
las abundancias químicas de elementos tales como el helio, nitrógeno, oxígeno, neón, argón 
o azufre que dispongan de líneas de emisión en el rango óptico correspondientes a diferentes 
estados de ionización. El conocimiento detallado de estas abundancias resulta esencial para 
comprender la evolución del gas y de las estrellas, y esta información permite ampliar el co- 
nocimiento que se tiene acerca de la evolución química de las galaxias en el Universo Local. 
Con el advenimiento de los grandes telescopios, es posible disponer de datos para objetos 
menos brillantes y/o más lejanos, que acrecienta el conocimiento que se tiene del Universo 
más lejano y por tanto más joven. 
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1.2. Motivación y Objetivos 

Realizando una inspección en la lista de GHiiRs conocidas a comienzos de los años 2000, 
habíamos encontrado que todas pertenecían al hemisferio norte, salvo 30 Doradus, en la Nube 
Mayor de Magallanes. Feinstein (1997) había estudiado la función de luminosidad de unas 10 
galaxias del cielo austral, entre las cuales nosotros pudimos identificar al menos 40 regiones 
con luminosidades comparables a las de una GHiiR. Ocho de esas galaxias presentan regiones 
Hii con brillos superiores a 10^*^ erg s~^, involucrando un total de 31 candidatos a RHliG. 
Si se extiende la muestra a las regiones con luminosidades mayores a 10^^ erg s~^, el número 
crece hasta unos 70 objetos. Para comenzar a hacer un estudio de GHiiRs con un sentido 
estadístico, esta muestra de regiones Hii era ideal como puntapié inicial. Es así que de la 
muestra de Feinstein, seleccionamos las regiones más luminosas y estas pasaron a ser nuestros 
primeros candidatos a estudiar. En Firpo et al. (2005) obtuvimos datos de 4 candidatos, todas 
regiones H ii muy luminosas, en las galaxias NGC 7552 y NGC 2997 utilizando el espectrógrafo 
Magellan Inamori Kyocera Echelle (MIKE) del Telescopio Magallanes, LCO. La alta relación 
señal-ruido sumada a la resolución del échelle nos permitió resolver el perfil de las líneas 
de emisión y calcular la dispersión de velocidades del gas ionizado. Esto se hizo midiendo 
el ancho observado del perfil, que luego fue corregido por ensanchamiento térmico, y por el 
ancho intrínseco del perfil instrumental, para obtener, así el ensanchamiento verdadero debido 
al comportamiento cinemático de la Región Hii. De ese análisis descubrimos la naturaleza 
gigante de tres de las cuatro regiones candidatas y pudimos confirmar que el mecanismo 
de excitación de las mismas se debía a fotoioniz ación por estrellas. Además, encontramos 
que las regiones Hii gigantes descubiertas se ubican, dentro de las incertezas, en la regresión 
esperada para este tipo de objetos, en el plano L-a. Uno de los principales resultados, presente 
en los datos publicados por Firpo et al. (2005) y del análisis de 5 regiones Hii observadas 
posteriormente, es que todas las Regiones Hii Gigantes de las galaxias de la muestra de 
Feinstein presentan signos de evolución química y es necesario observar regiones más jóvenes. 

Entonces, con el objetivo de detectar regiones más jóvenes que permitan realizar un estudio 
comparativo de GHiiRs a distintas metalicidades y estados evolutivos, seleccionamos brotes 
de formación estelar intensa a partir de fotometría de líneas de emisión de galaxias más 
jóvenes, tales como las enanas compactas azules estudiadas por Cairos et al. (2001) y Gil de 
Paz et al. (2003). Detectamos que varios de los brotes observados en dichas galaxias presentan 
características de GHiiRs. De esta manera seleccionamos los brotes más brillantes en algunas 
galaxias de baja metalicidad, siendo estos los nuevos candidatos a estudiar. 

Mediante espectroscopia échelle de alta resolución y espectroscopia de dispersión simple 
obtenidas en distintos telescopios del Observatorio de Las Campanas (Chile), estudiamos 
los brotes más brillantes en diversos objetos. La muestra incluye galaxias espirales de alta 
metalicidad como NGC 6070 (Figura 1.7), NGC 7479 (Figura 1.8) (Firpo et al., 2010) y la 
galaxia de baja metalicidad HarolS (Figura 1.9) del catálogo de Cairos et al. (2001) (Firpo 



1.2. MOTIVACIÓN Y OBJETIVOS 




Figura 1.7: Galaxia NGC 6070. Esta fue la primer imagen (imagen de prueba) tomada por el 
Sloan Digital Sky Survey el 10 de junio de 1998. 



et al., 2011). 

A partir del ajuste de un perfil simple Gaussiano a las líneas de emisión, observamos un 
residuo presente en las alas de los perfiles en varias líneas de emisión, evidente en la línea 
Ha y que se confirma en otras líneas de emisión. Para determinar el número óptimo de com- 
ponentes Gaussianas y ajustar cada perfil de línea, y con el fin de minimizar las incertezas 
en la dispersión de velocidades del gas ionizado, consideramos la información proporcionada 
en la literatura a partir de la gran variedad de estudios propuestos para interpretar la exis- 
tencia del ensanchamiento en las líneas de emisión. Además de la componente relativamente 
angosta y siempre que fue posible, evaluamos la eventual presencia de una componente ancha 
supersónica medida en el perfil de la línea de emisión de la GHiiR, o la presencia de dos 
componentes simétricas de baja intensidad en el ajuste de los anchos de los perfiles de línea 
de emisión observados. 

En esta Tesis, confirmando lo ya mencionado en la literatura, presento evidencia de la 
existencia de diferentes componentes cinemáticas en las líneas de emisión del gas. Además, en 
algunas regiones pudimos detectar la presencia de más de una sub-componente angosta, cuyas 
velocidades radiales relativas no pueden explicarse mediante fiuctuaciones en la rotación del 
disco galáctico (Firpo ct al., 2010, 2011). 

Para avanzar en el análisis de las propiedades físicas de este tipo de objetos y en las 
diferentes componentes encontradas, debemos poder determinar temperaturas, densidades y 
abundancias químicas del gas. El análisis de dichos parámetros físicos puede llevarse adelante 
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Figura 1.8: Galaxia NGC 7479. Tomada por Nik Szymanek y obtenida con el Telescopio 
Faulkes Norte, Maui, copyright FTLLC 




Figura 1.9: Galaxia Haro 15. Imagen Ha tomada con la cámara Wide Field Planetary Camera 
2, obtenida del Multimission Archive del Space Telescope Science Institute (MAST). 
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mediante el estudio de las intensidades relativas de las líneas de emisión del gas. Para ello 
necesitamos también una amplia cobertura espectral. La calidad única de nuestros datos 
permitió realizar un análisis detallado de las características del gas ionizado discriminado para 
cada componente y para el flujo global de las líneas. Para llevar esto a cabo, se ha aplicado 
el método propuesto por el Grupo de de Astrofísica del Departamento de Física Teórica de la 
Universidad Autónoma de Madrid, España, para la determinación de los parámetros físicos y 
abundancias. 

En esta Tesis presento resultados sobre la estimación de temperaturas, densidades y abun- 
dancias químicas discriminadas por las diferentes componentes cinemáticas del gas. 
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Giant Hll Regions in NGC 7479 & NGC 6070 

Firpo, V., Bosch, G., Hágele, G. F., and Morrell, N., 2010, MNRAS, Vol 406, Issue 2, 

pp. 1094 

2.1. Resumen de la publicación 

En este Capítulo se presentan nuevos resultados de la búsqueda de Regiones H ii Gigantes 
en dos galaxias visibles desde el hemisferio sur: NGC 7479 y NGC 6070. Usando espectros 
de alta resolución, obtenidos con diferentes instrumentos del Observatorio Las Campanas, 
Chile, se han podido resolver los anchos de los perfiles de las líneas de emisión y determinar 
la dispersión de velocidades intrínseca del gas ionizado. Se ha detectado que los anchos de los 
perfiles de emisión corresponden a dispersiones de velocidades supersónicas en las seis regio- 
nes Hii observadas. Se ha encontrado que todos ellos muestran, al menos, dos componentes 
cinemáticas diferentes: una componente relativamente angosta (entre 11 y 22 kms~^)yotra 
componente más ancha (entre 31 y 77 kms~^). Dos de las regiones estudiadas muestran un 
perfil complejo en las líneas de emisión angostas de los diferentes iones, el cual pudo ser dividi- 
da en dos componentes con distintas velocidades radiales. Mientras tanto, el ensanchamiento 
de las alas del perfil global pudo ser ajustado con una componente ancha en casi todos los 
perfiles de las líneas de emisión observadas. En una región este residuo pudo ser reproducido 
como perfiles tipo "shell" (cascaras) en expansión debido a la presencia de dos componentes 
separadas a ambos lados de las componentes principales. Además, se ha analizado el impac- 
to que tiene la presencia de varias componentes sobre la ubicación de la región Hii en el 
plano log(L) — log{a). Aunque la distribución global confirma la presencia de una regresión, 
la ubicación precisa de las regiones en el plano depende fuertemente de las propiedades de las 
componentes derivadas en el ajuste del perfil de emisión. 
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2.2. Fundamentos básicos 

Con el fin de verificar la naturaleza de nuestras regiones candidatas a Regiones Hii Gi- 
gantes se obtuvieron espectros de alta resolución para poder medir los anchos en los perfiles 
de las líneas de emisión y así estimar si la dispersión de velocidades es, en efecto, supersónica. 
Para determinar el número óptimo de componentes Gaussianas que mejor se ajusta a cada 
perfil de línea, y con el fin de minimizar las incertezas en la dispersión de velocidades del 
gas ionizado, se ha considerado la información disponible en la literatura. Algunos trabajos 
previos han tratado el hecho de que un perfil simple Gaussiano puede no ser suficientemente 
realista a la hora de ajusfar las líneas de emisión observadas. Varios estudios han favorecido 
la presencia de una componente ancha omnipresente que podría estar explicando las alas del 
perfil integral (Muñoz- Tuñón et al., 1996; Melnick et al., 1999; Díaz et al., 1987; Terlevich 
et al., 1996; Méndez & Esteban, 1997; Homeier & Gallagher, 1999; Hágele et al., 2007, 2009, 
2010, entre otros). En cambio, otros estudios han apoyado la existencia de perfiles tipo shell 
con dos alas, azul y roja, a ambos lados de la componente principal (Chu & Kennicutt, 1994; 
Relaño & Beckman, 2005; Rozas ct al., 2006). 

En este trabajo se presentan los datos échelle obtenidos en el Observatorio Las Campanas 
(LCO) de seis regiones candidatas a GHiiRs en las dos galaxias mencionadas. Se ha determi- 
nado la dispersión de velocidades y evaluado la posible presencia de una componente ancha 
o dos componentes simétricas a ambos lados de las componentes angostas en el ajuste de 
los anchos de los perfiles de líneas de emisión observados. En la publicación a continuación 
se presentan las observaciones y la reducción de los datos en la Sección Observations and 
Reductions; en la Sección Analysis OF line profiles los detalles del anáfisis realizado 
sobre los perfiles de las líneas de emisión y se analiza cada caso individualmente; en la Sección 
Relation between Ha Luminosities and Velocity Dispersión se anafiza el impacto 
que tiene la presencia de varias componentes sobre la ubicación de la región H ii en el plano 
log(L) — log(a"); por último se incluye un resumen y las conclusiones del trabajo en la Sección 
nombrada como SUMMARY AND Conclusions. 

2.3. Herramientas utilizadas 

De acuerdo a la denominación canónica de GHllRs, los candidatos han sido nombrados 
siguiendo el rango de brillo publicado en el catálogo de Feinstein (1997). La alta resolución 
espectral se obtuvo con el espectrógrafo échelle del telescopio de 100" du Pont, del Observa- 
torio Las Campanas (LCO), observados en 2006. El rango espectral abarcó desde 3800 hasta 
9500 A. Las condiciones de observación fueron buenas con un seeing de un segundo de arco. 
Con el fin de minimizar la contribución del ruido en la lectura del espectro final, se aplicó un 
binning 2x2 al CCD. La resolución espectral alcanzada en los espectros échelle (obtenido con 
una ranura de 1 segundo de arco de ancho) es de R~25000: AA=0.25A a los A 6000A, medida 
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a partir del ancho a potencia mitad (FWHM) de las líneas de la lámpara de comparación de 
Thorio- Argón. Esto se traduce en una resolución en velocidad radial de ~12 kms~^. 

Durantes las noches de observación se tomaron espectros de estrellas estándares, HR 7950, 
HR4468, y una estándar espectrofotométrica de CALSPEC, FeigellO, con el propósito de 
calibrar en flujo los espectros. Además, se tomaron todas las noches espectros de comparación 
de Th-Ar, milky ñats (flats de cielo obtenidos con un difusor durante la tarde), y bias. El 
tiempo de exposición para las estrellas estándares fue de cinco segundos para las estrellas 
brillantes HR7950 y HR4468, y 1200 segundos para FeigellO. 

También se han obtenido espectros de alta resolución para las tres regiones Hii de la 
galaxia NGC 6070 utilizando el espectrógrafo échelle doble, Magellan Inamori Kyocera Eche- 
lle (MIKE) del Telescopio de 6.5-m, Magellan II (Clay), de LCO, en julio de 2004. No se 
aplicó ningún binning al CCD de 2Kx4K y la ranura utilizada fue de 1 segundo de arco. 
La resolución espectral medida en los espectros de MIKE es de 11 kms~^, muy similar a 
la obtenida en los espectros de du Pont. El tiempo de exposición para la estrella estándar 
espectrofotométrica, BD+28D4211, fue de 900 segundos. También se tomaron lámparas de 
comparación de Th-Ar, milky flats y bias para calibrar los datos. 

El análisis de los datos se realizó con el software IRAF^. Luego de la resta de bias y 
correcciones por flat fields con los Milky Flats, las imágenes bidimensionales fueron corregidas 
por rayos cósmicos con la tarea cosmicrays de IRAF, la cual detecta y elimina los rayos 
cósmicos con un algoritmo de cociente de flujos. Los datos corregidos se redujeron con las 
rutinas de IRAF siguiendo el mismo procedimiento descripto en Firpo et al., (2005). 

2.4. Metodología en la cinemática 

Para determinar las velocidades radiales y las dispersiones de velocidades del gas ionizado, 
se ha medido la longitud de onda central y el ancho en varias líneas de emisión intensas. La 
velocidad radial de cada línea de emisión se ha obtenido de su longitud de onda central 
determinada ajusfando el perfil Gaussiano, y sus errores han sido proporcionados por el error 
obtenido en el ajuste realizado con la tarea ngaussf it de IRAF y teniendo en cuenta el error 
medio cuadrático de la calibración en longitud de onda. La dispersión de velocidades intrínseca 
(o'int) de cada línea de emisión ha sido calculada como se describe en Firpo et al., (2005) y 
detallada aquí en la Sección 2.5, donde se explica que la dispersión de velocidades intrínseca 
(cmí) es corregida por las contribuciones instrumental y térmica de cada línea de emisión. 

Para ajusfar las múltiples componentes en los perfiles observados se utilizó la tarea 
ngaussf it de IRAF. Esta tarea realiza un ajuste Gaussiano iterativo para varios perfiles 
de líneas espectrales. La tarea permite seleccionar distintos parámetros a ser ajustados y 
ofrece una estimación de la incerteza de los parámetros derivados. La tarea necesita valores 



^Image Rcduction and Analysis Facility, distribuido por NOAO, operado por AURA, Inc., bajo contrato 

con NSF. 

15 

Tesis de Doctorado Verónica Firpo 



16 CAPITULO 2. REGIONES Hll GIGANTES EN NGC 7479 & NGC 6070 

iniciales de los coeficientes, los cuales pueden ser especificados a través de tablas generadas 
previamente. 

En los casos en que más de una componente se hace evidente, hemos seguido un proce- 
dimiento iterativo en el que sólo se le permite a la tarea ajustar un subconjunto limitado de 
parámetros a la vez. Esto se hace para limitar el universo de posibles soluciones en el ajuste, 
y haciendo uso de la información disponible en cada línea de emisión. Luego de obtener una 
estimación inicial de los parámetros de la segunda componente, son fijados la amplitud, el 
centro y el ancho de cada componente, dejando los parámetros de la otra componente sin 
tocar. Una vez que la longitud de onda central de cada componente se conoce relativamente 
bien, se establece ésta como parámetro fijo y se permite que la tarea realice un ajuste final 
sobre la amplitud y el ancho en todas las componentes y al mismo momento. 

En todos los casos los ajustes de los perfiles han demostrado la presencia de un residuo en 
las alas de las emisiones, bien evidente en la línea intensa Ha. Por lo tanto, se ha evaluado la 
existencia de una componente ancha o dos componentes tipo shell en el perfil de la línea. En 
la primera hipótesis, se ha introducido una componente con una anchura inicial que es tres 
veces el ancho de la componente angosta (Hágclc ct al., 2007). Se ajustaron las posiciones 
de los picos para las componentes angostas y, a continuación, los anchos de los perfiles. Los 
parámetros de las componentes anchas fueron dados luego, antes de iniciar el procedimiento 
iterativo nuevamente. En el segundo escenario, el ajuste del residuo presente en las alas del 
perfil, fue realizado siguiendo un procedimiento iterativo similar, pero en este caso se añade 
una componente angosta corrida al azul y una componente angosta corrida hacia el rojo. 
La validación de la multiplicidad en el perfil y el ensanchamiento es controlado sobre las 
diferentes líneas de emisión para cada región, aunque el análisis más confiable proviene del 
perfil de las líneas más intensas, como la línea de recombinación Ha. 

De la muestra de seis GHiiRs se ha encontrado que todas ellas muestran evidencia de un 
residuo presente en los perfiles de las líneas de emisión. Sólo en una de las regiones, donde 
no pudo ser ajustada una componente ancha, se han ajustado dos componentes angostas 
simétricas separadas en velocidad respecto a las componentes principales; en el resto de las 
regiones el ajuste de una componente ancha dio resultados satisfactorios. En la publicación se 
discuten los resultados agrupándolos según los tipos de resultados de los perfiles y por región. 

2.5. Dispersión de velocidades 

El parámetro físico que se desea medir es la dispersión de velocidades intrínseca del gas io- 
nizado {(Tint)- Sin embargo, el ancho del perfil observado (o"o) está afectado por la contribución 
del movimiento térmico aleatorio {at) y del perfil instrumental (cj) dado por el espectrógrafo. 
Por lo tanto, dispersión de velocidades verdadera estará dada por: 

2 2 2 2 

(yint = (To - <Tj - (T^ 
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suponiendo que las regiones observadas tienen una temperatura cinética típica T = W^K 
y que el perfil instrumental (ai) se puede determinar a partir del ajuste Gaussiano de las 
líneas de emisión angostas de la lámpara de Th-Ar. Los errores en la determinación de la 
dispersión de velocidades son calculados utilizando los errores de observación en el perfil 
observado ao y considerando errores insignificantes en c7¿ y at- Debido a la alta metalicidad 
de las regiones Hii en galaxias espirales, las temperaturas son probablemente inferiores a 
lO^K (Díaz et al., 2007). En estas regiones también se ha hecho la suposición del caso más 
extremo donde la temperatura electrónica sea de unos 5000 K, para analizar como afecta este 
cambio de temperatura a la contribución de la dispersión de velocidades, pero los resultados no 
arrojaron grandes diferencias respecto de suponer T = lO^K. Para una temperatura cinética 
de T=10^K, la velocidad del sonido es de unos 13 kms~^. Se dice que la dispersión de 
velocidades es supersónica cuando es mayor o igual a ese valor. En este caso se está frente a 
una GHiiR. Valores menores, corresponden a regiones Hii clásicas. 
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ABSTRACT 

We present new results from our search for giant Hn regions in galaxies visible from the 
Southern tiemisphere. In this work we study two galaxies: NGC 7479 and NGC 6070. Using 
high-resolution spectra, obtained with different Instruments at Las Campanas Ohservatory, 
we are ahle to resolve the emission-line profile widths and determine the intrinsic velocity 
dispersión of the ionized gas. We detect profile widths corresponding to supersonic velocity 
dispersions in the six observed H n regions. We find that all of them show at least two distinct 
kinematical components: a relatively narrow feature (between 1 1 and 22 km s^' ) and a broader 
(between 3 1 and 77 km s^^) component. Two of the regions show a complex narrow profile in 
all ion lines, which can be further split into two components with different radial velocities. 
Whereas the wing broadening of the overall profile can be fitted with a low-intensity broad 
component for almost all profiles, in one región it was better reproduced by two sepárate 
shell-like wings. We have analysed the impact that the presence of múltiple components has 
on the location of the H ii regions in the log(L)- log((T ) plañe. Although the overall distribution 
confirms the presence of a regression, the precise location of the regions in the plañe is strongly 
dependent on the components derived from the profile fitting. 

Key words: H n regions - galaxies: individual: NGC 7479 - galaxies: individual: NGC 6070 - 
galaxies: starburst. 



1 INTRODUCTION 

The most extended and luminous objects which are observed on the 
discs of spirals, in irregulars and in starburst galaxies are called giant 
extragalactic H ii regions (GH nRs). GH iiRs are places of very active 
star formation and they are characterized by the large emission of 
ultraviolet photons from a large number of young and massive 
stars. Owing to the large number of ionizing photons most of the 
sunounding gas is ionized, and the recombination process provides 
a strong signature of recent or ongoing massive star formation. 

The emission-line profile widths of the giant H ii regions imply 
the existence of supersonic motions in the gas (Smith & Weedman 
1 970) . Furthermore, Terlevich & Melnick (1981) f ound a correlation 
between the gas velocity dispersión {a) and the total luminosity (L) 
emitted in the respective line. Hippelein (1986) and Arsenault & Roy 
(1988) confirmed the existence of such a regression, but there was 
no agreement on the valúes derived for its slope so no conclusión 
could be reached regarding the origins of this motion. Fuentes- 
Masip et al. (2000) and Bosch, Terlevich & Terlevich (2002) using 
different techniques on different samples agreed on a slope cióse 
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to a valué of 4, which seemed to favour the gravitational source of 
energy for the observed supersonic motion. Even if the nature of this 
behaviour is not fuUy understood, it still allows us to distinguish 
giant H II regions from an agglomeration of classical H ii regions 
by means of high-resolution spectroscopy. GH iiRs provide the link 
between small-scale star-forming regions, such as Orion in our 
Galaxy and violent star formation processes taking place in starburst 
galaxies. 

In Firpo, Bosch & Morrell (2005) we have confirmed the giant 
nature of three candidatos to giant H ii regions in the southem sky 
which were identified as very luminous Hii regions by Feinstein 
(1997). 

Here we continué our search and detailed analysis of GH iiRs 
in local Universe galaxies. We have selected from the Feinstein 
catalogue, a sample of the brightest Hii regions from two spiral 
galaxies: NGC 7479 and NGC 6070. 

NGC 7479 is a barred spiral galaxy, at a distance of 31.92 Mpc 
(Rozas et al. 1999), classified as SB(s)c by de Vaucouleurs et al. 
(1991). Devereux (1989) and Martinet & Friedü (1997) have clas- 
sified it also as a starburst galaxy, whereas Keel (1983) and Ho, 
Filippenko & Sargent (1997) have identified it as a low-ionization 
nuclear emission-line región (LINER) and a Seyfert 1.9 galaxy, 
respectively. Many studies have been performed on this galaxy ow- 
ing to the presence of a certain asymmetry in the spiral structure 
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(Quillen et al. 1995) that leads several authors to suggest this might 
be owed to a recent merger event (Laine 1996; Laine & Heller 
1999). This scenario can also account for perturbations in the ve- 
locity field (Laine & Gottesman 1998) and the latest minor merger 
model indicates that the remnant may be situated within the bar 
(Laine 2001) as no remnant of the merger could be identified in 
optical images (Saraiva & Benedict 2003). 

NGC 6070, at a distance of 29.8 Mpc (García-Gómez, 
Athanassoula & Barbera 2002), is classified as SA(s)cd by de 
Vaucouleurs et al. (1991) although Grosb0l, Patsis & Pompei (2004) 
State that a bar with two main arms can be identified. Pohlen & 
Trujillo (2006) further studied the bar structure classifying it as 
Type II.o.CT according to Erwin, Pohlen & Beckman (2008). The 
galaxy shows a normal arm distribution, with two arms that break 
into several on each side of the central región (see Fig. 2). The disc 
dynamics was studied by Márquez et al. (2002), and García-Gómez 
et al. (2002) studied the distribution of H ii regions in the galaxy. 

In order to verify the nature of our giant H ii región candidates, 
we obtained high-resolution spectra to measure the emission-line 
proñle widths and to estímate if the velocity dispersión is indeed 
supersonic. To determine the optimum number of Gaussian com- 
ponents to fit each Une profile in order to minimize the uncertainty 
in the velocity dispersión of the ionized gas, we have considered 
information from other studies available in the literature. Previ- 
ous papers have dealt with the fact that a single Gaussian profile 
might not be realistic enough when fitting the observed emission 
Unes. Several studies favour the presence of an omnipresent broad 
component which explains the integral profile wings. They include 
Melnick, Tenorio-Tagle & Terlevich (1999) who identified a broad 
component in the recombination lines in the central región of the 
30 Doradus nébula, Díaz et al. (1987) and Terlevich et al. (1996) 
in the M33 giant Hii región NGC 604, Méndez & Esteban (1997) 
in four Wolf-Rayet galaxies and Homeier & Gallagher (1999) in 
the starburst galaxy NGC 7673. More recently, Hagele et al. (2007, 
2009, 2010) and Hagele (2008) found that the best Gaussian fits 
in circumnuclear star-forming regions involved the presence of 
broad and narrow components for the emission lines of the ionized 
gas. 

A variety of line broadening mechanisms have been proposed 
to interpret the existence of the broad supersonic component mea- 
sured in the emission-line profile of GH iiRs. Dynamics of virialized 
systems (Tenorio-Tagle, Muñoz-Tuñón & Cox 1993), superposition 
of múltiple gas bubbles in expansión (Chu & Kennicutt 1994) or 
turbulence of the same interstellar gas (Medina Tanco et al. 1997). 
More recently, Westmoquette et al. (2007a,b) concluded that the 
narrow component represents the general disturbed ionized inter- 
stellar médium (ISM), arising through a convolution of the stirring 
effects of the starburst and gravitational virial motions. On the other 
hand, the broad component results from the highly turbulent veloc- 
ity field associated with the interaction of the hot phase of the ISM 
with cooler gas knots, setting up turbulent mixing layers. 

Other studies support the existence of shell-like profiles on 
the blue and red wings of the main feature as found by Chu & 
Kennicutt (1994), who state that the integrated line profiles of the 
entire 30 Dor región may be fitted with a broad Gaussian with 
low-intensity wings. Relaño & Beckman (2005) claim that a large 
fraction of the Hii regions in three barred spiral galaxies shows 
high-velocity and low-intensity features in its integrated line pro- 
files and Rozas et al. (2006) find that the great majority of their 
sampled H ii regions show that the brightest principal component is 
accompanied by low-intensity components symmetrically shifted in 
velocity. 
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Figure 1. Ha iraage with the GHiiR candidates observed in NGC 7479, 
identified by circles. Ha images were kindly made available by C. Feinstein. 



In this paper we present echelle data obtained at Las Campanas 
Observatory (LCO) of six GHiiRs candidates in the two men- 
tioned galaxies. We have determined the velocity dispersión and 
we have evaluated the possible presence of a broad component or 
two symmetric low-intensity components in the fit with the ob- 
served emission-line profile widths. In Section 2 we present the 
observations and the data reduction. Section 3 details the analysis 
performed on the emission-line profiles and discusses each conspic- 
uous case individually. The summary and conclusions of this work 
are included in Section 4. 



2 OBSERVATIONS AND REDUCTIONS 

According to canonical denomination for GHiiRs we have la- 
belled our candidates following their brightness ranking from Fe- 
instein catalogue. NGC 7479 1 is therefore the brightest región in 
NGC 7479, NGC 7479 II is the second brightest and so on. Ob- 
served regions from NGC 7479 and NGC 6070 are shown in Figs 1 
and 2, respectively. High-resolution spectra were obtained with an 
echelle spectrograph at the 100-inch du Pont Telescope, LCO, be- 
tween 2006 July 19 and 22. The spectral range covered by the 
observations was from 3800 to 9500 A. Observing conditions were 
good with 1 arcsec seeing and photometric sky. A 2 x 2 binning 
was applied to the CCD in order to minimize readout contribution 
to the final spectrum noise. The spectral resolution achieved in our 
du Pont Echelle data (obtained with 1 arcsec effective slit width) is 
i? ~ 25 000: Ak = 0.25 Á at A.6000Á, as measured from the full 
width at half-maximum (FWHM) of the Th-Ar comparison lines. 
This translates in a resolution of ~ 12 km s"' . 

The standard stars, HR7950 and HR4468, and the CALSPEC 
spectrophotometric standard star, Feige 110 (Bohlin, Dickinson & 
Calzetti 2001), were also observed for flux calibration purposes. In 
addition, Th-Ar comparison spectra, milky flats (sky flats obtained 
with a diffuser, during the afternoon) and bias frames were taken 
every night. The exposure time for standard stars was 5 s for the 
bright HR stars and 1200 s for Feige 1 10. 
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Figure 2. Same as Fig. 1 but for NGC 6070. 

Table 1. Journal of observations for the candidate GH iiRs and flux standard 
stars. The candidate regions have labelled following their brightness ranking 
from Feinstein catalogue (see Feinstein 1997). 



Telescope 


Región 


Date 


Exp. 


Secz 


Standard 


du Pont 


74791 


2006 July 19 


1800 


1.6 


FeigellO 




7479 II 


2006 July 20 


1800 


1.5 


HR7950 




7479 III 


2006 July 20 


1800 


1.4 


HR7950 




60701 


2006 July 21 


1800 


1.2 


HR4468 




6070 II 


2006 July 21 


1800 


1.2 


HR4468 




6070 IV 


2006 July 21 


1800 


1.2 


HR4468 


Clay 


60701 


2004 July 12 


1800 


1.3 


BD+28D42I1 




6070 II 


2004 July 12 


1800 


1.2 


BD+28D42I1 




6070 IV 


2004 July 12 


1200 


1.7 


BD+28D42I1 



We also obtained high-resolution spectra for the three Hii re- 
gions in the NGC 6070 galaxy using the echelle double spectrograph 
Magellan Inamori Kyocera Echelle (MIKE) at the 6.5-m Magellan 
II (Clay) Telescope, LCO, in 2004 July. No binning was applied to 
the 2k X 4k CCD detector and a I-arcsec slit was used. The spec- 
tral resolution measured on our MIKE spectra is 1 1 km s"' , which 
is quite similar than that obtained at du Pont. The exposure time 
for CALSPEC spectrophotometric standard star BD+28D42I 1 was 
900 s. Th-Ar comparison lamps, milky flats (infernal flat field lamp 
with a diffuser glass slide) and bias frames were used to calíbrate 
these data. Table 1 lists the journal of observations, together with 
the exposure times, airmasses and the standard stars that were used 
for flux calibrating each región. 

The data analysis was carried out with iraf' software. After bias 
subtraction and flat-field corrections with Milky Flats, the bidimen- 
sional images were corrected for cosmic rays with the task cosm- 
ICRAYS which detects and removes cosmic rays using a flux ratio 
algorithm. The corrected data were reduced by iraf routines fol- 
lowing similar procedures to those described in Firpo et al. (2005). 



' Image Reduction and Analysis Facility, distributed by NOAO, operated by 
AURA, Inc., under agreement with NSF. 



We also compared the red end of the wavelength-calibrated spec- 
tra with the night-sky spectrum by Osterbrock et al. (1996). This 
turned out to be a very reliable confirmation of the goodness of the 
wavelength solution, and we were able to check that differences 
between our wavelengths and the sky line wavelengths were below 
0.05 Á. 

We performed flux calibrations using the observed spectropho- 
tometric standard stars as described in Firpo et al. (2005). The 
observed standard stars had their fluxes tabulated every 16 A, and 
the amount of defined intervals within an echelle order ranged from 
four to 12, depending on the quality of the spectrum. 

For the Magellan spectra, and in the case of the NGC 7479 1 
región du Pont data, we were able to obtain CALSPEC spectropho- 
tometric standard stars (Bohlin et al. 2001), whose fluxes were tabu- 
lated every 2 A. These stars are ideal for calibrating high-resolution 
echelle spectra, although their relatively low brightness (Feige 110, 
V = 11.83 and BD+28D4211, V = 10.51) makes them time- 
consuming targets for 2-m class telescopes. 

3 ANALYSIS OF LINE PROFILES 

Making use of the known redshifts z = 0.006685 for NGC 6070 
(Spríngob et al. 2005) and z = 0.007942 for NGC 7479 (Haynes 
et al. 1998) we identify the hydrogen recombination lines, such as 
Ha and H/!, and coUisionalIy excited lines, such as [Nii] AA6548, 
6584 A, [S ii] AA67I7, 6731 A, present in the spectra. These strong 
lines are used to analyse the structure of the profiles as they allow 
US to verify the existence of more than one component. The adopted 
laboratory wavelengths are taken from work by García-Rojas et al. 
(2005). 

To determine the radial velocities and the velocity dispersions of 
the ionized gas, we measure the central wavelength and width of 
several emission lines. The radial velocity of each emission line is 
obtained from its central wavelength determined from the Gaussian 
profile fitting, and their errors are given by the fitting error provided 
by the resampling done within the ngauss task of iraf and taking 
into account the rms of the wavelength calibration. The true veloc- 
ity dispersión {a) of each emission line is calculated as described 
in Firpo et al. (2005). The width of the observed profile (ffo) is 
affected by the contribution of thermal random motions (crt) and 
the instrumental profile {a\). Then, we can obtain the true velocity 
dispersión as a^ = a^ — a^ ~ a} . As in Firpo et al. (2005), a typical 
kinetic temperature T = lO'' K is assumed, and the instrumental 
profile ((Ti) is very well approximated by a single Gaussian func- 
tion. The velocity dispersión errors have been calculated using the 
observational errors in the observed profile and assuming negligible 
errors in (Tí and ai. Owing to the generally high metallicity of these 
kinds of objects as a class, the temperatures are probably lower 
than 10"* K (see Díaz et al. 2007, and references therein). The low 
excitation of these objects estimated from their spectra makes any 
temperature sensitive line too weak to be observed. In fact, in most 
cases, the [O iii] A5007 A line, which is typically 100 times more in- 
tense than the auroral [O iii] A4363 A one, can barely be seen. In this 
way, we have also calculated the velocity dispersions assuming a 
kinetic temperature of 5000 K and we do not find great differences 
with the previous calculations derived using 10"* K. These differ- 
ences are between 0.02 and 3kms"' depending on the measured 
velocity dispersión since Oi is constant for a given temperature and 
subtracted in quadrature to a„. Obviously, the estimated a are al- 
ways higher for the lower assumed kinetic temperatures. Then, for 
consistency with Firpo et al. (2005) we use a kinetic temperature 
of 10"* K. Therefore, this effect owing to the kinetic temperature 
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uncertainty, would produce an underestimation in the derived true 
velocity dispersions. 

When fitting single Gaussian profiles to the emission lines we 
notice that, although it works well for most of the observed profiles, 
there is a residual present in the wings of several lines. This has 
been observed in GHiiRs before. Some authors (Melnick et al. 
1999; Hagele et al. 2007, 2009, 2010, among others) have shown 
that the profile wings are well fitted by single (broad) Gaussian of 
supersonic widths, whereas other authors (Chu & Kennicutt 1994; 
Relaño & Beckman 2005; Rozas et al. 2006) proposed that the 
profiles appear to be dominated by a few, well-resolved, expanding 
shells. Wherever possible, we have considered and applied both 
suggestions. 

In order to fit múltiple components to the observed profiles we 
make use of the ngaussfit task within iraf. ngaussfit performs 
iterative fitting of múltiple Gaussian profiles to spectral lines. The 
task allows us to select individual parameters to be fitted and delivers 
an estímate for the uncertainty of derived parameters. The task needs 
initial guesses for the function coefficients, and these are specified 
through previously generated tables. 

In those cases when more than one component is evident, we 
follow an iterative procedure in which we only allow the task to fit 
a limited subset of parameters at one time. This is done to constrain 
the universe of possible solutions to the fit by making use of the 
different Information available from each emission Une. Basically, 
after obtaining an initial guess for the second component parame- 
ters, we fit the amplitude, centre and width of each component at 
one time, leaving the parameters of the other one untouched. Once 
the central wavelength of all components is relatively well known, 
we set them as fixed parameters and allow the task to perform 
a final fit of amplitude and width of all components at the same 
time. 



The profile fittings show the presence of a residual in the 
emission-line wings for all cases, becoming more evident in the 
strong Ha emission line. Therefore, we evalúate the presence of 
either a broad component or shell features. For the former assump- 
tion, we introduce a broad component with an initial width which 
is three times the width of the narrow component(s). This guess 
comes from an overall view of the valúes found by Hagele et al. 
(2007) when performing a similar analysis. We fit the peak valúes 
for the narrow components and then their profile widths. The broad 
component parameters are then fitted before we start the iterative 
procedure again. In the latter scenario, a similar iterative procedure 
takes place, but in this case we add a blue and a red narrow com- 
ponent to account for the additional emission. The validity of the 
profile multiplicity and broadening is checked over the different 
emission lines available for each región, although the most reliable 
profile analysis comes from the strongest lines, such as Ha. 

Out of our sample of six GH iiRs we find that all of them show 
evidence of wing broadening. For these regions we attempt to fit an 
overall broad component or two narrow components symmetrically 
shifted in velocity with respect to the intense component. In the 
foUowing subsections we will discuss our findings grouping them 
by the outcome profiles. 

3.1 NGC 7479 1 

Rozas et al. (1999) undertook a complete study to date of the Hii 
regions in NGC 7479 galaxy, determining positions, angular sizes 
and absolute fluxes of over 1000 Hii regions and they constructed 
the luminosity function for the regions over the whole galaxy. We 
have cross-correlated these results with our sample and find that 
the three NGC 7479 regions selected in our work can be identified 
in the Rozas catalogue: NGC 7479 1 as RZH99 = 1, NGC 7479 II 



Table 2. Results of Gaussian profiles fitting to the observed emission lines in NGC 7479. Each emission line is identified by its ion laboratory wavelength and 
ion ñame in columns 1 and 2. According to the different fits performed on each line, column 3 identifies each 'narrow component' (A and B, where applicable), 
a broad component and/or blue and red wings (b wing and r wing, respectively). Radial velocities (Vr) and intrinsic velocity dispersions (a^a) together with 
their respective errors are expressed in km s^ ' . The intrinsic velocity dispersions are corrected for the instrumental and thermal widths. EMs are included and 
shown as a percentage of the component flux relative to the total EM of the región. 
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Figure 3. ngauss fits with three Gaussian components in the NGC 7479 1 emission-line profiles, two narrow principal components and a broad one. In order: 
Ha, [Nii] A,6584Á and [S ii] A,6731 Á. 



as RZH99 = 7 and NGC 7479 III as RZH99 = 74 (RZH99 is the 
number identification used by Rozas et al. 1999). 

We identify and fit the Gaussian profiles to the Ha, [N ii] A6584 A 
and [Sil] XX6111, 6731 Á Unes in NGC 7479 1 región with splot 
routines. The Gaussian fits reveal the presence of two distinctly 
separated kinematic components (labelled A and B) in all analysed 
profiles. Component A shows a much broader profile than compo- 
nent B, and both components show little spread in their individual 
radial velocities among different emission lines present in the spec- 
trum. The reliability of these valúes is confirmed when we improve 
the profile fitting using ngaussfit which in turn yields valúes for the 
profile width of each component. 

Even considering the presence of two distinct components with 
different radial velocities, the overall fit cannot account for a resid- 
ual emission in the profile wings. Following the procedures outlined 
in the previous section, we are able to fit a broad component, with 
a velocity dispersión of about 75 km s"' from the Ha emission line 



Figure 4. 2D Ha velocity iraage of NGC 7479 1 región with the [N ii] 
A6584Á (red) contours overlapped. Three sections are numbered in the 
2D image. 



Table 3. Gaussian fitting results for individual sections of the Ha emission- 
line profile in NGC 7479 1. Kinematical parameters [radial velocity (V¡) 
with errors (kms^ )], the intrinsic velocity dispersions ((Tint) corrected for 
the instrumental and thermal widths (kms^') and EM of each component 
as a fraction of the total EM of the sarapled section (per cent) data are 
shown. 
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and about 60 km s^' from the other ones. An attempt to fit twin nar- 
row profiles did not converge to a reliable solution. Table 2 shows 
parameters for the three components that fit the global profile. In- 
dividual Gaussian component fluxes are listed as emission measure 
(EM) relative to the total line flux following the work by Relaño & 
Beckman (2005). The overall Ha flux, uncorrected for reddening, 
is found to be 3.37 x 10^'"* erg s^' cm^^. In Fig. 3 we show the 
NGAUSS fitting done with three different Gaussian components in the 
emission lines which have enough signal to provide an accurate fit. 
The emission lines with a low signal-to-noise ratio for which the 
fits do not provide accurate results for the broad component are 
not listed in the table. The velocity dispersión for the ionized gas 
is derived by taking cr; = 5.1 kms^' as instrumental width, and a 
kinetic temperature of T ~ 10** K as mentioned above. 

Since part of the spatial Information along the slit is lost through 
the aperture extraction process we genérate 2D velocity images in 
Ha, [Nii] A6584Á and [Sii] A6717Á lines. Taking into account 
that the slit size is 1 x 4arcsec^, the data are binned by 2pixels 
in both directions and the seeing was 1 arcsec, we take aperture 



sizes of 0.5 arcsec ~ 1 pixel, therefore, dividing the slit into eight 
sections. We extract traces of the spectrum at eight representative 
cuts along the slit, preserving the same spatial resolution. This is 
done for each echelle order where a line of interest is present and 
each spectral section is rebinned to velocity space. In this reference 
frame, individual slices are stacked together using the imstack task, 
which allows us to genérate 2D velocity images. Fig. 4 shows the 
2D Ha velocity image with the [N ii] A6584 A contours overlapped. 
It can be readily seen that individual knots can only be identified 
when the kinematical Information is considered in this V^ versus 
slit section plañe. This suggests that 3D spectroscopy will be able 
to provide a powerful insight on the structure of this type of star- 
forming knots. 

We have performed a detailed analysis on the slit stack, by fit- 
ting Gaussian profiles at each slit section. The resulting parameters 
are listed in Table 3 and the graphics output can be seen in Fig. 5, 
where shown 2, 4, 5 and 7 Ha sections, and 4 and 5 [Nii] A.6584Á 
sections. Focusing on the slit edges first (sections 1 to 3 and 6 to 
8) we identify two distinct kinematical components at 2500 and 
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Figure 5. ngauss fits with three (two) Gaussian components in the NGC 7479 1 sections, two (one) narrow principal components and a broad one. In order: 
Ha (sections 2, 4, 5 and 7), [Nii] A.6584Á (sections 4 and 5). 
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2540 kms^' labelled A and B, respectively. Each of these relatively 
narrow features has a broader component associated with it. Al- 
though the narrow component can still be deblended in the central 
región, we are not able to disentangle the broad ones, which results 
in a single broad profile with a velocity similar to the average be- 
tween components A and B. Regarding the velocity dispersión of 
each component, we have listed them in Table 3 as a reference, but 
they should not be taken as a definite valué as the signal-to-noise 
ratio is variable and somewhat low for the slit edges. 

If we compare these results with the radial velocity map derived 
for the ionized gas in NGC 7479 (Laine 1996) we find, from visual 
inspection, that the expected velocity at the position of NGC 7479 1 
ranges between 2520 and 2535 km s^' cióse to an average valué 
of our components A and B. A similar comparison, using the Hi 



map by Laine & Gottesman (1998) which has a better spatial and 
spectral resolution, yields a valué between 2530 and 2540 km s^' 
much closer to the valué measured for component B. Although the 
latter map is sampling neutral rather than ionized hydrogen, the 
former might be affected by the complex structure we find within 
the región. At this stage we can only suggest that component A 
shows an odd kinematic behaviour and that a detailed study of 
NGC 7479 1 and its surroundings is needed to clarify this issue. 

3.2 NGC 7479 II 

A similar procedure is used to analyse the profiles of the emission 
Unes of NGC 7479 II. In this case, two components are identiñed 
with different radial velocities, but as the overall signal obtained 
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Figure 6. ngauss fits with four Gaussian components in the NGC 7479 II emission-line profiles, two narrow principal components and two low-intensity 
components. In order: Ha and [Nii] 6584 A. 




6785 6790 

Wavclcngth (A) 



Figure 7. ngauss fits with two Gaussian components in the NGC 7479 III emission-line profiles. In order from left to right and from top to bottom: Ha, [N ii] 
A6584Á, [Sil] A6717Áand [Sii] A.6731Á. 
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Figure 8. ngauss fits with two Gaussian components in the NGC 6070 1 emission-line profiles. du Pont data on the left-hand side and MIKE data on the 
right-hand side. In order frora top to bottom panels: [Oiii] >.5007 A, Ha, [Nii] X6584 A and [Su] A.6717A. ngauss fits with two Gaussian components in the 
NGC 6070 1 emission hne profiles. du Pont data on the left-hand side and MIKE data on the right-hand side for [S ii] A.673 1 A line. 
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Figure 8 - continued 



for this región is lower than for NGC 7479 1 it is not possible to 
disentangle them with 2D spectra. These two narrow components 
are found in the Gaussian profiles of the Ha, H/S, [Nii] AA6548, 
6584 Á and [S ii] AA,67 17, 673 1 Á lines. In this región, however, the 
wings of the brightest component profile could not be fitted with one 
single broad component. On the other hand, we are able to fit it with 
two wings only for the two lines with the strongest signal, namely 
Ha and [Nii] A.6584A (see Fig. 6). The final fitting parameters are 
shown in Table 2. The overall Ha flux, uncorrected for reddening, 
is found to be 1 .06 x 10"'"* erg s"' cm"^. The emission lines with a 
low signal-to-noise ratio for which the fits do not provide accurate 
results for the broad and the two wings components are not Usted 
in the table. 

We compare the velocity separation between the redshifted (and 
blueshifted) components with the average of the central features 
and the Ha luminosity of the región. The velocity separation of the 
red- and blueshifted components in this work is in agreement with 
the regression found by Relaño & Beckman (2005) for their Hii 
regions sample. 

We also analyse the radial velocity behaviour within the galaxy, 
and observe a similar situation to that of NGC 7479 1. There is no 
difference, within the errors, when we compare our average valué 
derived from the Ha emission line with the Ha velocities field map 
(Laine 1996). If, however, the comparison is done against the Hi 
rotation velocity map (Laine & Gottesman 1998), knot B seems to 
follow galactic rotation whereas knot A does not. 

3.3 NGC 7479 III 

The observed spectrum of this región shows, although weak, the 
[O iii] XX4959, 5007 A lines. This is probably owed to a difference 
in metallicity, NGC 7479 III being slightly less metallic than the 
brightest two regions discussed above, or a difference in the ioniza- 
tion structure. The latter might be due to a relatively small difference 
in the hardness of the ionization radiation field (and, what is equiva- 
lent, in the ionizing stellar population) responsible for the excitation 
observed in the spectra or a difference in the nebular geometry (see 
a detailed discussion in Hágele et al. 2006, 2008; Díaz et al. 2007; 
Pérez-Montero et al. 2010, and references therein). Unfortunately, 
as we said above, since no temperature sensitive lines can be ob- 
served, we are not able to derive their metallicities, and distinguish 
the origin of this small difference in the excitation present in the 
spectra of these regions of NGC 7479. 

We identify and fit the Gaussian profiles to the [O iii] A5007 A, 
Ha, [N II] X6584 Á and [S ii] AA.67 17, 673 1 Á lines in NGC 7479 III 



región with ngaussfit. In this región we have found one supersonic 
component, with no evidence of múltiple contributions but it is 
possible to detect flux excess on the wings of the emission lines. We 
therefore fit a broad Gaussian component together with the narrow 
component (see Fig. 7). This procedure works well for almost all 
detected lines. In the case of [Oiii] A,5007A, however, there are 
no prominent residuals when we fit only one Gaussian component 
which has a dispersión of about 20kms"' and we are therefore 
not able to fit a broad Gaussian component. For the other listed 
lines the derived radial velocities, the velocity dispersions (using 
a^ = 5.4 km s"' and T ~ 10"* K) and errors are listed in Table 2. 
The overall Ha flux, uncorrected for reddening, is found to be 
5.43 X 10"'^ erg s"' cm"^. It is worth mentioning that there seems 
to be some systematic offset between the derived profile centres of 
the narrow and broad profiles which is, in average, about 10 km s"' . 
This happens for all lines, although it becomes more evident for [S ii] 
X673 1 A. An undiscovered blueshifted narrow component could 
account for the observed behaviour, although better signal-to-noise 
ratio data are needed to confirm this hypothesis. 

The radial velocity average derived for the single component 
is in agreement with the expected radial velocity at the position 
of NGC7479III, showing only a slight offset of (lOkms"'), wefl 
within the uncertainties. 

3.4 NGC 6070 1, II and IV 

All regions observed in NGC 6070 show one narrow component 
together with the underlying broad component. We identify and fit 
the Gaussian profiles to the H^, [O iii] A5007 Á, Ha, [N ii] A.6584 Á 
and [Su] AA6717, 6731 Á lines in almost all regions (see Figs 8-10). 
Table 4 shows the derived kinematical parameters (using a^ = 5.3, 
5.6 and 5.3 km s"' for regions I, II and IV, respectively; T is set to 
~ 10"* K for all regions). 

For these three regions, we are also able to compare the du Pont 
data with higher resolution spectra data obtained with the double 
echelle spectrograph MIKE in 2004 July. The agreement in the 
kinematical parameters can be seen in Table 4, which shows the 
reliability of the data here presented. The overall Ha flux, uncor- 
rected for reddening, for du Pont data are 1.68 X 10"'^, 6.26 X 10"'"* 



tively, and for MIKE data are 1.74 x 10"", 1.29 x 10"'^ and 
1.88 X 10"" erg s"' cm"^. Only a slight bias towards lower veloc- 
ity dispersión can be seen for NGC 6070 1, which can be probably 
attributed to differences in the centring of the slit during acquisition 
of the target. There seems to be some systematic offset between 
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Figure 9. ngauss fits with two Gaussian components in the NGC 6070 II emission-line profiles. du Pont data on the left-hand side and MIKE data on the 
right-hand side. In order from top to boltom panels: Hfi, Ha, [Nn] A.6584A and [Sii] A6731 A. 
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Figure 10. ngauss fits with two Gaussian components in the NGC 6070 IV emission-line profiles. du Pont data on the left-hand side and MIKE data on the 
right-hand side. In order from top to boltom panels: [O iii] A.5007 Á, Hü', [N ii] A.6584 Á and [S ii] A.6717 Á. 
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Table 4. Results of Gaussian profiles fitting to the observed emission lines in NGC 6070. Each emission line is identified by its ion laboratory wavelength and 
ion ñame in columns 1 and 2. According to the different fits performed on each line, column 3 identifies each 'narrow component' (A and B, where applicable), 
a broad component and/or blue and red wings (b wing and r wing, respectively). Radial velocities (V,) and intrinsic velocity dispersions ((Tí„i) together with 
their respective errors are expressed in km s^ ' . The intrinsic velocity dispersions are corrected for the instrumental and thermal widths. EMs are included and 
shown as a percentage of the component flux relative to the total EM of the región. Upper panel: du Pont data; lower panel: MIKE data. 
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the derived profile centres of the narrow and broad profiles. The- 
ses valúes are, in average, about — 251cms^' for du Pont data and 
- 10 km S-' for MIKE data in the case of NGC 6070 II, and about 
10 and 5 km S-' for NGC 6070 IV. 

Regarding radial velocities of the three observed regions in 
NGC 6070, there is no detailed velocity map in the literature. We 
are able to compare them with rotation curves derived from long 
slit observations made by IVIárquez et al. (2002). Although the ori- 
entation of the slits does not match exactly the location of all our 
regions, the velocities are in fair agreement considering the angular 
distance of the H ii regions to the centre of NGC 6070. 

3.5 Relation between Ha luminosities and velocity dispersión 

We have also analysed the impact that the presence of múltiple 
components has on the location of H ii regions in the log(L)- log(cr ) 
plañe. To evalúate this we have plotted in Fig. 1 1 the luminosities 
and velocity dispersions derived directly from our spectrophoto- 
metric data. Luminosities were derived from the fluxes measured 
directly from the component fitting to our echelle spectra (uncor- 
rected by reddening), and using distances as published by Rozas 
et al. (1999) for NGC 7479 and by García-Gómez et al. (2002) 
for NGC 6070. We are aware that slit spectra might not provide the 
most realistic sample of the overall flux of the ionized nébula, but the 



accuracy of acquisition and guiding is not good enough as to allow us 
to perform a proper identification with published luminosity valúes 
in fields relatively 'crowded' with Hii regions by Feinstein (1997) 
and Rozas et al. (2006). Velocity dispersión valúes and its errors are 
taken directly from Tables 2 and 4. The luminosity errors have been 
estimated taking into account the errors in the amplitude (A) and 
the FWHIvl in the component fitting (F = 1.0645 x A x FWHM), 
and the distance errors. As a reference valué we have plotted a few 
giant H II regions from Bosch et al. (2002) together with their linear 
fit to their 'young' GH ii regions. 

Inspection of Fig. 1 1 reveáis that the distribution of points in the 
log(L)-log((Ti„t) plañe are strongly dependent on the components 
derived from the profile fitting. Individual components, labelled A 
and B (when applicable) as in their respective tables, and global 
valúes are those derived from a single Gaussian fitting to the line 
profile. As expected, individual components have smaller fluxes and 
velocity dispersions than the global profile and points are therefore 
shifted in the diagram. The incidence of this analysis has, however, 
a different outcome for regions in NGC 6070 and NGC 7479. In 
NGC 6070 the valúes derived from the global profile seem to lie, 
within observational errors, in the expected location if they follow 
the relation expected for virialized systems. When the individual 
components are plotted in the same diagram, only NGC 6070 II lies 
cióse to the regression, and NGC 6070 1 and IV seem to be too 
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Figure 11. log(L)-log((T) relation for our Hii regions. Luminosities and 
velocity dispersions are derived from our spectrophotometric data. The plot 
includes results from individual components labelled A and B (where ap- 
plicable) together with velocity dispersions for the overall profile (labelled 
as 'global') which were derived from a single Gaussian fitting to the Une 
proñle (colour solid error bars). The narrow components in NGC 6070 were 
labelled A, for which we only plot the MIKE measured data. A few giant 
H II regions from Bosch et al. (2002) (blue dashed error bars) together with 
their linear ñt to their 'young' giant H ii regions are plotted as a reference 
valué. The luminosities are not corrected for extinction. 

luminous for their measured velocity dispersión. On the other hand, 
two out of three regions in NGC 7479 he weü away this regression 
if parameters are derived from a single Gaussian, but four out of five 
individual components show a very good agreement when plotted 
in the same diagram. 

The interpretation of this behaviour is not straightforward. Def- 
initely, the presence of more than one Gaussian component rules 
the final position of the Hii regions in the log(L)-log(o') plañe. 
Our high-resolution spectroscopic data have allowed us to dis- 
entangle múltiple kinematic components, which reflects in a big 
improvement in the x-axis of the plot. However limitations in spa- 
tial resolution prevent us from identifying múltiple components, 
if present, which do not become evident in the spectral signature. 
This has been a problem in previous attempts to derive relations in 
the log(L)-log((7) as discussed by Bosch et al. (2002) when they 
were able to obtain photometry for individual knots within giant 
Hii regions in MlOl. 



4 SUMMARY AND CONCLUSIONS 

From new high-resolution spectra of the H ii regions NGC 7479 1, 
NGC 7479 11, NGC 7479111, NGC 6070 1, NGC 6070 II and 
NGC 6070 IV obtained at the 100-inch du Pont Telescope, LCO, 
we have confirmed the giant nature of all these regions. We have 
also found that all of them show evidence of wing broadening evi- 
dent mainly in the Ha line and confirmed in other emission lines. In 
five of them we are able to fit a broad component which explains the 
integral profile wings. Only in NGC 7479 II can we fit the integrated 
profile wings with two narrow components symmetrically shifted in 
velocity with respect to the average of the two intense components. 
NGC 7479 1 and NGC 7479 II reveal the presence of two sepa- 
rated kinematic components with relatively narrow profiles in all 
analysed emission lines. For NGC 7479 1 región we can also gen- 
érate 2D velocity images in Ha, [Nii] A,6584Á and [Sii] A6717Á 



lines. These images show that the two narrow components are spa- 
tially resolved and the kinematical information helps to split the 
knot in the velocity-distance plañe. There is also evidence of the 
presence of a broad component in all the narrower extracted sec- 
tions, at least in the Ha emission line. For NGC 7479 II región we 
can also fit two relatively narrow Gaussian components. In both 
regions the radial velocity average valué of our components A and 
B derived from the Ha emission line matches the global velocity in 
the galaxy when we project our valué over the isovelocity contours 
in Ha velocity field map. A similar comparison, using the H i map 
yields a valué much closer to the valué measured for component B. 
Then, we could only suggest that component A shows an odd kine- 
matic behaviour. New observations with better spatial and spectral 
resolution are needed to clarify this point. 

In the rest of the studied H ii regions, we find one narrow Gaussian 
component together with an underlying broad component, with 
no evidence of múltiple narrow contributions. The radial velocity 
average derived for the narrower component corresponds to the 
one expected from galaxy rotation curves, within the observational 
errors, in all of these cases. 

In most cases, the velocity dispersions of the broad components 
for these H ii regions derived by optimal Gaussian fit are in agree- 
ment with Hágele et al. (2009, 2010) for circumnuclear star-forming 
regions with cri„i in H/S line ~34 to 65 km s^', and those derived by 
Chu & Kennicutt (1994) and Melnick et al. (1999) for 30 Doradus 
nébula with cri„, in Ha line ~45 km s^' . 

The estimated offsets between the narrow and broad components 
in NGC 7479 III, NGC 6070 II and NGC 6070 IV are in complete 
agreement with the valúes found in NGC 2903 and NGC 3310 by 
Hágele et al. (2009, 2010), respectively, which are between —25 
and 35kms"'. 

The distribution of the regions in the log(L)-log(cr) plañe are 
strongly dependent on the components derived from the profile 
fitting. Individual components have smaller fluxes and velocity dis- 
persions than the global profile and points are therefore shifted in 
the diagram. In NGC 6070 the valúes derived from the global pro- 
file seem to lie, within observational errors, in the expected location 
if they foUow the relation expected for virialized systems. When 
the individual components are plotted in the same diagram, only 
NGC 6070 II lies cióse to the regression, and NGC 6070 1 and IV 
seem to be too luminous for their measured velocity dispersión. On 
the other hand, two out of three regions in NGC 7479 lie well away 
this regression if parameters are derived from a single Gaussian, but 
four out of five individual components show a very good agreement 
when plotted in the same diagram. Definitely, the presence of more 
than one Gaussian component rules the final position of the Hii 
regions in the log(L)-log((T) plañe. 
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3.1. Resumen de la publicación 

En este Capítulo, haciendo uso de espectroscopia échelle obtenida en el Observatorio Las 
Campanas, se presenta un estudio detallado de la cinemática interna del material nebular 
en múltiples brotes de formación estelar de la galaxia enana campacta azul Haro 15. Un 
análisis detallado de los complejos perfiles de líneas de emisión muestran la presencia de una 
componente ancha subyacente en casi todos los brotes, y la región de formación estelar más 
brillante muestra indudables signos de la presencia de dos componentes cinemáticas angostas 
diferentes. Estudiando la información que este análisis proporciona sobre el movíento de los 
brotes individuales en el potencial de la galaxia Haro 15, se confirma que siguen la rotación 
galáctica, aunque las componentes angostas del brote A tienen velocidades relativas demasiado 
grandes para ser explicada por rotación galáctica. Por último, se examina la relación entre la 
dispersión de velocidades y la luminosidad, encontrando que casi todos los brotes siguen la 
relación de sistemas virializados. Esto es válido para las componentes angostas identificadas 
en los perfiles complejos, en cambio los perfiles ajustados con una componente Gaussiana 
simple muestran una pendiente más plana. 

33 



CAPITULO 3. ESPECTROSCOPIA DE ALTA RESOLUCIÓN DE LA GALAXIA BCD 
34 H ARO 15: CINEMÁTICA INTERNA 

3.2. Fundamentos básicos 

Durante años, varios estudios en diferentes frecuencias han sido llevados a cabo sobre la 
galaxia enana compacta azul Haro 15. 

Desde el punto de vista morfológico, la imagen Ha de Cairos et al. (2001) muestra que la 
galaxia presenta una compleja morfología resuelta en un gran número de brotes dispersos por 
toda la galaxia (ver Figura 1 de la publicación). Haro 15 ha sido clasificada como una galaxia 
peculiar (R)SBO (de Vaucouleurs et al., 1991) y, de hecho, en las imágenes profundas tomadas 
con el telescopio de 2.2m CAHA se puede apreciar una morfología en espiral (López-Sánchez 
& Esteban, 2008). López-Sanchez et al., (2010) encontraron una cinemática y abundancias 
químicas diferentes entre las regiones A y B lo cual podría estar indicando que la galaxia 
probablemente esté experimentando una fusión menor. 

En la publicación se hace mención a que Haro 15 cumple con los criterios de una galaxia 
luminosa compacta azul (LCBGs, sigla en inglés de Luminous Compact Blue galaxies), y 
es que el nombre de estas galaxias fue asignado en Jangren et al., (2004) a las galaxias 
luminosas (Mb<-17.5), compactas (;^b<21.0 mag arcsec^) y azules (B-V<0.6) las cuales 
están sometidas a un importante estallido estelar. Esta definición intenta incluir las galaxias 
BCGs locales más luminosas, como así también a una variada familia de brotes de formación 
estelar a redshifts intermedios. Generalmente, estos objetos son poblaciones aisladas de alto 
brillo superficial pudiendo ser observados a altos redshifts y tienen propiedades muy similares 
a las regiones de formación estelar de alta luminosidad en las galaxias H ii (Hoyos et al., 2004). 

Como se vio en el Capítulo anterior, la espectroscopia échelle proporciona un medio para 
buscar componentes cinemáticamente diferentes en las líneas de emisión del gas ionizado, ya 
que llega a la resolución espectral necesaria para resolver la presencia de estructuras en el 
perfil de líneas de emisión, generalmente enmascarados por el gran ancho supersónico y la 
baja resolución espectral instrumental. 

En el Capítulo 2 (Firpo et al., 2010) se observó que al hacer un ajuste Gaussiano simple 
en el perfil de la línea, existe un residuo presente en las alas del perfil. Basado en la literatura, 
y siempre que fuera posible, se evaluó la posible presencia de una componente ancha o dos 
componentes simétricas a ambos lados de la componente angosta en el ajuste de los anchos 
de las emisiones observadas. En el presente trabajo, también se ha encontrado que todas las 
regiones de la galaxia Haro 15 muestran este ensanchamiento siendo más evidente en el perfil 
de la línea Ha y luego se confirma en el resto de las líneas de emisión. Por otro lado, se analiza 
la relación entre la dispersión de velocidades en la línea Ha y su luminosidad estudiando cómo 
repercute la existencia de diferentes componentes cinemáticas en el plano L(Ha) vs. a. 

En este Capítulo se presenta espectroscopia échelle obtenida con el telescopio de 100" du 
Pont, del Observatorio Las Campanas (LCO), en cinco brotes de formación estelar de Haro 15 
con una resolución en velocidad de 12 kms~^. Esta resolución permite la identificación de las 
diferentes componentes cinemáticas del gas y la medición de sus correspondientes dispersiones 
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de velocidades. La identificación de los brotes sigue la nomenclatura utilizada por Cairos et 
al., (2001). Dos brotes más fueron observados en este trabajo, a los cuales se los ha denominado 
E y F, continuando la notación de Cairos. 

Este estudio forma parte de un proyecto en el que, con la obtención de datos de alta 
resolución espectral échelle, se determina la naturaleza de las Regiones Hil Gigantes visibles 
desde el hemisferio sur y se analizan, además, las condiciones físicas del gas ionizado de las 
regiones en esta galaxia enana compacta azul. En la publicación a continuación se presentan 
las observaciones y la reducción de los datos en la Sección Observations and Reductions; 
en la Sección Results and Discussion se presenta el análisis de los perfiles de la línea 
de emisión y se analizan los resultados de cada brote individual. Además se estudian las 
velocidades radiales y se analiza el impacto que tiene la presencia de varias componentes 
sobre la ubicación de cada región Hii en el plano log(L) — log{a); por último se incluye 
resumen y conclusiones del trabajo en la Sección nombrada como Summary. 

3.3. Herramientas utilizadas 

Los espectros alta resolución de Haro 15 se obtuvieron utilizando el espectrógrafo échelle 
del telescopio de 100" du Pont, del Observatorio Las Campanas (LCO), Chile, en julio del 
2006. El rango espectral cubierto por las observaciones fue de A3400 a AlOOOO A. Las condi- 
ciones de observación fueron buenas, con un seeing de 1 segundo de arco. Se aplicó un binning 
2x2 al CCD el cual tiene una resolución espacial de 0.259"px^^. La ranura utilizada tenía un 
ancho de 1" y una longitud de 4". Considerando que las exposiciones son largas, la resolución 
espacial esta limitada no solo por seeing sino también por el guiado, y teniendo en cuenta 
la distancia a Haro 15, la escala espacial de los espectros nos dice que la menor estructura 
que se puede ver es de al menos 0.43 kpc. La resolución espectral alcanzada en los espectros 
échelle es de R~25000: AA=0.25A a los A 6000A, medida a partir del ancho a potencia mitad 
(FWHM) de las líneas de la lámpara de comparación de Thorio-Argón. Esto se traduce en 
una resolución de ~ 12 kms~^. 

El tiempo de exposición para cada región fue 1800 segundos. También se observó para 
la calibración en fiujo una estrella estándar espectrofotométrica, FeigellO, con un tiempo 
de exposición de 1200 segundos. Además, durante todas las noches se tomaron espectros de 
lámpara de comparación Th-Ar, Milky Fíats (flats de cielo obtenidos con un difusor durante 
la tarde), y bias. 

El análisis de los datos se realizó con el software IRAF^. Luego de la resta de bias, 
correcciones por ñat fields con los Milky Flats, las imágenes bidimensionales fueron corrigidas 
por rayos cósmicos con la tarea cosmicrays de IRAF, la cual detecta y elimina los rayos 
cósmicos con un algoritmo de cociente de flujos. Los datos corregidos se redujeron con las 



^Image Reduction and Analysis Facility, distribuido por NOAO, operado por AURA, Inc., bajo contrato 
con NSF. 
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rutinas de IRAF siguiendo el mismo procedimiento descrito en Firpo et al., (2005). 

3.4. Metodología en la cinemática 

En el espectro échelle calibrado, se cortó el rango de longitud de onda en una dada línea 
de emisión y este ha sido transformado al plano de velocidad aplicando la corrección Doppler. 
Midiendo la velocidad central (longitud de onda) y el ancho en varias líneas de emisión se 
han determinado las velocidades radiales y la dispersión de velocidades del gas ionizado en 
las diferentes regiones de formación estelar. La dispersión de velocidades intrínseca (amt), 
fue corregida por las contribuciones instrumental y térmica de cada línea de emisión. Los 
pequeños cambios en las temperaturas electrónicas, de unos pocos cientos de grados, no mo- 
difican sustancialmente esta corrección, como ya se ha discutido en el Capítulo anterior (Firpo 
et al., 2010). En el Capítulo 4 se presentará un tratamiento detallado en la determinación de 
temperaturas y densidades electrónicas de cada uno de estos objetos. 

Como fue comentado en el Capítulo anterior, en el presente Capítulo los brotes de for- 
mación estelar de la galaxia BCD, HarolS, muestran evidencias de residuo presente en las 
alas de las emisiones, el cual es observado en todas las líneas de emisión. Haciendo uso del 
proceso iterativo, explicado en Firpo et al., (2010), para el ajuste de múltiples componentes 
Gaussianas en el perfil de línea, se evaluó la presencia de una componente ancha y la presencia 
de más de una componente angosta. En este caso, se ha podido ajusfar una variada gama de 
componentes, explicando así las alas del perfil integrado en cada una de las regiones. 

Todas las regiones, a excepción del brote F, muestran múltiples componentes supersónicas. 
En el brote F el ensanchamiento del perfil integrado solamente pudo ser detectado en la línea 
Ha. En esta línea de emisión se pudo ajusfar una componente ancha supersónica y una angosta 
subsónica de 8 kms~^, dispersión de velocidades típica de regiones Hll clásicas, aunque la 
presencia de una componente ancha sea casi exclusiva de las Regiones H ii Gigantes. 

En la publicación se discuten los resultados agrupándolos según los resultados de los 
perfiles y por brote. Este trabajo ha sido aceptado el día V° de marzo de 2011 para su 
publicación en la revista internacional con referato Monthly Notices of the Royal Astronomical 
Society. 
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ABSTRACT 

Using echelle spectroscopy, obtained at Las Campanas Observatory, we present a detailed study 
of the internal kinematics of the nebular material in múltiple knots of the blue compact dwarf 
galaxy Haro 15. A detailed analysis of the complex emission-line profiles shows the presence 
of an underlying broad component in almost all knots, and the brightest star-forming región 
shows unmistakable signs for the presence of two distinct narrow kinematical components. 
We also study the Information that our analysis provides regarding the motion of the individual 
knots in the Haro 15 galaxy potential, confirming that they follow galactic rotation. Finally, we 
examine the relation between their velocity dispersión and luminosity, finding that almost all 
knots follow the relation for virialized systems. This holds for the strong narrow components 
identified in complex fits and for single profile fits, although the latter show a flatter slope. 
In agreement with previous ñndings, in this paper we show that the existence of múltiple 
kinematical components among massive starbursts cannot be overlooked, as it has a noticeable 
effect on any subsequent analysis that relies on basic parameters. 

Key words: Hii regions - galaxies: individual: Haro 15 - galaxies: starburst. 



1 INTRODUCTION 

Already in the original list of blue galaxies with emission Unes 
by Haro (1956) who described this object as a 'minute cometary 
nébula', Haro 15 was later included in severa! compilations of blue 
compact galaxies (BCGs). Several studies at all frequencies have 
been carried out of this object. In particular, optical spectroscopy 
was analysed by Hunter & Gallagher (1985), Mazzarella, Bothun 
& Boroson (1991), Kong et al. (2002) and Shi et al. (2005), among 
others. Owing to the detection of the Heii X4686 emission line 
by Kovo & Contini (1999), Schaerer, Contini & Pindao (1999) 
have classified this galaxy as a Wolf-Rayet (WR) galaxy. López- 
Sánchez & Esteban (2010) also found a blue WR bump supporting 
the WR nature of Haro 15. The Ha image shown in Cairos et al. 
(2001a) shows a knotty morphology with the starburst región re- 
solved into a large number of clumps which appear scattered over 
the entire galaxy. The integrated Ha luminosity yields a star forma- 
tion rate (SFR) of 3.3 M© yr"' (López-Sánchez & Esteban 2008). 
López-Sánchez (2010) analysed the SFR in Haro 15 using multi- 
wavelength data. The derived valúes using Ha, far-infrared (FIR) 
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and radio data are very similar, but the valué derived using the far- 
ultraviolet (FUV) luminosity is twice the others, confirming that a 
young stellar population is dominating the light of the galaxy and 
suggesting that the starburst phenomenon in Haro 15 started some 
time ago (at least 100 Myr ago). From the morphological point of 
view, Haro 15 has been classified as an (R)SBO peculiar galaxy 
by de Vaucouleurs et al. (1991) and in fact, in the deep images 
taken with the 2.2-m CAHA telescope shown by López-Sánchez 
& Esteban (2008), the spiral morphology of the galaxy can eas- 
ily be appreciated. Two high surface brightness concentrations can 
be distinguished, named A and C by the authors. Both show blue 
colours and high FUV emission, indicative of recent and ongoing 
star formation activity, further supported by the presence of the WR 
features. 

Long slit optical spectroscopy of Haro 15 has been presented in 
López-Sánchez & Esteban (2009) covering four main regions in the 
galaxy: the centre (named C by the authors) and three bright regions 
located ESE (named A), WNW (named D) and NE (named B). The 
spectrum of the central región shows strong nebular emission lines 
and prominent stellar absorption wings in the H i Balmer lines, evi- 
dence of an underlying stellar population which the authors estímate 
to be around 500 Myr oíd. Only emission lines are seen in regions 
A and B and the spectrum of región D is too noisy to be analysed. 
These authors derived an oxygen abundance of 12 + log(0/H) = 
8.37 ± 0.10 and 8.10 ± 0.06 for regions A and B, respectively. 
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Based on these data the authors conclude that all the observed knots 
can be classified as typical H ii regions. 

Haro 15 has an absolute magnitude Mj¡ = —20.69, a surface 
brightness fin = 18.56magarcsec"^ and a colour B — V = 0.33 
(Cairos et al. 2001a). At a distance of 86.6 Mpc (de Vaucouleurs 
et al. 1991), Haro 15 meets the criteria for a luminous compact blue 
galaxy (LCBG; Hoyos et al. 2004) despite the fact that its bright- 
ness distribution shows an exponential profile with a scalelength of 
1.37 kpc. The nature of this type of galaxy is not yet clear since 
they probably constitute a mixed population of starburst galaxies. 
Some authors suggest that LCBGs might represent the final out- 
come of a merger between a dwarf elliptical and a gas-rich dwarf 
galaxy or Hi cloud (Ósthn 1998; Cumming et al. 2008). The in- 
teractions taking place during the merging process would act as 
the starburst trigger (see e.g. López-Sánchez, Esteban & Rodríguez 
2004; López-Sánchez, Esteban & García-Rojas 2006). The images 
of Haro 15 showing two separated nuclei surrounded by a more 
regular, roughly elliptical envelope with twisted isophotes, together 
with different features and faint extensions as reported by Cairos 
et al. (2001a), give support to this picture. Other galaxy properties 
like different kinematics and chemical abundances between regions 
A and B may also indícate that Haro 15 is probably experiencing a 
minor merger (López-Sánchez 2010). On the other hand, a reliable 
determination of the mass of this kind of object is needed in order to 
decide on the evolutionary path they may follow: becoming dwarf 
spheroidal galaxies (Hoyos et al. 2007) or the bulges of small spirals 
(Hammer et al. 2001). However, interactions and mergers, as well 
as feedback processes resulting from intense star formation, might 
be reflected in a peculiar gas kinematics that prevent the derivation 
of their dynamical masses from the gas velocity widths. 

Echelle spectroscopy provides a means to look for kinematically 
different components in the emission Unes of the ionized gas since 
it reaches the spectral resolution needed to resolve the presence of 
structures within the emission-line profile, usually masked by its 
large supersonic width. 

In Firpo et al. (2010) we observed a residual presence in the 
wings of several lines when fitting single Gaussian profiles to the 
emission lines observed in the giant H ii regions within the galaxies 
NGC 7479 and 6070. Based on the several studies that have been 
proposed in the literature to interpret the existence of the broad su- 
personic component measured in the emission-line profile of giant 
H II regions, and whenever possible, we have evaluated the possible 
presence of a broad component (Muñoz-Tuñón et al. 1996; Mel- 
nick, Tenorio-Tagle & Terlevich 1999; Hágele et al. 2007; Hagele 
et al. 2009, 2010, among others) or two symmetric low-intensity 
components in the fit to the observed emission-line profile widths 
(Chu & Kennicutt 1994; Relaño & Beckman 2005; Rozas et al. 
2006). In the present work, we have also found that all Haro 15 
regions show evidence of wing broadening evident mainly in the 
Ha line and confirmed in other emission lines. Another common 
feature of giant H ii regions is the supersonic width of their emission- 
line profiles, although the origin of these velocities is not yet clear 
(Terlevich & Melnick 1981; Hippelein 1986; Roy, Arsenault & 
Joncas 1986). A relation between Ha luminosity and velocity dis- 
persión of the form L(Ha) a a'^ is expected under the assumption 
of a gravitational origin of the ionized gas dynamics. Therefore, a 
study of the implications that the existence of several kinematically 
distinct components has on the L(Ha) versus a relation is of great 
valué. 

In this paper we present echelle spectroscopy obtained with the 
100-inch du Pont telescope at Las Campanas Observatory (LCO) 
of five different positions across Haro 15 which provide a velocity 



resolution of about 1 2 km s" ' allowing the Identification of different 
kinematical components of the gas and the measurement of their 
corresponding velocity dispersions. This study is part of a project 
to obtain high spectral resolution echelle data to determine the 
nature of giant H ii regions visible from the southern hemisphere 
and analyse the physical conditions of the ionized gas of these 
regions and blue compact dwarf galaxies. In Section 2 we present the 
observations and the data reduction. Section 3 presents the analysis 
of the emission-line profiles and discusses the results. Finally, the 
summary and conclusions of this work are given in Section 4. 

2 OBSERVATIONS AND DATA REDUCTION 

High-resolution spectroscopy of Haro 15 was obtained using an 
echelle spectrograph attached to the 100-inch du Pont telescope, 
LCO, on 2006 July 19 and 20. The spectral range covered by the 
observations was from 3400 to 10 000 A. Observing conditions were 
good, with an average of 1 arcsec seeing and photometric nights. A 
2x2 binning was applied to the CCD in order to minimize the read- 
out contribution to the final spectrum noise. With 1 arcsec effective 
slit width and 4 arcsec slit length, the spectral resolution achieved in 
our du Pont echelle data was i? ~ 25 000: AA. = 0.25 Á at X6000 Á, 
as measured from the fuU width at half-maximum (FWHM) of the 
Th-Ar comparison lines taken for wavelength calibration purposes. 
This transíales into a velocity resolution of ~ 12 km s"' . Consider- 
ing the distance to Haro 15 and that the spatial resolution is limited 
by seeing, the smallest structure that can be resolved is 0.43 kpc 
(~1 arcsec ~ 1.93 pixel). 

Five different regions in Haro 15 were observed as shown in 
Fig. 1. The spectra were obtained as single exposures of 1800 s. 








10 arsec = 4.2 kpc 

-* >- 
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Figure 1. Finding chart for the giant Hii región candidates observed in 
Haro 15, identified by circles and labelled foUowing Cairos et al. (2001b). 
The Wide Field Planetary Camera 2 Ha iraage was obtained from the 
Multimission Archive at the Space Telescope Science Institute (MAST). 
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Table 1. Journal of observations. The identifications in the first column 
correspond to the nomenclature used in this paper, foüowing and extending 
that in Cairos et al. (2001b); columns 2 and 3 list the observation date 
and the airmass, respectively; columns 4 and 5 show the distance to the 
galactic centre in arcsec and kpc, respectively, and finally, column 5 lists 
the identiñcation of the regions given in López-Sánchez & Esteban (2008) 
referred to as LS notation. The exposure time for each región was 1 800 s. 



Knots 



Date 



sec z Dist. (arcsec) Dist. (kpc) LS 



Haro 15 A 2006Julyl9 1.2 O O C 

HarolSB 2006Julyl9 1.1 11 4.62 A 

Harol5C 2006Julyl9 1.1 20 8.4 B 

HarolSE 2006 July 20 1.2 16 6.72 

Harol5F 2006July20 1.1 17 7.14 



CALSPEC spectrophotometric standard star Feige 110 (Bohlin, 
Dickinson & Calzetti 2001) was also observed for flux calibra- 
tion purposes with an exposure time of 1200 s. In addition, Th-Ar 
comparison spectra, milky flats (sky flats obtained with a diffuser, 
during the afternoon) and bias frames were taken every night. A 
joumal of observations is shown in Table 1. 

The data analysis was carried out using the iraf' software. After 
bias subtraction and flat-field corrections by means of milky flats 
the bidimensional images were corrected for cosmic rays using the 
task cosMiCRAYS which detects and removes cosmic rays using a flux 
ratio algorithm. The corrected data were reduced by iraf routines 
following procedures similar to those described in Firpo, Bosch & 
Morrell (2005). 

At the end of the process, we compared the red end of the 
wavelength-calibrated spectra with the night-sky spectrum by Oster- 
brock et al. ( 1 996). This turned out to be a very reliable confirmation 
of the goodness of the wavelength solution, and we checked that 
differences between our wavelengths and the sky line wavelengths 
were below 0.04 Á. 

Flux calibration was achieved by observing the CALSPEC spec- 
trophotometric standard star Feige 110 (Bohlin et al. 2001) whose 
flux was tabulated every 2 A. Despite its relatively low brightness 
(y = 11.83), Feige lio is ideal for calibrating high-resolution 
echelle spectra. The amount of defined intervals within an echelle 
order ranged from four to 12, depending on the quality of the spec- 
trum. The flux calibrations were performed as described in Firpo 
et al. (2005). 

3 RESULTS AND DISCUSSION 

3.1 Line proflle analysis 

We identified the hydrogen recombination Unes, such as Ha and H/í, 
and colUsionally excited Unes, such as [Nii] A.A6548, 6584, [Sii] 
AA.67 1 7, 673 1 , present in the spectra, making use of the known red- 
shift, z = 0.021 371, for Haro 15 (de Vaucouleurs et al. 1991). The 
strong Unes were used to analyse the structure of velocity profiles as 
they allow us to verify the existence of more than one component as 
described in Firpo et al. (2010). The adopted laboratory wavelengths 
were taken from the work of García-Rojas et al. (2005). 

From the echelle calibrated spectrum, we cut the wavelength 
range where a given emission line is and we transformed from wave- 
length to velocity plañe using the Doppler correction. By measuring 



' Image Reduction and Analysis Facility (iraf), distributed by NOAO, op- 
erated by AURA, Inc., under agreement with the NSF. 
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the central velocity (wavelength) and width of several emission Unes 
we determined the radial velocities and velocity dispersions of the 
ionized gas in the different star-forming regions of Haro 15. The 
radial velocity and the intrinsic velocity dispersión (ffin,), corrected 
for the instrumental and thermal contributions of each emission 
line, were also derived. For these observations we considered a^ = 
5.2 km s^' as the instrumental width. The thermal contribution was 
derived from Boltzmann's equation (cti = 2kT/m.^), where k is 
Boltzmann's constant, T the kinetic temperature (T ~ 10"* K) and 
rria the atomic mass of the corresponding element. Although a de- 
tailed determination of the electrón density and temperatures, and 
chemical ionic and total abundances for each región will be pre- 
sented in Hagele et al. (in preparation), small changes of a few 
hundred degrees do not noticeably modify this correction, as previ- 
ously discussed in Firpo et al. (2010). The fluxes were derived from 
the amplitude (A) and the FWHM of the Gaussian proflle obtained 
in the component fitting (F = 1.0645 x A x FWHM), and the 
corresponding errors were estimated taking into account the errors 
in these two parameters. 

As already reported by Firpo and collaborators, in the present 
work we have also found that all Haro 15 emission knots show 
evidence of wing broadening, which is always found in the Ha 
line profiles and is usually also observed in the profiles of the 
bright emission Unes. Making use of the iterative fitting of múltiple 
Gaussian profiles, we evaluated the presence of a broad component 
and more than one narrow component present in the emission-line 
profile. In this case, we fitted a broad component, explaining the 
integral profile wings for all regions. 

In the following subsections we will discuss our findings for each 
knot resulting from the profile fits. 



3.1.1 Haro 15 A 

In our high-resolution spectra Haro 15 A shows a complex struc- 
ture which, although evident in radial velocity space, could not be 
spatially resolved. 

We identified and fitted Gaussian profiles to the H/!, [O iii] A5007, 
[Nii] A6548, Ha, [N ii] A6584 and [S ii] A.6717 Unes in this knot with 
NGAUSSFIT routines. The Gaussian fits of these emission line profiles 
are shown in Figs 2 and 3. The Gaussian fits in the profiles reveal 
the presence of two distinctly separated kinematical components 
labelled narrow 1 and 2 (ni and n2). Component ni shows a pro- 
file slightly broader than component n2 (cTim ~ 28 and 24kms"', 
respectively), and both components show a relatively large spread 
in individual radial velocities among the different emission Unes 
present in the spectrum. The reliability of these valúes is confirmed 
when we improve the profile fitting using ngaussfit which in tum 
yields valúes for the profile width of each component. It is worth 
noting that, although the profile fitting to the [Oiii] Unes shows 
similar overall results, the n2 components show the broader profile, 
opposite to what is found for the other emission Unes. This could be 
related to a different kinematic behaviour of the highly ionized gas, 
although this needs to be confirmed for other high excitation Unes. 

Always considering the presence of two distinct components 
with different radial velocities, the overall fit continúes to show 
the presence of a residual in the emission-line wings. Following 
the procedures outlined in Firpo et al. (2010), we are able to fit 
a broad component, with a velocity dispersión of about 78kms^' 
from the Ha emission line and slightly lower from the rest of the 
Unes. Table 2 shows the parameters for the three components that fit 
the global profile. Individual Gaussian component fluxes are Usted 
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Figure 2. Strongest emission Unes frora the Haro 15 A spectrum. Each panel includes a flux-calibrated spectrum, where the individual A:-axis has been 
normalized to the observed radial velocity for comparison purposes. To enhance details at low-luminosity levéis, the y-axes are shown in logarithmic scale. 
From top to bottom and left to right: H^, [O iii] 15007, Ha, [N ii] A6584 and [S ii] 16717. 



as fractional emission measures (EMf) relative to the total line flux 
following the work by Relaño & Beckman (2005). The sum of 
these individual fluxes, which we will hereafter refer to as overall 
Ha flux, uncorrected for reddening, is found to be (5.49 ± 0.05) x 
10"'"* ergs"' cm"^. 

In Fig. 2 we show the ngaussfit fitting done with three different 
Gaussian components in the emission lines which have enough sig- 
nal to provide a reliable fit. The vahdity of the profile multipücity 
and broadening is checked over the different emission Unes, be- 
coming more evident for the strongest emission lines. Fig. 3 shows 
the excellent agreement among individual fits for the most intense 
emission lines in knot A. 

The top panel in Fig. 4 shows the spatial profile of the Ha emis- 
sion line in knot A. Two distinct zones are clearly distinguished, 
separated by 1.72 arcsec ~ 0.72 kpc. We have performed additional 
extractions of these two zones (labelled A I and A II) individually. 
Inspection of the individual plots displayed in the bottom panels 
show that the kinematical structure is dominated by región AI, 



somewhat expected as it is brighter than región A II, which exhibits 
a more simple emission profile in turn. High spatial resolution in- 
tegrated field spectroscopy is needed to disentangle this complex 
behaviour. 



3.1.2 Haro 15 B 

The spectrum of región B is dominated by nebular emission lines, 
where we detected and measured [O iii] 14363 and He ii 14686. We 
identified and fitted Gaussian profiles to the H/!, [O ni] 114959, 
5007, [Nii] 16548, Ha, [Nii] 16584 and [Sii] 116717, 6731 lines 
in this región using ngaussfit. It is interesting to note that, al- 
though from inspection of the Wide Field Planetary Camera 2 im- 
age shown in Fig. 1 a complex structure of región B is apparent, 
this is not reflected in the kinematical behaviour of the región, since 
it is possible to fit a single Gaussian component with a velocity 
dispersión of 20kms"'. Nevertheless, the profile fittings show the 
presence of a residual in the emission-line wings that we were able to 
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Velocity (km s ) 

Figure 3. Same set of Unes as plotted in Fig. 2 now overlapped in the same 
diagram. To enhance details at low-luminosity levéis, the v-axes are shown in 
logarithmic scale. The evident similarity among different lines belonging to 
different elements confirms the presence of a complex kinematical structure. 
The results of fitting of individual components for the [S ii] Xdlll forbidden 
Une are included as a reference. 



explain by fitting a broad component with a velocity dispersión of 
about 43kms^'. Results of the fitting procedure can be seen in 
Fig. 5; it is worth noting that the logarithmic scale chosen for the 
y-axis magnifies the 1 per cent residual in the profile wings for 
Ha and [O iii] X5007. For the listed lines the derived radial veloci- 
ties, velocity dispersions and their corresponding errors are listed in 



Table 2. The overall Ha flux, uncorrected for reddening, is found 
to be (9.88 ± 0.3) x lO"''' erg s"' cm-^. 

As concluded by López-Sánchez & Esteban (2009), knot B has 
a much higher ionization degree than the rest of the regions in 
Haro 15. These authors suggest that this may be a consequence of 
the extreme youth of this knot. This fact and the low metallicity 
found in knot B (lower than for the other regions) indícate that it 
has a different nature, probably being the remnant of a dwarf galaxy 
which is experiencing a minor interaction with Haro 15 (López- 
Sánchez 2010). 



3.1.3 Haro 15 C 

Although the observed spectrum of región C is faint, we were able 
to identify and fit Gaussian profiles to the H/í, [Oiii] A.5007, Ha, 
[Nii] >i6584and [Sii] AA.6717,6731 lines. Knot C shows one naiTow 
component together with the underiying broad component. How- 
ever, due to the poor signa! of the spectrum, all emission lines were 
fitted using the Ha ngaussfit solution as a témplate, fixing the pro- 
file centres and widths, and only allowing the task to fit the profile 
amplitudes. In Fig. 6 we only show Ha since the rest of the intense 
lines are very noisy due to the low signal-to-noise (S/N) ratio of 
the spectrum. Table 3 shows the derived kinematical parameters. 
The overall Ha flux, uncorrected for reddening, is (2.67 ± 0.07) x 
' ergs^' cm" 



3.1.4 Haro 15 E 

We detect the presence of a broad profile, suggesting the presence of 
a double peak together with a flux excess in the strongest emission 



Table 2. Results of Gaussian profile fitting to the observed emission lines in Haro 15 A and B. Each emission Une is identified by its ion laboratory wavelength 
and ion ñame in columns 1 and 2. According to the different fits performed on each Une, column 3 identifies each nan'ow component ( 1 and 2, where applicable), 
and a broad component. Radial velocities (V,-) and intrinsic velocity dispersions (aini) together with their respective errors are expressed in km s^ . The intrinsic 
velocity dispersions are corrected for the instrumental and thermal widths. Fractional emission measures (EMf) in per cent. 
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Figure 4. Top panel: intensity distribution along the slit at the peak of the Ha emission. Bottom panels: result of profile extraction at the locations of regions 
AI and A II, together with the result of individual Gaussian profile fitting. To enhance details at low-luminosity levéis, the y-axes are shown in logarithmic 
scale. Two narrow components and a broad one are fitted for A I and only a single narrow and broad component are identiñed in A II. 



Unes. We identified and fitted Gaussian profiles to the H/i, [Oiii] 
A5007, Ha, [Nii] A.6584 and [Sii] AX6717, 6731 Unes. The results 
of the profile fitting procedure yield two components of similar 
width. Although both are supersonic profiles, it is not possible to 
distinguish a broad and narrow component and therefore refer to 
the identified components as narrow 1 and 2 (see Fig. 6). For the 
Usted Unes the derived radial velocities, the velocity dispersions 
and the corresponding errors are Usted in Table 3. The overall Ha; 
flux, uncorrected for reddening, is found to be (4.99 ± 0.14) x 
10"" ergs"' cm"^. 

Owing to the peor S/N ratio spectrum we fitted some of the emis- 
sion Unes, such as [Nn] A6584 and [Sii] A.X6717, 6731, using the 
Ha NGAUSSFIT solution, and iterating separately for individual sets of 
parameters, foUowing the procedure described in Firpo et al. (20 10) . 



3.1.5 Haro 15 F 

The observed spectrum of knot F is the weakest one, and henee 
we only identified and fitted the Gaussian profiles of the H/S, [O iii] 
A,5007 and Ha Unes. H/í and [Oiii] A5007 are in fact very noisy 
and narrow, rendering impossible the task of fitting múltiple Gaus- 
sian components (narrow plus broad). The spectral Une broadening 
in their integrated spectra is only observed in the Ha Une where 
we could fit a narrow Gaussian component together with the broad 
one, with velocity dispersions of about 8.3 and 22kms"', respec- 
tively (see Fig. 6). The estimated velocity dispersión for the narrow 
component in the Ha Une is subsonic (8kms"'), typical of classic 
H II regions, although the presence of a broad component is almost 
exclusive of giant H ii regions. 



For Ha, the derived radial velocity, the velocity dispersions and 
the corresponding errors are Usted in Table 3. The overall Ha 
flux, uncorrected for reddening, is found to be (1.06 ± 0.16) x 
10"'^ ergs"' cm"^. 



3.2 Radial velocities 

The radial velocities for each component of the different knots of 
Haro 15 are given in the corresponding tables where the results 
for the profile fits are shown. In Table 4 we list the average ra- 
dial velocity for each component, together with the corresponding 
errors. The single component is the result of a single Gaussian 
fit to the emission-line profiles. In order to compare our results 
with those of López-Sánchez & Esteban (2009) we have plotted 
our radial velocities in Fig. 7 relative to the velocity of the galac- 
tic centre 6415 km s"' (López-Sánchez 2006), after correcting to 
the Galactic standard of rest (GSR) using the rvcorrect task by 
IRAF. Fig. 7 shows our results together with data presented by 
López-Sánchez & Esteban (2009) for their long slits with posi- 
tion angles (PA) of 41° and 117°. Regarding the single component 
velocity, it can be clearly seen that the average radial velocities 
derived from the emission Unes foUow nicely those derived by 
López-Sánchez & Esteban (2009). However, it is worth mentioning 
that individual velocities of narrow components (ni and n2) for 
knot A differ by about 60kms"'. This cannot be easily explained 
by scatter within galactic rotation, suggesting that their relative 
velocities are due to mutual orbital motion around their gravity 
centre. 
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Figure 5. ngaussfit fits with two Gaussian components in the Haro 15 B strongest emission-line profiles, a narrow componen! and a broad one. In order: H/i, 
[Oiii] A,5007, Ha, [Nii] 16584 and [Su] A6717. Note that the logarithmic >-scale raagnifies the actual low (1 per cent) error in the overall fits. 



3.3 Relation between Ha luminosities and velocity dispersión 

We also explored the influence that the detection of múltiple com- 
ponents has on the location of individual giant H ii regions in the 
Ha huminosity versus velocity dispersión (crim) diagram. Table 5 
lists the velocity dispersions (for a single Gaussian profile) and Ha 
luminosities derived from our data. These are plotted, together with 
data from Tables 2 and 3, in Fig. 8. Each knot is labelled from A to 
F as used throughout this paper. Narrow components are identified 
with 'n', while 's' refers to hnewidths measured by fitting a single 
Gaussian component to the Une. Also plotted are data on NGC 6070 
and 7479 from Firpo et al. (2010) and giant H ii regions from Bosch, 
Terlevich & Terlevich (2002) together with the regression line fitted 
by those authors to their 'young' giant H ii regions. Ha luminosities 
were derived from the fluxes measured directly from the component 
fitting to our echelle spectra, and using distances as published by de 
Vaucouleurs et al. (1991) for Haro 15. The correction for reddening 
could change flux valúes by a factor up to 1.72 (0.24 in logarithm; 
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López-Sánchez & Esteban 2009) but, as we cannot determine it for 
every región, we prefer to list the uncorrected ones. 

Similar to what was found by Firpo et al. (2010), the location of 
the studied regions is far from random. Haro 15 Hii regions show 
a correlation between their luminosities and velocity dispersions. If 
we focus on the fit of single components they seem to foUow the 
regression found for virialized systems. The outstanding exception 
to this trend is región A, but this can be attributed to the fact that 
the profile can be split into two narrow components, resembling 
the behaviour found for regions NGC 7479-1 and 7479-11 (Firpo 
et al. 2010). The narrow components of the composite fits show, as 
expected, relatively smaller luminosities and velocity widths, but 
they still lie around the same linear regression. In this case, it is 
región F that fails to follow the trend but, as mentioned earlier, its 
low velocity dispersión makes it a very unlikely candidate for a 
giant H II región. It can be seen from the figure that while narrow 
components tend to cluster around the L(Ha) oc. a"^ relation, sin- 
gle component fits would provide a flatter power-law exponent. 
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Figure 6. ngaussfit fits with three Gaussian components in the Haro 15 C, E and F Ha emission-line profiles, a narrow component and a broad one. 



Table 3. Results of Gaussian profile fitting to the observed emission Unes 
in Haro 15 C, E and F. Each emission Une is identified by its ion laboratory 
wavelength and ion ñame in columns 1 and 2. According to the different ñts 
performed on each Une, column 3 identifies the narrow and broad compo- 
nents. Radial velocities {Vj) and intrinsic velocity dispersions (a ¡ni) together 
with their respective errors are expressed in km s^ ' . The intrinsic velocity 
dispersions are corrected for the instrumental and thermal widths. Fractional 
emission raeasures (EMf ) in per cent. 
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This is probably due to the broad component contributing a 
substantial part of the total Hcc flux. Although not included in the 
figure, the broad components are located, as expected, in a parallel 
sequence shifted towards lower luminosities and/or higher velocity 
dispersions. 



4 SUMMARY 

From high-resolution spectroscopy of Haro 15 we have performed 
a thorough analysis of the emission-line profiles of several knots 
considering múltiple component fits to their profiles. Our results 
can be summarized as foUows. 

(i) Giant H ii regions of Haro 1 5 show a complex structure within 
the profile of all their emission Unes, detected in both recombination 
and forbidden Unes. 

(ii) The brightest emission Unes can be split into at least two 
strong narrow components plus an underlying broad component. 

(iii) Although regions tend to foUow the galaxy kinematics, the 
narrow components of knot A have relative velocities that are too 
large to be explained by galactic rotation. This behaviour can be 
explained if we are observing the orbital motion around their gravity 
centre. 

(iv) Almost all knots foUow the relation found between luminosi- 
ties and velocity dispersions for virialized systems, when consider- 
ing either single profile fitting or the strong narrow components in 
more complex fits. Among these, the single one shows a relatively 
flatter slope. 

The presence of more than one component in the Gaussian fits to 
the emission-line profiles, such as those analysed in this paper, has 
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Haro 15: internal kinematic 

Table 4. Average LSR radial velocities for Haro 15 emission-line knots. Column 1 indicates the fitted feature and the remaining 
columns list the average radial velocity for each knot and its uncertainty, both in km s^ ' . 
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Figure 7. Position versus velocity diagrams for the observed knots in Haro 15. Ha radial velocity results, relative to knot A, are superimposed on the data by 
López-Sánchez & Esteban (2009) for slit positions PA 41° (left-hand panel) and PA 117° (right-hand panel). The Ha and [O m] A5007 data of López-Sánchez 
& Esteban (2009) in Haro 15 are identiñed by blue circles and red triangles, respectively. 



Table 5. Intrinsic velocity dispersions (aim) and the Ha luminos- 
ity (L(Ha)s) for a ñtted single Gaussian component, together with 
their associated errors, expressed in kms^' and ergs^', respec- 
tively. We add the Ha luminosity (L(Ha)overall) obtained from the 
overall Ha flux, uncorrected for reddening, and its error, in erg s^ ' . 



Knot 



L(Ha)s 



L(Ha)overall 



Haro 15 A 
Haro 15 B 
Haro 15 C 
Haro 15 E 
Haro 15 F 



47.3 ± 0.9 

34.4 ± 0.8 

19.3 ±0.8 

24.4 ± 0.9 
11.6 ±0.8 



4.9 ± 0.7 X 10* 
8.6 ± 1.2 X lO'*" 
2.3 ± 0.3 X 10^*^ 
4.5 ± 0.6 X 10^*^ 
9.1 ± 1.3 X 1038 



4.9 ± 0.7 X 10''° 

8.9 ± 1.3 X 10''° 

2.4 ± 0.3 X lO^*^ 

4.5 ± 0.6 X lO^*^ 
9.5 ± 1.3 X 1038 



been discussed in several previous studies. Hágele et al. (2007, 2009, 
2010) showed fits that involved the existence of broad and narrow 
components for the emission Unes of the ionized gas in circumnu- 
clear star-forming regions. Furthermore, Firpo et al. (2010) were 
able to detect two distinct components within the narrow feature 
of the emission Unes in the brightest giant H ii regions NGC 7479-1 
and 7479-11. Although we still lack the number of objects needed 
to support a generalization of this trend, these data show that the 
presence of múltiple kinematical components among extragalac- 
tic star-forming regions cannot be overlooked. They might have a 
strong impact on subsequent analyses that rely on basic parameters, 



such as the velocity dispersión and chemical abundance of the ion- 
ized gas, the inferences about the nature and strength of the source of 
ionization or the classification of the activity in the central regions of 
galaxies. 

It is necessary to map these regions with high spectral and spatial 
resolution and much better S/N ratio to disentangle the origin of 
these different components. 2D spectroscopy performed with the 
Integral Field Unit is the ideal tool to tackle this issue. 
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Capítulo 4 

Determinación de abundancias del 

gas ionizado para las diferentes 

componentes cinemáticas 



El análisis de la distribución de densidad y temperaturas electrónicas en el interior de 
una nebulosa de gas ionizado permite analizar las propiedades físicas del medio interestelar 
y determinar las abundancias químicas iónicas y totales de dicho gas. La determinación de 
la abundancia de cada ion está asociada a su temperatura electrónica. En aquellos casos 
en que estas temperaturas no son medibles directamente a través de las líneas de emisión 
aurórales sensibles a la temperatura se suelen utilizar relaciones con otras temperaturas iónicas 
determinadas empíricamente o derivadas utilizando modelos de fotoionización. El análisis de 
los parámetros físicos tales como densidad, temperatura, composición química, etc., puede 
llevarse adelante mediante el estudio de las intensidades relativas de las líneas de emisión del 
gas. En este Capítulo se procede a determinar las propiedades físicas del gas en los distintos 
brotes de formación estelar de Haro 15 en cada componente cinemática ajustada en los perfiles 
de las líneas de emisión. Si bien este tipo de análisis desglosado para las diferentes componentes 
cinemáticas no es la primera vez que se realiza, este trabajo presenta una mayor completitud 
frente a otros trabajos publicados dada la calidad de nuestros datos y la gran cantidad de líneas 
analizadas, entre las que se encuentran las débiles líneas aurórales, sumado a la posibilidad 
de descomponer sus perfiles en más de una componente angosta y una componente ancha. Es 
allí donde radica el carácter relevante de este Capítulo para esta Tesis. 

4.1. Introducción 

Este Capítulo se centra en el estudio de regiones Hll de baja metalicidad como lo son las 
regiones ubicadas en las galaxias Enanas Compactas Azules (BCDs, sigla en inglés de Blue 
Compact Dwarf galaxies), que pertenecen a la clase de galaxias más numerosa del Univer- 
so: las galaxias enanas. Los espectros de las BCDs están dominados por la emisión del gas 
asociado a cúmulos estelares formados por estrellas jóvenes y calientes que ionizan el medio 
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interestelar, y estas galaxias se caracterizan por su color azul, aspecto compacto, alto conte- 
nido de gas, fuertes líneas de emisión nebular y abundancias químicas bajas (Kunth &¿ Ostlin, 
2000; Stasiñska et al., 2001). Gran parte de los objetos con espectros de emisión en el régimen 
de alta excitación, siendo este régimen de gran importancia para el estudio de la metalicidad 
y de las propiedades de la población estelar, son galaxias Hii. Las galaxias Hil son galaxias 
irregulares de baja masa con al menos un episodio de formación estelar violento concentrado 
en unos pocos parsecs cerca a su núcleo (Melnick et al., 1985aa, 1985bb). Presentan un alto 
brillo superficial, intenso color azul y baja metalicidad con espectros de emisión debidos al 
gas ionizado. Actualmente, se ha multiplicado el interés en estudiar las propiedades de es- 
te tipo de objetos menos metálicos, dado que existen indicios de que pudieran tratarse de 
objetos jóvenes y, por lo tanto, proporcionarían importantes pistas acerca de la formación 
estelar en estas condiciones, al margen de la trascendencia que tienen para la búsqueda de 
helio primordial. 

El término galaxia H ii se utiliza cuando se hace referencia a objetos seleccionados a partir 
de placas obtenidas mediante la técnica del prisma-objetivo cuyo criterio de búsqueda son las 
líneas de emisión (por ejemplo, Terlevich et al., 1991), mientras que las BCDs son objetos 
seleccionados en base a su color muy azul y su aspecto compacto en las placas fotográficas. De 
todos modos ambos términos suelen designar a los mismos objetos, ya que sus propiedades 
observables en el óptico están dominadas por el espectro de la región Hii en la zona de 
formación estelar. 

Los primeros trabajos sobre estos objetos compactos son los de Haro (1956) y Zwicky et 
al. (1965) que distinguen estas galaxias de gran brillo superficial en sus placas fotográficas. 
Sargent & Searle (1970) encuentran que casi todos estos objetos presentan un espectro domi- 
nante de líneas de emisión como el de las regiones Hii en los discos de las galaxias espirales 
y las denominan Regiones Hii Gigantes Extragalácticas aisladas. Generalmente se componen 
de una región central de unos pocos cientos de parsecs de diámetro con un elevado brillo 
superficial en la cual se están formando miles de nuevas estrellas en un proceso de formación 
estelar violento (starburst, de su significado en inglés) (Searle & Sargent, 1972) y de una 
galaxia de cientos de parsecs de diámetro de baja luminosidad (M„ r¿-17), pudiendo alcanzar 
unos pocos kpc cuando se la observa en los filtros V y R, excediendo el tamaño de la región 
de formación estelar (Cairos et al., 2001). 

En general, las líneas de emisión observadas en los espectros tanto de las galaxias Hii 
como de las BCDs son muy similares a las de Regiones Hii Gigantes Extragalácticas. Por 
lo tanto, mediante la aplicación de las mismas técnicas de medición que en las regiones Hii, 
pueden ser derivadas las temperaturas, densidades y composición química del gas interestelar 
en los brotes de formación estelar de este tipo de galaxias deficientes en metales (Sargent & 
Searle, 1970; Frcnch, 1980; Terlevich et al., 1991). Sometidas a intensos brotes de formación 
estelar, espectros con intensas y angostas líneas de emisión, extinción de polvo relativamente 
pequeña, entornos de baja metalicidad y diferentes historias de formación estelar, hacen que 
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las galaxias BCDs sean los más apropiados laboratorios para estudiar la metalicidad de las 
galaxias (Kmith & Óstlin, 2000). 

En el Capítulo anterior se presentó un estudio detallado de la cinemática interna del 
material nebular en múltiples regiones de formación estelar de la galaxia BCD HarolS. De 
la espectroscopia de alta resolución, obtenida en el Observatorio Las Campanas, Chile, se ha 
realizado un análisis exhaustivo de sus líneas de emisión, incluyendo los ajustes de múltiples 
componentes a sus perfiles de emisión. Los resultados obtenidos han demostrado que las 
regiones Hii de Harol5 presentan una estructura compleja en los perfiles de sus líneas de 
emisión, detectada tanto en las líneas de recombinación del hidrógeno y helio como también 
en las líneas prohibidas. La emisión de la región más brillante puede ser dividida en, al 
menos, dos componentes angostas más una componente ancha. Aunque las regiones tienden 
a seguir la cinemática de la galaxia, las componentes del brote A tienen velocidades relativas 
demasiado grandes para ser explicadas por la rotación galáctica. La mayoría de las regiones 
siguen la relación encontrada entre la luminosidad y la dispersión de velocidades para sistemas 
virializados, ya sea cuando se considera un ajuste de perfil simple o un ajuste más complejo 
con varias componentes, excepto cuando se considera la componente ancha del ajuste. Entre 
todas estas componentes, la componente simple muestra una pendiente relativamente más 
plana que el resto. 

En general, la calibración en ñujo de datos échelle es difícil de comprobar y las líneas 
débiles, sensibles a la temperatura, tienen una relación señal-ruido relativamente baja para 
permitir una medición confiable. Los datos espectroscópicos de ranura larga con resolución 
intermedia son más sensibles para realizar un análisis de las abundancias. Dichas abundancias 
dependen de la precisión de las temperaturas estimadas. La combinación del échelle y de 
ranura larga brinda la oportunidad de interpretar el resultado de la abundancia a la luz de la 
cinemática interna de las regiones. 

Shi et al. (2005) encuentran un promedio de 8.56 en el logaritmo de las abundancias 
totales de oxígeno de Haro 15, para una temperatura electrónica de 8330 K. Sin embargo, se 
sabe por un trabajo anterior (Cairos et al., 2001) que la imagen Ha muestra una morfología 
nudosa con la región estelar resuelta en un gran número de grupos que aparecen dispersos en 
toda la galaxia. En López-Sánchez &: Esteban (2009b) los autores muestran espectroscopia 
de los diferentes nudos espacialmente resueltos y determinan densidades, temperaturas y 
abundancias químicas de los brotes A, B y C de Haro 15 (siguiendo la nomenclatura de 
Cairos et al., 2001). Las abundancia derivadas de oxígeno fueron de 12+log(0/H) = 8.37ib0.10 
y 8.10ib0.06 para las regiones A y B, respectivamente. 

Esteban &¿ Vilchez (1992) introducen el término "chemodynamics" al conbinar por primera 
vez, información química y cinemática en orden de realizar un estudio químico-dinámico de 
la nebulosa anillo W-R, NGC 6888. Los autores combinan observaciones de alta resolución 
espectral con observaciones de alta resolución espacial para las diferentes zonas de la nebulosa, 
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lo cual les permite analizar las diferentes componentes en velocidad de la región estudiada. 

Más tarde, James et al., (2009), realizaron un estudio en la galaxia BCD Mrk 996 basado 
en alta resolución espectral con VLT VIMOS en modo IFU. Encontraron que la galaxia Mrk 
996 presenta multicomponentes en los perfiles de sus líneas de emisión, descomponiendo los 
mismos en una componente angosta y una componente ancha, llevando a cabo un análisis 
de las propiedades físicas separado por componentes. Según los autores esta peculiar galaxia 
BCD presenta en el núcleo zonas extremadamente densas de gas, a través del cual los outfiows 
estelares y posibles frentes de choque contribuyen a la excitación de la emisión ancha. Ellos 
detectaron que las líneas de Balmer muestran una estructura común en velocidad de multi- 
componentes anchas que no es observable en las líneas prohibidas. Las líneas prohibidas con 
excepción de [Oiii] A4363 y [Nii] A5755, pudieron ser ajustadas con una componente angosta 
y una componente ancha. En cambio el perfil de emisión de ninguna de las líneas aurórales 
sensibles a la temperatura, como por ejemplo las líneas de [Oiii] y [Nii] que acabínamos de 
mencionar, pudo ser descompuesto en más de una componente, sólo en estas dos líneas au- 
rórales los perfiles de emisión pudieron ser ajustados con una componente ancha. Si bien el 
análisis de abundancias por componentes presentes en los perfiles de las líneas de emisión 
es llevado a cabo por James y colaboradores, la diferencia fundamental con nuestro traba- 
jo radica en que nosotros pudimos descomponer algunas de las líneas aurórales sensibles a 
la temperatura electrónica en sus diferentes componentes cinemáticas suponiendo que dicha 
solución de multicomponentes cinemáticas es la misma que se encuentra en las líneas más 
intensas del espectro con igual grado de ionización. 

En este Capítulo se realiza un análisis de las características del gas ionizado discriminado 
para cada componente y para el fiujo global de las líneas. Para esto se estimaron temperatu- 
ras y densidades electrónicas, y se obtuvieron abundancias iónica y totales de las diferentes 
especies como ser: O, S, N, Ne, Ar. 

En lo que sigue, se explicarán las observaciones y la reducción de los datos, y a continuación 
se mostrarán los resultados obtenidos. En la Sección 4.4 se explicarán los diferentes métodos 
utilizados para determinar las condiciones físicas del gas ionizado (densidades y temperaturas 
electrónicas), como las abundancias químicas derivadas y sus incertezas. En la Sección 4.5 
se discuten los datos dando las respectivas conclusiones. Y finalmente, en la Sección 4.6 se 
presenta un resumen y las conclusiones del Capítulo. 

4.2. Observaciones y Reducción de los datos 

Se obtuvieron espectros de resolución intermedia de ranura larga usando la cámara Wide- 
Field CCD (WFCCD) montada en el telescopio du Pont de 100" del Observatorio Las Campa- 
nas (LCO) en Chile (28 de septiembre de 2005). Las observaciones obtenidas corresponden a 
las regiones B y C de Haro 15 (ver Figura 1 del Capítulo 3). El detector TEK5 que se utilizó cu- 
bre todo el rango de longitud de onda óptica, desde 3800 a 9300 A (centrado en Xc = 6550 Á). 
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La resolución espectral alcanzada en los espectros es de Rjt'vkíím — 900: AAfvfíím = 7.5 A a 
los A 6700 A, medida a partir del ancho a potencia mitad (FWHM) de las líneas de la lámpara 
de comparación de Thorio- Argón tomadas con el fin de calibración de longitud de onda. Las 
observaciones se realizaron en ángulo paraláctico, para evitar los efectos de refracción dife- 
rencial en el ultravioleta, en excelentes condiciones fotométricas, con un seeing de 1 segundo 
de arco. Al comienzo de la noche se tomaron bias (nivel de pedestal o cuentas que tiene de 
base el detector aún sin exponer) y ñat-field (del inglés, imagen de campo plano) de cúpula. 
Además, en esta misma noche se tomaron espectros de la lámpara de comparación Th-Ar. Las 
imágenes fueron procesadas y analizadas con las rutinas de IRAF^ mediante el procedimiento 
usual. El procedimiento incluye la eliminación de rayos cósmicos, resta de bias, división por 
un ñat-field normalizado y calibración en longitud de onda. Se observó una estrella estándar, 
EG131, para la calibración en flujo y el tiempo de exposición para esta estrella estándar de 
flujo fue de 180 segundos. Los espectros fueron corregidos por extinción atmosférica. 

Se ha encontrado un cierto grado de contaminación de segundo orden en estos espectros de 
ranura larga. La misma comienza a estar presente en A > 6000 A y es casi imposible de eliminar 
por completo debido a la amplia cobertura espectral usada. Por este motivo, los flujos medidos 
para las líneas nebulares que se encuentran más allá de los 6000 A han sido sistemáticamente 
subestimados, aunque nunca esta contribución es importante para los cocientes de líneas (ver 
Terlevich et al., 1991). 

Por otro lado, se obtuvieron espectros de alta resolución usando el espectrógrafo échelle 
en el telescopio du Pont de 100" de LCO, entre el 19 y el 20 de julio de 2006. Las observaciones 
obtenidas corresponden a varios brotes de formación estelar de la misma galaxia Haro 15 (A, B, 
C, E y F) adquiridos en el modo échelle (ver Capítulo 3). El desplazamiento al rojo calculado 
con estos datos es de z = 0.021306. El rango espectral abarcado en estas observaciones es 
desde 3400 hasta 10000 A. Este rango espectral garantiza la medición simultánea de las líneas 
nebulares de [Olí] AA 3727, 3729 Á a las del [Siii] AA9069, 9532 Á en ambos extremos del 
espectro y en la misma región de la galaxia. Las condiciones de observación fueron buenas 
con un seeing de 1 segundo de arco. La estrella espectrofotométrica estándar observada para 
la calibración en flujo fue FeigellO (Bohlin et al., 2001), la cual tiene sus flujos tabulados 
cada 2 A, y la cantidad de intervalos definidos dentro de un orden del échelle varía de cuatro a 
doce, dependiendo de la calidad del espectro. El tiempo de exposición de la estrella estándar 
fue de 1200 segundos. La calibración en flujo obtenida en estos datos fue muy buena. Para 
los detalles de reducción y el proceso de análisis ver Capítulo 3. 

La Tabla 4.1 lista las observaciones obtenidas, junto con los tiempos de exposición y las 
masas de aire para cada región y su respectivo modo de observación. 



^Image Reduction and Analysis Facility, distribuida por NOAO, operada por AURA, Inc., bajo acuerdo 

con NSF. 
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Tabla 4.1: Observaciones de los diferentes brotes en HarolS. La primera columna indica el 
modo de observación, la segunda columna muestra la nomenclatura utilizada en el Capítulo 
3 siguiendo Cairos et al., (2001), la tercera columna corresponde a la fecha de observación, 
la cuarta columna exhibe el tiempo de exposición de los brotes, la quinta columna señala las 
masas de aire. 



Modo 


brote 


fecha 


t.exp 


sec z 


ranura larga 


B,C 


2005 Sep 28 


2x1200 + 1x900 


1.5 


échelle 


A 


2006 Jul 19 


1800 


1.2 




B 


2006 Jul 19 


1800 


1.1 




C 


2006 Jul 19 


1800 


1.1 




E 


2006 Jul 20 


1800 


1.2 




F 


2006 Jul 20 


1800 


1.1 



4.3. Resultados 

En la Figura 4.1 se han marcado algunas de las líneas de emisión más relevantes observadas 
en los espectros WFCCD de ranura larga de los brotes B y C de Haro 15. 

Los flujos de las líneas de emisión para los datos de ranura larga se midieron utilizando 
la tarea splot de IRAF (que integra la intensidad de cada línea sobre un continuo ajustado 
localmente) siguiendo el mismo procedimiento descripto en Hágele et al. (2006). Se utilizaron 
dos formas de integrar el flujo de la línea de emisión: (i) en el caso de una línea aislada o dos 
líneas fusionadas y sin resolver, la intensidad es calculada integrando entre dos puntos dados 
por la posición del continuo local colocados a ojo. A este flujo, que incluye todas las diferentes 
componentes cinemáticas, se lo denominará, de ahora en más, como "medida global"; (ii) si 
dos líneas están fusionadas, pero se pueden resolver, se utiliza el procedimiento de múltiples 
ajustes Gaussianos para estimar los flujos individuales (esto se aplicó solamente para los datos 
de dispersión simple). Ambos procedimientos están explicados detalladamente en Hágele et 
al., (2006; 2008). 

Siguiendo Pérez-Montero &: Díaz (2003), los errores estadísticos asociados a los flujos de 
las emisiones observadas se han calculado utilizando la expresión 

ai = acN^^^[l + EW/{NA)]^/^ (4.1) 

donde ai es el error en el flujo de la línea observada, ctc representa la desviación estándard en 
una caja cerca de la línea de emisión medida y representa el error en la posición del continuo, 
N es el número de píxeles en la medida del flujo de la línea, EW es el ancho equivalente de la 
línea el cual está dado por el flujo del continuo cerca de la línea y el flujo de la línea, y A es 
la dispersión de la longitud de onda en angstroms por pixel (González-Delgado et al., 1994). 
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Figura 4.1: Espectros de líneas de emisión de ranura larga obtenidos para los brotes B y C 
(paneles superior e inferior, respectivamente) sin corrección Doppler. 
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Cada intensidad medida debe ser corregida de enrojecimiento. La constante de enrojeci- 
miento c{H/3) se puede calcular a partir del decremento de Balmer de las líneas de recombi- 
nación de hidrógeno más intensas. Una población estelar subyacente es fácilmente apreciable 
debido a la presencia de alas en las absorciones que deprimen las líneas de emisión de Balmer 
en los espectros de ranura larga y échelle. Para medir las intensidades de las líneas y minimizar 
los errores introducidos por la población subyacente, se define un pseudo-continuo en la base 
de las líneas de emisión del hidrógeno (ver discusión en Hágele et al., 2006). La presencia de 
alas en las líneas de absorción implica que, a pesar de que se ha utilizado un pseudo-continuo, 
aún existe una fracción absorbida del flujo emitido que no es posible medir con precisión (ver 
la discusión al respecto en Díaz, 1998). Esta fracción no es la misma para todas las líneas, 
ni lo son los cocientes entre las fracciones absorbidas y la emisión. Hágele et al. (2006) esti- 
maron que la diferencia entre los valores obtenidos usando el pseudo-continuo deflnido o un 
ajuste múltiple Gaussiano de las componentes de absorción y emisión es la misma, dentro 
de los errores, para todas las líneas de Balmer. En cualquier caso, siguiendo este trabajo, 
para las líneas de emisión de Balmer se ha duplicado el error derivado, ai, como un enfoque 
conservador para incluir la incertidumbre introducida por la presencia de la población estelar 
subyacente. 

La ley que permite obtener la constante de enrojecimiento a partir del decremento de 
Balmer de las líneas de recombinación más intensas del hidrógeno y que una vez obtenidas, 
permite corregir el resto de las líneas de emisión, tiene la siguiente forma: 

-^o(^) = ^(^) . in-c(/í/3)-[/{A)-/{íí/3)] /. 2^ 

Io{H/3) I{HI3) ^ • ^ 

donde /(A) es la intensidad de la línea que se quiere corregir e /o(-^) es el valor corregido en 
relación a H(3, c{H¡3) es la constante de enrojecimiento y /(A) es la función de extinción. Los 
valores teóricos de las líneas de recombinación del hidrógeno que permiten calcular c{Hf3) son 
función de la temperatura electrónica y de la densidad, lo cual implicará que se debe realizar 
un procedimiento iterativo para calcularlo. 

El coeficiente de enrojecimiento ha sido calculado asumiendo la ley de extinción galáctica 
de Miller &: Mathews (1972) con i?v=3.2, que es el valor para la Galaxia, y la cual toma la 
forma: 

f ÍA) - f(HB) = í °'^^^^"' " ^'^°^' ^"' - ^'^^ 
^ ^ I 0,342A-i - 0,645, A^^ > 2,29 

donde A es la longitud de onda de la línea a calcular en unidades de micrones. 

Luego, la constante de enrojecimiento, C(H/3), fue obtenida mediante la realización de 
un ajuste de mínimos cuadrados a la relación entre F(A) y F(H/3) para los valores teóricos 
calculados por Storey & Hummer (1995), usando el método iterativo que aplicó Hágele (2008) 
para estimar Te y n^ según sea el caso. Los valores iniciales tomados, cuando éstos están 
disponibles, son aquellos derivados de la medida de los flujos de las líneas [Sil] AA 6717,6731 Á 
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y [Olll] AA4363, 4959, 5007 A. Debido a los grandes errores introducidos por la presencia de 
la población estelar subyacente, solo se han tenido en cuenta las cuatro líneas de emisión de 
Balnier más intensas (Ha, H/3, H7 y H¿). 

Por otro lado, los flujos asociados a cada componente cinemática identificadas en los datos 
de alta resolución espectral (échelle) para las líneas de emisión de mayor cociente señal/ruido 
(S/N, sigla en inglés de signal-to-noise ratio) han sido medidos como se detalla en el Capítulo 
3. Para aquellas líneas donde el cociente S/N no es suficiente, el proceso iterativo utilizado 
para deconvolucionar las diferentes componentes de los perfiles de las líneas no ha dado un 
resultado significativo. Sobre éstas líneas más débiles y utilizando la tarea ngaussfit de 
IRAF, se ha copiado la solución encontrada del ajuste realizado en las líneas más intensas 
con un grado de ionización comparable. Luego, fijando los centros (a la longitud de onda 
correspondiente de la línea que se quiere ajustar) y anchos de las componentes, sólo se le 
permite a la tarea variar simultáneamente la amplitud de cada componente. 

Se han supuesto dos zonas de ionización: una zona de baja ionización, cuya temperatura 
característica sería la de [Olí] y de donde se emiten las líneas de recombinación de hidrógeno 
y las líneas prohibidas de [Olí], [Nii], [Sil], y otra de alta ionización caracterizada por la 
temperatura de [Oiii] de donde procederían las emisiones de las líneas prohibidas de [Oili], 
[Neiii] y las líneas de recombinación del Hei. Se supone que tanto las líneas del [Slll] como las 
del [Ariii] son emitidas en una zona intermedia a las dos primeras por lo que se han realizado 
los ajustes a estas líneas utilizando ambas aproximaciones iniciales y así poder ver cuál es 
la que produce errores más pequeños. En los paneles superiores de las Figuras 4.2 y 4.3 se 
muestran el plano del Flujo- Velocidad de la línea de emisión más intensa [Oiii] A 5007 A en los 
brotes A y B, respectivamente, con las componentes ajustadas al perfil de línea superpuestas. 
Para mostrar la buena correlación de los ajustes entre los perfiles de línea de emisión intensas 
contra los perfiles de las líneas con menor señal-ruido, se graficaron en el panel inferior de 
cada figura un ejemplo del ajuste obtenido para una línea de emisión débil con similar grado 
de ionización que [Oiii] A 5007 A. En el brote A, se compara [Oiii] A 5007 A con Hei 5876 A, 
mostrando el buen ajuste realizado en esta línea de emisión débil. En el brote B, es notable 
la buena correlación entre las dos líneas de emisión del oxígeno. Gracias a este proceso es 
posible ajustar y deconvolucionar la importante y débil línea auroral [Oiii] A 4363 A, la cual 
es una línea sensible a la temperatura electrónica del medio donde se forma. 

Los flujos de las líneas de emisión de los cinco brotes observados son mostrados en la Tabla 
4.2 para los datos de ranura larga, y en las Tablas 4.3, 4.4, 4.5, 4.6 y 4.7 para los datos échelle. 
Cada Tabla lista la intensidad de las líneas de emisión para cada componente corregidas por 
enrojecimiento, junto a la constante de enrojecimiento (C(H/3)) y su correspondiente error 
considerado como la incerteza del ajuste de mínimos cuadrados y la intensidad de H/3 corregida 
por enrojecimiento. La columna 1 muestra la longitud de onda y el nombre de la línea medida. 
En la columna 2, la curva de enrojecimiento adoptada, f(A), normalizada a H/3. La columna 3 
el ancho equivalente (EW) en A. En la columna 4 la intensidad de la línea (/(A)) corregida por 
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Haro 15 A 
[OIII]5007 



flujo "global" 




Velocidad (km s ) 




Velocidad (km s ) 

Figura 4.2: Comparación entre las soluciones de distintos ajustes encontrados para una línea 
intensa y otra débil del brote A 
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Haro 15 B 
[OIII]5007 




Velocidad (km s ) 




Velocidad (km s ) 

Figura 4.3: Comparación entre las soluciones de distintos ajustes encontrados para una línea 
intensa y otra débil del brote B 
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enrojecimiento y relativa a 10000 H/3. Los errores relativos en las líneas de emisión listados en 
la columna 5 se obtuvieron mediante la propagación cuadrática de los errores observacionales 
en los flujos de las líneas y la incertidumbre en la constante de enrojecimiento. No se han 
tenido en cuenta los errores de las intensidades teóricas ya que son mucho menores que los 
errores de observación. 

En el espectro de ranura larga para el caso del brote B, todos los parámetros físicos que 
dependen de la intensidad de la línea [Siii] A 9532 A fueron calculados utilizando la relación 
teórica entre esta línea y la línea [Siii] A9069Á, 1(9069) ss 2.44 x 1(9532) (Osterbrock, 1989), 
ya que esta línea está fuera del rango espectral de observación. En el brote C fue imposible 
cuantiflcar una medida de estas líneas. En este caso, no se incluyeron estas líneas ni en las 
tablas ni en nuestros cálculos. También para este mismo brote, es imposible medir la línea 
[Nii] A 6548 A, ya que se solapa con Ha. Así que se tuvo que corregir el flujo de Ha de la 
contaminación proveniente de esta línea de nitrógeno estimando su contribución a partir de 
[Nll] A6584Á utilizando su relación teórica, 1(6584) « 2.9 • 1(6548). 

Debido al desplazamiento al rojo que presenta esta galaxia, en todos los espectros échelle, 
la línea [Sil] A 6731 A está afectada por la emisión de líneas telúricas de cielo de OH. Entonces, 
se procedió a utilizar la estrella estándar FeigellO (sdOB) como estándar telúrica, dado que 
la misma no presenta líneas estelares en ese rango, que fue observada con la misma masa de 
aire que la galaxia y que pudo ser muy bien calibrada en flujo. Así es que, haciendo uso de 
la rutina telluric de IRAF, se pudieron corregir los espectros de los brotes observados en 
el modo échelle por el efecto que producen las líneas telúricas. Luego esta línea fue ajustada 
con la tarea ngaussf it copiando la solución del ajuste obtenido para la línea [Sil] A 6717 A, 
y usando el mismo procedimiento que se ha usado para las líneas de bajo S/N. Luego, este 
será un importante punto a tener en cuenta al momento de analizar las densidades electrónicas 
derivadas para las diferentes componentes en estos brotes. Aquellas líneas que presentaban 
errores relativos en flujo mayores al 45 % no fueron tenidas en cuenta para los cálculos. La 
única excepción ha sido la línea [Olll] A 4363 A en la componente angosta del brote B que 
tiene un error mayor al 45 % pero se la consideró para poder calcular de temperatura del 
[Oiii]. En el brote F, solo se ha podido obtener la medida global dado que las líneas no tienen 
el suficiente S/N para realizar un ajuste por componentes. En ningún caso fue posible medir 
las líneas del [Oiii] AA 7319, 7330 Á. Por último, en los brotes C, E y F, debido al bajo S/N de 
estos tres espectros, las líneas [Siii] AA9069, 9532 A no pudieron ser medidas y, por otro lado, 
el flujo de la línea [Oiii] A 4959 A fue calculado a partir de la línea [Oiii] A 5007 A de acuerdo 
a la relación teórica entre ambas líneas, 1(5007) ~3x 1(4959) (Osterbrock, 1989). 

Para el brote A, no se tuvo en cuenta la línea Há debido a que la misma es muy débil y no 
puede ser descompuesta en las diferentes componentes cinemáticas. En el brote E, la segunda 
componente angosta y la medida global presentan un valor de C(H/3) un poco menor a cero, 
pero compatible con éste dentro de los errores observacionales, por eso fue tomado igual a 
cero. Además en el brote F, la línea de recombinación H/3 tiene un bajo cociente S/N, por lo 
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Tabla 4.2: Líneas de emisión medidas en el espectro de ranura larga en los brotes B y 
C obtenidas para la medida global. Intensidad de la línea corregida por enrojecimiento 
[F(/í'/3)=/(ií/3)=10000] 









B 






C 




A(Á) 


f(A) 


-EW 


/(A) 


Error 


-EW 


/(A) 


Error 






(A) 




(%) 


(A) 




(%) 



3727 [Olí]" 


0.271 


109.1 


11943 ± 117 


1.0 


3835 H9 


0.246 


3.8 


465 ± 23 


5.0 


3868 [Neiii] 


0.238 


34.0 


4494 ± 11 


0.2 


3889 Hei+H8 


0.233 


9.4 


1263 ± 65 


5.1 


3970 [Neiii]+He 


0.215 


21.1 


2535 ± 4 


0.2 


4068 [Sil] 


0.195 


0.7 


87 ± 17 


19.6 


4102 RS 


0.188 


21.1 


2395 ± 31 


1.3 


4340 H7 


0.142 


41.5 


4483 ± 17 


0.4 


4363 [Oiii] 


0.138 


6.3 


694 ± 42 


6.0 


4471 He: 


0.106 


2.8 


311 ± 16 


5.2 


4658 [Feíii] 


0.053 


0.3 


36 ± 6 


17.8 


4686 Heii 


0.045 


1.5 


154 ± 19 


12.1 


4713 [Ariv]+Hei 0.038 


0.6 


64 ± 7 


11.4 


4740 [Ariv] 


0.031 


0.9 


87 ± 12 


14.1 


4861 H/3 


0.000 


114.1 


10000 ± 25 


0.2 


4959 [Oiii] 


-0.024 


221.3 


21220 ± 77 


0.4 


5007 [Oiii] 


-0.035 


539.1 


51342 ± 88 


0.2 


5876 He: 


-0.209 


22.5 


1235 ± 28 


2.3 


6300 [01] 


-0.276 


6.8 


366 ± 16 


4.3 


6312 [Sill] 


-0.278 


2.3 


118 ± 18 


15.5 


6364 [01] 


-0.285 


2.2 


117 ± 11 


9.6 


6548 [Nii] 


-0.311 


8.3 


416 ± 20 


4.8 


6563 Ha 


-0.313 


565.8 


27618 ± 36 


0.1 


6584 [Nii] 


-0.316 


23.4 


1175 ± 22 


1.8 


6678 He: 


-0.329 


5.5 


264 ± 13 


4.7 


6717 [Su] 


-0.334 


32.1 


1453 ± 24 


1.7 


6731 [Su] 


-0.336 


13.7 


617 ± 49 


8.0 


7065 He: 


-0.377 


6.3 


251 ± 18 


7.2 


7136 [Ariii] 


-0.385 


19.4 


707 ± 21 


3.0 


7281 Hei^ 


-0.402 


1.1 


42 ± 5 


11.9 


7319 [Oll]= 


-0.406 


4.0 


148 ib 6 


3.9 


7330 [011]"^ 


-0.407 


3.1 


116 ± 4 


3.2 


7751 [Ariii] 


-0.451 


5.8 


186 ± 13 


7.0 


9069 [Slll] 


-0.561 


69.9 


1202 ± 118 


9.8 


I(H/3)* 






4.49 X 10-14 




c(H/3) 






0.29 ± 0.01 





45.3 25406 ± 2104 



6.0 3853 ± 657 



3.8 2842 ± 566 



1.6 1261 ± 330 



14.6 10000 ± 541 
8.3 6682 ± 209 
27.9 22042 ± 467 



2.5 



1542 ± 295 



52.6 28500 ± 824 

5.8 3230 ± 372 

10.2 5209 ± 613 

5.1 2609 ± 493 



0.11 X 10-14 
0.24 ± 0.12 



" [Olí] AA 3726 + 3729; " 
'^ [Olí] AA 7318 + 7320; '^ 
* ergseg-i cm"^ 



17.1 
19.9 
26.1 



19.1 



2.9 

11.5 

11.8 
18.9 



posible doblete con una línea desconocida 
[Olí] AA 7330 + 7331. 
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Tabla 4.4: Líneas de emisión medidas en el espectro échelle en el brote B obtenidas en la 
medida global, y en la componente angosta y en la ancha. Intensidad de la línea corregida 
por enrojecimiento [F{Hl3)=I{H ¡3) =10000] 









global 






angosta 






ancha 




A(Á) 


f(A) 


-EW 


/(A) 


Error 


-EW 


7(A) 


Error 


-EW 


7(A) 


Error 






(Á) 




(%) 


(A) 




(%) 


(Á) 




(%) 


3727 [Olí]" 


0.271 


20.4 


9049 ± 236 


2.6 


8.1 


24096 ± 854 


3.5 


8.0 


5877 ± 310 


5.3 


3835 H9 


0.246 


1.0 


358 ± 16 


4.4 














3868 [Neiii] 


0.238 


12.4 


5585 ± 184 


3.3 


1.6 


4214 ± 649 


15.4 


8.4 


5845 ± 271 


4.6 


3889 Hei+H8 


0.233 


5.1 


1571 ± 40 


2.6 














3968 [Nein]+H7 


0.216 


5.2 


1498 ± 30 


2.0 


0.5 


1016 ± 375 


36.9 


2.8 


1501 ± 163 


10.8 


3970 [Nein]+He 


0.215 


3.8 


1062 ± 57 


5.3 


1.2 


1626 ± 242 


14.9 


3.4 


1260 ± 104 


8.2 


4102 Há 


0.188 


9.7 


2766 ± 118 


4.3 


1.4 


2219 ± 338 


15.2 


5.7 


2534 ± 113 


4.5 


4340 H7 


0.142 


37.8 


4827 ± 122 


2.5 


5.6 


5305 ± 159 


3.0 


16.7 


4434 ± 56 


1.3 


4363 [Oiii] 


0.138 


4.4 


783 ± 127 


16.3 


0.4 


465 ± 213 


45.8 


2.4 


720 ± 90 


12.5 


4471 Hei 


0.106 


3.0 


471 ± 17 


3.6 














4713 [Ariv]+Hei 


0.038 


2.6 


336 ± 28 


8.4 














4861 H/3 


0.000 


81.3 


10000 ± 98 


1.0 


14.1 


10000 ± 157 


1.6 


45.6 


10000 ± 65 


0.7 


4959 [Oiii] 


-0.024 


130.0 


22714 ± 175 


0.8 


16.7 


14958 ± 286 


1.9 


91.0 


25569 ± 150 


0.6 


5007 [Oiii] 


-0.035 


470.0 


69812 ± 213 


0.3 


35.0 


37507 ± 795 


2.1 


233.8 


79550 ± 336 


0.4 


5876 Hei 


-0.209 


12.4 


1053 ± 9 


0.9 


1.7 


571 ± 37 


6.5 


10.2 


1237 ± 20 


1.6 


6300 [01] 


-0.276 


2.7 


181 ± 8 


4.4 


2.0 


501 ± 38 


7.5 


0.8 


73 ± 22 


30.3 


6312 [Siii] 


-0.278 


1.8 


125 ± 25 


20.2 


0.7 


195 ± 43 


22.2 


1.0 


103 ± 25 


24.1 


6548 [Nii] 


-0.311 


2.4 


185 ± 5 


2.6 


1.6 


422 ± 24 


5.7 


1.6 


158 ± 21 


13.0 


6563 Ha 


-0.313 


358.0 


28618 ± 40 


0.1 


73.5 


29349 ± 267 


0.9 


184.8 


28008 ± 130 


0.5 


6584 [Nn] 


-0.316 


7.7 


623 ± 12 


1.9 


6.0 


1573 ± 53 


3.4 


4.0 


404 ± 25 


6.2 


6678 Hei 


-0.329 


4.8 


320 ± 3 


1.1 


0.5 


129 ± 27 


20.7 


4.6 


427 ± 16 


3.8 


6717 [Sil] 


-0.334 


10.6 


728 ± 10 


1.3 


7.4 


1921 ± 46 


2.4 


3.6 


362 ± 29 


8.1 


6731 [Sil] 


-0.336 


8.4 


585 ± 22 


3.8 


5.2 


1405 ± 53 


3.8 


2.6 


270 ± 34 


12.6 


7065 Hei 


-0.377 


4.9 


293 ±6 


1.9 














7136 [Ariii] 


-0.385 


11.0 


705 ±9 


1.3 


4.0 


976 ± 52 


5.3 


6.9 


666 ± 23 


3.4 


7751 [Ariii] 


-0.451 


3.2 


196 ± 13 


6.4 














9069 [Siii] 


-0.561 


13.9 


1193 ± 111 


9.3 


3.1 


1007 ± 136 


13.5 


9.2 


1325 ± 90 


6.8 


9532 [Siii] 


-0.592 


21.3 


2510 ± 141 


5.6 


6.9 


2885 ± 185 


6.4 


13.1 


2488 ± 125 


5.0 


I(H/?)* 






3.16 X 10"'* 






0.74 X 10"'" 






2.40 X 10"'* 




c(H/3) 






0.13 ± 0.01 






0.36 ± 0.02 






0.07 ± 0.01 





" [Olí] AA 3726 4- 3729 
* ergseg"' cm~^ 
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Tabla 4.5: Líneas de emisión medidas en el espectro échelle en brote C obtenidas en la 
medida global, y en la componente angosta y en la ancha. Intensidad de la línea corregida 
por enrojecimiento [F(ií/3)=/(ií/3)=10000] 









global 






angosta 






ancha 




A(Á) 


f(A) 


-EW 


/(A) 


Error 


-EW 


/(A) 


Error 


-EW 


/(A) 


Error 






(A) 




(%) 


(A) 




(%) 


(A) 




(%) 


4861 H/? 


0.000 


3.2 


10000 ± 2188 


21.9 


1.0 


10000 ± 1464 


14.6 


1.1 


10000 ± 2349 


23.5 


4959 [Oiii] 


-0.024 


5.6 


6879 ± 391 


5.7 


3.4 


9818 ± 508 


5.2 


2.0 


5245 ± 612 


11.7 


5007 [Oiii] 


-0.035 


5.6 


20467 ± 1265 


6.2 


3.4 


29371 ± 1584 


5.4 


2.0 


15582 ± 1859 


11.9 


6563 Ha 


-0.313 


18.1 


28500 ± 1148 


4.0 


6.5 


28500 ± 709 


2.5 


8.0 


28500 ± 799 


2.8 


6584 [Nn] 


-0.316 


1.3 


2334 ± 740 


31.7 


0.9 


3852 ± 942 


24.4 


0.4 


1412 ± 588 


41.7 


6717 [Sil] 


-0.334 


2.6 


5015 ± 1591 


31.7 


1.2 


5400 ± 1357 


25.1 


1.4 


5316 ± 1890 


35.6 


6731 [Sil] 


-0.336 


1.8 


2891 ± 1071 


37.0 


0.7 


3186 ± 1295 


40.6 


1.0 


3891 ± 1606 


41.3 


I(H/3)* 






0.08x10-''* 






0.04 X 10"*'' 






0.04 X lO"'** 




c(H/3) 






0.21 ± 0.40 






0.14 ± 0.25 






0.38 ± 0.41 




* ergseg"^ 


cm^^ 





















Tabla 4.6: Líneas de emisión medidas en el espectro échelle en brote E obtenidas en la medida 
global, y en la componentes angostas. Intensidad de la línea corregida por enrojecimiento 
[F(ií/3)=/(ií/3)=10000] 









global 






angosta 1 






angosta 2 




A (A) 


f(A) 


-EW 

(A) 


/(A) 


Error 

(%) 


-EW 

(A) 


/(A) 


Error 


-EW 

(A) 


/(A) 


Error 

(%) 


3727 [Olí]" 


0.271 


8.2 


32110 ± 9745 


30.3 














4861 H/? 


0.000 


7.5 


10000 ± 1530 


15.3 


1.8 


10000 ± 1267 


12.7 


1.9 


10000 ± 1816 


18.2 


4959 [Oiii] 


-0.024 


2.9 


5747 ± 729 


12.7 


2.2 


3443 ± 351 


10.2 


4.4 


6900 ± 383 


5.6 


5007 [Oiii] 


-0.035 


6.9 


12160 ± 921 


7.6 


2.2 


10157 ± 1047 


10.3 


4.4 


20671 ± 1214 


5.9 


5876 Hei 


-0.209 


1.6 


1733 ± 561 


32.4 














6300 [01] 


-0.276 


2.1 


2017 ± 563 


27.9 














6548 [Nii] 


-0.311 


1.4 


1419 ± 468 


33.0 














6563 Ha 


-0.313 


32.0 


28500 ± 1087 


3.8 


15.9 


28500 ± 637 


2.2 


11.0 


28500 ± 964 


3.4 


6584 [Nii] 


-0.316 


5.8 


5268 ± 1157 


22.0 


1.7 


2810 ± 620 


22.0 


3.7 


9025 ± 2224 


24.6 


6717 [Sil] 


-0.334 


6.0 


6212 ± 1443 


23.2 


1.4 


2509 ± 512 


20.4 


4.3 


11233 ± 2878 


25.6 


6731 [Sil] 


-0.336 


5.2 


4448 ± 1400 


31.5 


0.8 


1875 ± 789 


42.1 


1.8 


6636 ± 2630 


39.6 


I(H/?)* 






0.18 X lO^"-* 






0.07 X 10"" 






0.07x10-'* 




c(H/3) 






0.00 






0.59 ± 0.22 






0.00 




" [Olí] AA 3726 + 3729 
* ergseg^^ cm^^ 
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Tabla 4.7: Líneas de emisión medidas en el espectro échelle en brote F obtenidas en la medida 
global. Intensidad de la línea [F{H /3) = I {H (3) =10000]. ^Las líneas no fueron corregidas por 
enrojecimiento, ver texto 









global 




A(Á) 


f(A) 


-EW 

(Á) 


F(A) 


Error 

(%) 


4861 H/3 


0.000 


-2.5 


10000 ± 3213 


32.1 


4959 [Olll] 


-0.024 


-1.2 


7274 ± 597 


8.2 


5007 [Oiii] 


-0.035 


-3.6 


21797 ± 1981 


9.1 


6563 Ha 


-0.313 


-10.6 


28500 ± 2593 


9.1 


6584 [Nii] 


-0.316 


-1.1 


3415 ± 2153 


63.0 


6717 [Sil] 


-0.334 


-2.6 


7317 ± 4024 


55.0 


6731 [Sil] 


-0.336 


-3.1 


9011 ± 5317 


59.0 


F(H/3)* 






0.04 X 10~i^ 




c(H/3) 






t 





erg seg ^ cm ^ 



tanto las intensidades de las líneas en este brote no fueron corregidas por enrojecimiento. 

4.4. Análisis de Abundancias 

La espectrocopía es una poderosa herramienta que permite la determinación de las abun- 
dancias de regiones Hii en el Universo Local haciendo uso de los métodos basados en la 
medición de las intensidades de las líneas de emisión y en la física atómica. Esto es bien co- 
nocido como el método "directo" . En el caso de las galaxias más distantes o intrínsecamente 
débiles, el bajo cociente S/N obtenido con los telescopios actuales se contrapone a la aplica- 
ción de este método y por ello es necesario la aplicación de métodos empíricos, los cuales se 
basan en las líneas de emisión más intensas. La base fundamental de estos métodos empíricos 
se conoce bastante bien (ver por ejemplo Pérez-Montero &: Díaz 2005). La exactitud de los 
resultados dependerá de la bondad de su calibración, que a su vez depende de un conjunto 
de abundancias derivadas con precisión a travéz del método "directo" para que el proceso de 
interpolación sea fiable. 

Sin embargo, derivar de forma precisa las abundancias de cada elemento no es un asunto 
sencillo. En primer lugar, porque es necesario medir las líneas de emisión con exactitud. Y en 
segundo lugar, porque se requiere de un cierto conocimiento de la estructura de ionización de la 
región con el fin de obtener las abundancias iónicas de los diferentes elementos y, en algunos 
casos, se necesita de los modelos de fotoionización para corregir los estados de ionización 
inobservables en el rango espectral en cuestión. Para ello es necesario cubrir un amplio rango 
espectral que incluya desde las líneas del ultravioleta, el doblete de [Olí] AA 3727,3729 A, hasta 

63 
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Tabla 4.8: Cocientes utilizados para derivar las densidades y temperaturas electrónicas 



Cociente de líneas 

ne[Sii] Rs2 = I(6717)/I(6731) 

Teioiii] Ros = I(4959,5007)/I(4363) 

Te [Olí] Ro2 = I(3727)/I(7319,7330) 

Teisiii] Rs3 = I(9069,9532)/I(6312) 

Te [Sil] R's2 = I(6717,6731)/I(4068,4076) 

Te[Nii] Rn2 = 1(6548,6584) /I(5755) 

el infrarrojo, con el doblete [Siii] AA 9069,9532 A. Para tener un diagnóstico preciso del medio 
se requiere de la medición de las líneas débiles aurórales. Y además, la precisión en la medición 
de los cocientes entre estas líneas aurórales y las líneas de recombinación de Balmer debe ser 
superior al 5 %, ya que en general, estas líneas son débiles, alrededor del 1 % de la intensidad 
H/3. Estos requerimientos permiten derivar las diferentes temperaturas de líneas: Te([Oiii]), 
Te ([Olí]), Te([Siii]), Te ([Sil]), Te([Nll]), necesarias para estudiar la estructura de temperatura 
y la ionización de cada brote o galaxia H II considerados como una región ionizada multizona. 
En esta sección serán explicados los diferentes métodos utilizados para estimar las densi- 
dades y temperaturas electrónicas con las que luego serán derivadas las abundancias iónicas y 
totales de las principales especies químicas presentes en las regiones Hii de Harol5. Las con- 
diciones físicas del gas como las abundancias químicas son derivadas y analizadas para cada 
brote de formación estelar (medida global) y para cada componente cinemática estudiada en 
el Capítulo 3. Se realizará todo este análisis siguiendo el procedimiento descripto en Hágele 
et al. (2008) y que se describirá a continuación. 

4.4.1. Propiedades físicas del gas ionizado 

Las propiedades físicas del gas ionizado, como las densidades y temperaturas electrónicas, 
pueden ser determinadas a partir de los cocientes de los flujos de las líneas observadas de 
algunos iones. Según Osterbrock (1989), la temperatura puede ser calculada a partir del 
cociente de transiciones colisionales con energía similar pero originadas en diferentes niveles. 
El ajuste entre este cociente y la temperatura electrónica se ha obtenido utilizando la tarea 
temden de IRAF. Dicha tarea está basada en un modelo de equilibrio estadístico de un átomo 
de 5 niveles (De Robertis et al., 1987; Shaw & Dufour, 1995). Las probabilidades de transición 
y niveles de energía son los que dicho programa trae por defecto. Como fuente de error de las 
densidades y temperaturas se han tomado las incertezas asociadas a la medición de los flujos 
de las línea de emisión y la corrección de enrojecimiento, y han sido propagados a través de 
nuestros cálculos. Los coeflcientes atómicos de colisión adoptados son los mismos que en el 
trabajo de Hágele et al. (2008) (ver Tabla 2.5 de dicho trabajo). En la Tabla 4.8 se muestran 
los cocientes de líneas de emisión necesarios para calcular las densidades y cada temperatura. 
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Densidades electrónicas 

El valor de la densidad electrónica, Ne, fue determinado a partir del cociente Rs2 de las 
líneas [Sil] AA6717, 6731 A, el cual es el más usado para determinar la densidad electrónica 
en la zona de baja excitación del gas ionizado. Aquí abajo se muestra la expresión para 
calcular las densidad a partir del mejor ajuste encontrado entre el cociente y la temperatura 
electrónica: 

„(i5„]) = io3. f^-°°¡;'+°'¡;> (4,3) 

Rs2 ■a2[t) + as{t) 

donde 

ao(í) = 2,21 - 1,3/í - l,25í + 0,23^^ 

ai(í) = -3,35 + 1,94/í + l,93í - 0,36*^ 

a2{t) = -4,33 + 2,33/í + 2,72í - 0,57*^ 

asit) = 1,84 - 1/í - l,14í + 0,24í2 

donde n([Sii]) = Ue = 10~'*Ne y t suele ser generalmente, te([Oiii]), donde te = 10~^Te, aunque 
se puede hacer un proceso iterativo para calcular con te ([Sil]), ya que esta temperatura, al 
igual que ocurre a con te ([Olí]), al ser un ion de tipo np^ depende de la densidad. En este 
caso el cociente que se debe utilizar es el K'g2- 

Las densidades derivadas en los brotes, tanto para las medidas globales (datos de ranura 
larga y del échelle), como para cada componente (datos échelle) pueden ser enmarcadas en 
el límite de baja densidad (ne<500cm~^), excepto en unos pocos casos donde la densidad 
presenta un error grande. Estos casos serán discutidos más adelante. 

En la Sección 4.5 se tratará cada caso en particular. 

Temperaturas electrónicas 

Dado que las temperaturas electrónicas son cantidades físicas importantes para la determi- 
nación de las abundancias químicas del gas ionizado, es necesario derivar su valor para poder 
obtener una estimación de la metalicidad del medio donde se forman las líneas de emisión. Sin 
embargo, en ocasiones, las líneas aurórales sensibles a la temperatura no son detectadas en el 
espectro de emisión de la nebulosa como ocurre, por ejemplo, en las regiones ricas en oxígeno, 
donde la línea auroral del [Oiii] A 4363 A es muy débil y difícil de detectar (Díaz ct al., 2007). 
En general, este hecho puede deberse a una mala relación señal-ruido del espectro o a una 
elevada metalicidad del gas que produce un enfriamiento tan eficaz que la temperatura es 
muy baja. En este último caso, las líneas no son detectables debido a que sus flujos son in- 
versamente proporcionales a la exponencial de la temperatura y son intrínsecamente débiles. 
Existen alternativas que consisten en la determinación de la abundancia del oxígeno haciendo 
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uso de las calibraciones empíricas con las líneas de emisión intensas fácilmente observables. En 
este trabajo se ha encontrado que varias de las regiones observadas en Haro 15 no presentan 
una o varias de las líneas aurórales necesarias para derivar las temperaturas con el método 
directo. Es por eso que a continuación se detallan separadamente las diferentes vías por las 
que uno puede optar para derivar las temperaturas, tanto por el método directo o recurriendo 
a los diferentes modelos teóricos o relaciones empíricas disponibles en la literatura, los cuales 
predicen relaciones entre temperaturas electrónicas y distintos parámetros empíricos y las 
temperaturas. 

La temperatura electrónia está definida como la temperatura de la distribución de Maxwell- 
Boltzmann asociada a las velocidades de los electrones que siguen dicha distribución. Entonces, 
por convención, cada vez que se escriba el término "temperatura" o "temperatura de línea" 
estaremos refiriéndonos al valor aproximado de la temperatura determinada por la energía 
cinética media de los electrones libres en la región donde se forma las líneas de emisión del 
ion estudiado. 

Método directo 

■ Oxígeno 

Para el cálculo de la temperatura de [Oiii] se usa el cociente de líneas Ros- En este caso, 
el ajuste obtenido, al igual que para el resto de iones de tipo np^, no depende de la densidad 
electrónica y es, por lo tanto, un buen indicador de temperatura. El ajuste entre cociente y 
temperatura es: 

47 77 
te{[Oni]) = 0,8254 - 0,0002415 • Ros + -^^ (4.4) 

Líos 

donde te es la temperatura en unidades de 10^ K 

En el caso del [Olí], la temperatura electrónica se calcula a partir del cociente Ro2- En 
este caso hay que tener algo de cuidado dado que las líneas aurórales [Olí] AA 7319, 7330 A 
suelen tener algo de contaminación por emisión producida por recombinación. Dicha emisión, 
no obstante, se puede cuantificar y corregir. Según Liu et al. (2000) dicha contribución puede 
ser ajustada en el rango 0.5 < t < 1.0 mediante la función: 

1^(7319 + 7330) _ 0,44 0'+ .4 .^ 

J^jj^^ -9,36.í --^ (4.5) 

expresión válida sólo en el rango de temperaturas entre los 5000 y los 10000 K. Además, el 
cociente de las líneas de [Olí] depende fuertemente de la densidad electrónica. La situación 
ideal es tener la densidad de [Olí] mediante el cociente I(3726Á)/I(3729A), pero generalmente 
no hay bastante resolución para resolver este doblete, por lo que suele ser necesario recurrir 
a la densidad de [Sil], también representativa de la zona de baja excitación. En este trabajo, 
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a pesar de poder resolver las líneas del [Olí] en el espectro échelle, ambas se encuentran en 
el borde del orden ultravioleta del échelle con lo cual se prefirió calcular la densidad con las 
líneas prohibidas del azufre dos veces ionizado. El ajuste obtenido es: 

OQ O 

íe([Oii]) = 0,23 + 0,0017 ■Ro2 + ^ + /iK) (4.6) 

iÍ02 

donde fi(ne) es una función dependiente de la densidad. 

■ Azufre 

La temperatura de [Sil] se calcula con el ajuste: 

teilSll]) = 1,92 - 0,0375 • R's2 " ^ + ^^ + f2Íne) (4.7) 

donde f2(ne) es una función dependiente de la densidad. 

En el caso de las líneas aurórales de [Sil] basta con poder medir una de las dos, ya que 
entre ellas hay una relación fija teórica, de modo que I(4068A) ~ 3 • I(4076A). De esta manera, 
es posible calcular la temperatura de [Sil] . 

En relación a la temperatura de [Siii] su medida directa se ha comenzado a utilizar cuando 
han empezado a ser accesibles las líneas colisionales en el infrarrojo cercano. El cociente de 
líneas que se usa en este caso es el R53. En ese caso, la temperatura de [Siii] se puede calcular 
mediante el ajuste: 

íe([5lll]) = ^^^ + ^^'^ (4.8) 

^^^ ^' 1,8- «53 -3,01 ^ ' 

■ Nitrógeno 

La temperatura de [Nii] puede ser calculada mediante el cociente de las líneas de nitrógeno, 
RAr2- Las líneas nebulares de [Nil] están muy próximas a Ha por lo que a veces aparecen 
fusionadas con ésta y no es posible medirlas o medir la más débil de ellas. En ese caso se suele 
recurrir a la relación teórica entre ellas, de tal manera que 1(6584) ~ 2.9 • 1(6548). Asimismo, 
la emisión de la línea auroral de [Nii] A 5755 A está afectada por procesos de recombinación, 
que puede ser corregida mediante la siguiente expresión propuesta por Liu et al. (2000): 

M^I^ = 319.í°-30.^ .49) 

en el rango de temperaturas entre 5000 y 20000 K. Hágele et al. (2008) usando el cálculo 
de temperatura electrónica [Olll], estimaron las contribuciones por recombinación tanto en 
las líneas aurórales del [Olí] como en la del [Nii] y encontraron que en todos los casos son 
menores al 4%. De esta forma si se opta por no corregir esta contribución, es necesario saber 
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que se cuenta con un pequeño error adicional. En estas condiciones, es posible determinar la 
temperatura de [Nii] mediante la expresión: 

/I O 1 or* 
te{[Nn]) = 0,537 + 0,000253 • Rn2 + -¡^ (4.10) 

que no depende de la densidad electrónica debido a la naturaleza np'^ del N+. 

En ninguno de los espectros obtenidos en este trabajo, y para ninguna de las regiones 
observadas fue posible medir la línea auroral del [Nii] en 5755A, por lo cual se tuvo que recurrir 
a modelos o relaciones entre temperaturas para poder derivar un valor de esta temperatura. 

Las expresiones anteriores son válidas en el rango de temperaturas electrónicas entre 7000 
y 23000 K y los errores involucrados en los ajustes son siempre menores que los errores 
observacionales por factores entre 5 y 10. Los errores en las temperaturas son los errores 
calculados a partir de las incertezas en las medidas de las intensidades de las líneas, aplicando 
la fórmula de propagación de errores, sin asignarle error alguno a la expresión del ajuste con 
el que se determina la temperatura electrónica. 

Relaciones entre temperaturas 

En ocasiones es posible medir sólo alguna de las líneas de emisión aurórales de los iones 
más representativos de la parte óptica del espectro. En esos casos, es necesario suponer una 
estructura interna de ionización de la región Hii de estudio. Las primeras aproximaciones, 
consideraban a la región H ii como isoterma, pero algunos estudios posteriores demostraron la 
existencia, en su interior, de una estructura interna de ionización que apunta a una variación 
espacial de la temperatura electrónica, además de la posible presencia de fluctuaciones de 
temperatura (Peimbert, 1967). Esto implica que, en la determinación de cada abundancia 
iónica, hay que tener en cuenta la temperatura media de la zona en la que el ion se encuentra. 
Una mejor aproximación es considerar dos zonas: una de alta ionización, la más interna y 
cercana al cúmulo ionizante, cuya temperatura característica sea la de [Oiii], y otra de baja 
ionización, más externa, caracterizada por la temperatura de [Olí]. Con la temperatura de la 
primera zona, se calculan las abundancias de iones como O^^, S^+ o Ne^^ y en la segunda 
las de O^, S^ o N"*" (Pagel et al., 1992). Los trabajos de Garnett (1992) refinaron un poco 
esta visión al demostrar por un lado, que la temperatura de [Slll] ocupa un valor intermedio 
entre ambas zonas y, por otro, que las temperaturas deducidas a partir de las líneas ópticas 
de emisión sobreestiman el valor de la teinperatura electrónica para toda la nube, ya que hay 
un efecto de selección hacia las zonas más calientes. 

Pérez-Montero & Díaz (2003; 2005) examinaron las relaciones entre las temperaturas que 
se utilizan para la determinación de la abundancia de cada ion, adaptándola a los resultados 
de los últimos modelos de fotoionización y la utilización de los coeficientes atómicos más 
recientes. Siguiendo Hágele et al., (2006; 2008), estas modificaciones han sido tenidas en 
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cuenta en el presente trabajo, y fueron utilizadas en la medida que eran requeridas para 
derivar temperaturas. Cabe mencionar que en las relaciones entre temperaturas electrónicas 
no es posible obtener una estimación realista de las incertidumbres asociadas a los modelos. 
Los errores de los modelos no se pueden cuantificar y por ello el error que aparecerá en las 
tablas será una cota inferior dada por las líneas de emisión. 

■ te[0lll] VS. te [Olí] 

En la mayoría de las regiones, las líneas aurórales de [Olí] no son observadas con suficiente 
relación señal-ruido como para derivar la temperatura utilizando el método directo. Por el 
contrario, las líneas aurórales del [Oiii] son generalmente "muy intensas" en este tipo de 
objetos de baja metalicidad y se las suele medir con una complicación menor pudiendo así de- 
terminar la Te([Oiii]) observacionalmente a través del método directo. En estos casos, se suele 
utilizar alguna relación basada en modelos de fotoionización para deducir Te ( [Olí] ) a partir de 
Te([Oiii]). Una relación que se usa comúnmente en la bibliografía y con bastante aceptación 
es la expresión dada por Pagel et al., (1992): 

teilOu]) = ,,^ ,^ , (4.11) 

'^^ ^' íe([OlIl])-l + 0,8 ^ ' 

basada en los modelos de Stasinska (1990). Sin embargo dicha relación obvia la dependen- 
cia de Te ([Olí]) con la densidad (Luridiana et al., 1999), hecho confirmado en parte por la 
gran dispersión observada al representar los objetos cuyas dos temperaturas han podido ser 
medidas, aunque no se ha encontrado una clara tendencia observacional con la densidad Ne 
(Hágele, 2008). El ajuste obtenido por Pérez-Montero & Díaz (2003) es: 

1,2 -F 0,002 -ne-h^^ 

íe([Oii]) = "^ ^ (4.12) 

^^ ^' í([Oiii])-i+ 0,08 + 0,003 -ne + lf 

y es la que se usa en este trabajo cuando la situación lo requiera, y se tenga una estimación de 
Te[Oiii], y siempre teniendo en cuenta que en esta relación Te ([Olí]) depende explícitamente 
de la densidad electrónica. 

■ te[SlIl] VS. te[Oin] 

Existen trabajos que relacionan la temperatura de [Oiii] con la temperatura de [Siii]. La 
motivación de dicha relación aparece con Garnett (1992) quien sugirió que la Te([Siii]) ocupa 
un valor medio entre la de [Oiii] y la de [Olí], y calculó la abundancia de S^"*" sólo con la línea 
auroral [Siii] A 6312 A. La relación que da Garnett es: 

íe([5'iii]) = 0,83- íe([Oin])+ 0,17 (4.13) 
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En los últimos años, la importancia de esta relación ha tomado más relevancia, dado que en 
objetos de muy alta metalicidad se ha detectado con mayor facilidad la línea auroral de [Siii] 
que las líneas aurórales de otros iones y elementos. Sin embargo, la aparición de coeficientes 
atómicos más nuevos para las transiciones del azufre (Tayal &: Gupta, 1999) hizo que la 
relación subestime la Te([Slll]). En su lugar, con los nuevos coeficientes de los modelos de 
fotoionización descritos en Pérez-Montero &: Díaz (2005) la relación tomó la forma: 

íe([5iii]) = 1,05 • íe([Oin]) - 0,08 (4.14) 

Dado que la relación entre temperaturas electrónicas no permite obtener una estimación 
realista de las incertidumbres asociadas a los modelos y con la aparición de un cada vez mayor 
número de observaciones en el rango espectral donde se puede observar esta temperatura, ha 
sido posible obtener un nuevo ajuste empírico que tiene la forma: 

teiiSin]) = (1,19 ± 0,08)íe([Oiii]) - (0,32 ± 0,10) (4.15) 

definida por Hágele et al. (2006), de ahora en más H06. Este última relación empírica es la 
que se ha optado por aplicar en este trabajo de Tesis cuando la situación lo requiera. 

■ Análisis clásico: te[Oii] vs. te [Sil] vs. te[Nii] 

En general, es comúnmente aceptado que los tres iones. O"*", S"*" y N"*", se forman en la zona de 
baja excitación de una nebulosa y, por lo tanto, deberían tener temperaturas muy similares: 
Te ([Olí]) « Te ([Sil]) sá Te([Nii]). Dado que es en las zonas más externas de la nebulosa donde 
se producen las variaciones locales más significativas de la temperatura, ha sido necesario 
profundizar en las relaciones entre temperaturas para mejorar las determinaciones de las 
abundancias de estos tres iones, teniendo siempre en consideración la dependencia de Te ([Olí]) 
con la densidad. 

En aquellos casos en que las líneas aurórales de [Sil] no están accesibles en el espectro, 
se suele suponer que Te([Sii]) wTe([Oii]). Sin embargo, existen evidencias que apuntan a que 
esta Te ([Sil]) es levemente menor a la Te ([Olí]). De los modelos utilizados por Pérez-Montero 
(fe Díaz (2003), se obtiene un ajuste lineal que resulta ser: 

tei[Su]) = 0,71 • íe([Oii]) + 0,12 (4.16) 

para una densidad electrónica de 10 partículas por cm^. Para densidades electrónicas de 100 
partículas por cm'^ el ajuste da: 

tei[Su]) = 0,86- íe([On]) + 0,08 (4.17) 

Las observaciones confirman que en general la temperatura electrónica de [Sil] es menor 
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a la de [Oil] para la mayor parte de las galaxias Hll , sin embargo, muchas de las galaxias 
Hii presentan un valor de Te([Sii]) menor que el predicho por los modelos. Esto podría 
deberse a la contribución del gas difuso ionizado por un posible escape de fotones o a una 
mayor contribución del [Sil] en las zonas de baja excitación (Pérez-Montero & Díaz, 2003). A 
pesar de estas evidencias, no existe una certeza de dicha dependencia desde el punto de vista 
observacional (ver Figura 2.17 de Hágele (2008). Luego, la Te([Sii]) se ha obtenido a partir 
de la primera relación arriba mencionada. 

Por otro lado, como se mencionó más arriba, la línea [Nii] A 5755 A lamentablemente suele 
tener muy poca señal-ruido o es inobservable en nuestros espectros, por lo que es habitual 
considerar la aproximación Te([Nii]) ~Te([Oii]) como válida. Esta relación, confirmada por 
los modelos, es sensible a la densidad y a la estructura de ionización interna de la nebulosa, 
pudiendo alcanzar valores más cercanos a Te([Sii]), por lo que, en todo caso, debe ser tomada 
con su grado de incerteza. De los modelos de fotoionización propuestos por Pérez- Montero & 
Díaz (2003) resulta que la Te([Sii]) es sensiblemente menor que Te([Oli]) y que Te([Nii]), lo 
cual puede deberse, según los autores, a la contribución de las zonas de baja excitación. Este 
efecto puede no ser importante en objetos de alta excitación, en que las fracciones de S^ son 
menores, pero sí puede ser relevante a la hora de determinar la abundancia total de azufre en 
objetos menos excitados (ver discusión en Pérez- Montero et al., 2010). 

Con las relaciones hasta aquí expuestas, se puede determinar las temperaturas faltantes 
debido a la no detección de alguna de las líneas aurórales correspondientes. En los casos en 
que ninguna línea auroral sea detectable, y con sólo la medida de las líneas más intensas, se 
puede recurrir al uso de los dos métodos empíricos abajo descriptos. 

■ Método empírico D07 

El método utilizado para derivar la temperatura de [Siii] cuando no se tienen las líneas auróra- 
les de este mismo elemento ni las de ningún otro ion o elemento, es el desarrollado por Díaz et 
al. (2007) basado en el parámetro empírico S023^ definido por Díaz ¿¿ Pérez-Montero (2000) 
y recientemente reivindicado como un buen indicador de metalicidad para regiones Hii de 
alta metalicidad. Dicho parámetro es la combinación de los parámetros empíricos 023^ (Pagel 
et al., 1979) y del 823^ (Vilchez &: Esteban, 1996) usados para la determinación de las abun- 
dancias de oxígeno, de los cuales se hablará en la Sección 4.5.4. La relación empírica entre la 
temperatura y el parámetro empírico esta dada por el siguiente ajuste: 

íe([5'iii]) = 0,596- 0,283/0^5023 + 0,199(logS'O23)^ (4.18) 

Esta relación es válida para regiones de alta metalicidad. Pero de la Figura 9 de Díaz et 
al., (2007) se desprende que esta relación es válida también para el extremo de la zona donde 



2 
SO23 — S23/O23 

^023 = 1(3727, 29A) + /(4959A) + /(5007A)//{ií^) 

''S23 = -r(6717A) + i(673lA) + I(9069A) + /{9532Á)//(ír,9) 
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se encuentran las galaxias Hii, donde la dispersión es baja, y que es el rango de valores del 
parámetro empírico SO23 donde se encuentra Haro 15. 

Para poder estimar así la Te([Siii]), es necesaria la detección y medida de las líneas 
intensas del oxígeno y del azufre. Luego, haciendo uso de la relación empírica de H06 se 
puede estimar la temperatura del [Oiii] y a partir de esta temperatura la de [Olí] a través 
de los modelos, y posteriormente la de [Sil] y la de [Nii] igualándolas a Te ([Olí]). El método 
para derivar las abundancias con esta temperatura y usando las líneas intensas se denomina 
método semi-empírico. 

■ Método empírico P07 

En ocasiones, las líneas intensas del [Siii] no son observadas con buena calidad, o están dema- 
siado hacia el rojo para el rango espectral obtenido. Pilyugin (2007) propuso una calibración 
para la relación de intensidades entre las líneas intensas y la auroral del [Nii] en términos de 
las líneas del oxígeno, [Olí] AA3727, 3729 Á y [Oiii] AA4959, 5007 Á, y de la línea de recom- 
binación H/3. Entonces, la temperatura electrónica del [Nii], que aparentemente caracteriza a 
la zona de baja ionización de la nebulosa, puede ser derivada (ver Pilyugin, 2007). 

1,111 

Í9 = f4 19) 

log Qmi - 0,892 - 0,144 • log Í2 + 0,023 • Í2 

donde t2 es la temperatura representativa de la zona de baja excitación y Qath es la aproxi- 
mación al cociente de líneas de [Nii] y sale de la relación empírica 

log Qjvii = 2,619 - 0,609 • log R2 - 0,010 • [log R2f + 1,085 • log(l - P) + 0,382 • [log(l - P)]^ 

donde P es el parámetro de excitación que se expresa como 

P = Rs/iR^ + ^2) 

con R2 y R3 definidos como 

R2 = /[Olí] AA 3727,3729/1(^/3) R3 = /[Oiii] AA4959, 5007/(i//3) 

Siguiendo el método descripto por Pilyugin, Díaz et al., (2007) encuentran para las regiones 
circumnucleares de formación estelar (CNSFR) una correlación entre las temperaturas del 
[Siii] y del [Nii], las cuales son muy similares, siendo en promedio Te([Nii]) 500K mayor que 
Te([Slli]) (ver Figura 15 de Díaz et al., 2007). Una vez estimada la temperatura del nitrógeno, 
es posible aproximar un valor a la temperatura del [Siii] con la aproximación Te([Nii]) ~ 500K 
+ Te([Siii]). Y así, una vez más haciendo uso de la relación empírica de H06, se puede estimar 
la temperatura del [Oiii]. Luego, con esta temperatura y los modelos se determina Te([Oii]) y 
posteriormente la de Te ([Sil]) suponiendo que es igual a la temperatura electrónica de [Olí]. 
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- T=10^K 

Cuando no sólo no se pueden medir las líneas aurórales, sino que además no se cuenta 
con algunas de las líneas intensas para el cálculo de temperaturas electrónicas a través de los 
métodos empíricos que se acaban de describir, y considerando que para esta clase de objetos 
la temperatura suele estar alrededor de 10000 K, se tomará este valor para la temperatura 
electrónica del [Oiii], y con esta temperatura y los modelos se determinará Te([Oii]), y poste- 
riormente Te([Sii]) y Te([Nii]), suponiendo ambas igual a la temperatura electrónica de [Olí], 
y de Te([Siii]) se determinará a través de su relación con Te([Oiii]) derivada por H06. 

4.4.2. Derivación de las abundancias químicas 

Para estudiar las abundancias químicas globales y discriminadas por componentes ci- 
nemáticas en cada brote de formación estelar de Harol5, se han derivado las abundancias 
de las diferentes especies iónicas usando las líneas de emisión más intensas disponibles y de- 
tectadas en los espectros analizados. Para ello se utilizó la tarea ionic del paquete STSDAS 
de IRAF, tal como se describe en Hágele et al. (2008). Este paquete también se basa en la 
aproximación de un átomo de cinco niveles en equilibrio estadístico (De Robertis et al., 1987; 
Shaw & Dufour, 1995). Uno de los problemas para la determinación de abundancias es no 
poder medir las líneas de emisión de todos los estados de ionización de los elementos cuya 
abundancia se quiere calcular. En general, una manera de solucionar este problema es recu- 
rrir a la cantidad llamada factor de corrección de ionización (ICE sigla en inglés de ionization 
correction factors) que da cuenta de la contribución a la abundancia de las especies iónicas 
que no han sido observadas. En este Capítulo las abundancias totales han sido obtenidas 
teniendo en cuenta, cuando fuese necesario, los estados de ionización inobservables de cada 
elemento, recurriendo a los factores de corrección de ionización ampliamente aceptados para 
cada especie, X/H = ICE(X+*) X+VH+ (Pérez-Montero et al., 2007; Hágele et al, 2008). 

A continuación se describen las expresiones de los ajustes usados para cada ion, una vez 
que se dispone de las densidades y temperaturas electrónicas derivadas adecuadamente para 
cada componente y medida global en cada uno de los brotes estudiados. 

Helio 

Las líneas de helio, al igual que las de hidrógeno son líneas de recombinación. Son bas- 
tante intensas y numerosas, aunque muchas de ellas suelen estar fusionadas con otras líneas. 
En general, están afectadas por la absorción de la población estelar subyacente, por fluores- 
cencia, y tienen una contribución por excitación colisional (para un completo tratamiento 
de estas contribuciones ver Olive & Skillman, 2001, 2004). Generalmente, para estimar las 
abundancias de helio (una y dos veces ionizado) se usan las líneas de HeiAA4471, 5876, 6678 
y 7065A, y HeiiA4686Á, respectivamente. En este trabajo se ha tomado la temperatura 
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electrónica Te([Olii]) como representativa de la zona en que estas líneas son emitidas. Usando 
las ecuaciones dadas por Olive & Skillman se puede derivar el valor de He"'"/H+, usando la 
emisividad teórica escalada a H/3 de Benjamín (1999) y las expresiones para los factores de 
corrección colisional de Kingdon & Ferland (1995). Para calcular la abundancia de helio dos 
veces ionizado, He^^/H"*", se ha utilizado la ecuación (9) de Kunth & Sargent (1983). Dado 
que los objetos observados tienen densidades electrónicas bajas, la leve dependencia con la 
profundidad óptica que tienen tres de las líneas medidas del He será ignorada y por eso no se 
ha hecho ninguna corrección por fluorescencia. Tampoco se ha hecho ninguna corrección por 
la posible presencia de alguna población estelar subyacente dado que las líneas de absorción 
de este elemento son más angostas y no se pueden ver sus alas a los lados de la emisión. 
Con todo esto, es válida la aproximación para la abundancia total del helio: 

He ^ He +Me ^^ ^0) 

Según sea el caso, los resultados obtenidos para cada línea del helio junto a sus corres- 
pondientes errores son presentados en las Tablas 4.10 y 4.11. Se adjunta el valor "adoptado" 
para He"'" / H"*" que es el promedio ponderado por los errores de las diferentes abundancias 
iónicas derivadas para cada línea de emisión de He I. Este valor es conocido como "adop- 
ted" (del inglés, valor considerado o adoptado). De aquí en adelante se utilizará el término 
valor "adoptado" teniendo en cuenta esta definición. 

Abundancias químicas iónicas y totales determinadas con las líneas prohibidas 

Se han calculado las abundancias iónicas y totales de O, S, N, Ne, Ar y Fe, cuando las 
líneas están presentes, utilizando las expresiones dadas por Hágele et al., (2008) y que son 
incluidas en la Tabla 4.9. Muchas veces suele ser útil expresar el cociente de algunos elemen- 
tos respecto del oxígeno, por ejemplo, N/0, S/0, Ne/0, Ar/0, en unidades logarítmicas. 
Esto tiene la ventaja que como las abundancias iónicas se han determinado a partir de dos 
líneas de excitación colisional (por ejemplo, N"^/0"^), estos cocientes no dependen tanto de 
la temperatura electrónica. En los casos en los que no existen mediciones, ni de la línea de 
[Sil] A 4068 A ni de la línea de [Nii] A 5755 A (esta línea no se ha podido medir en ninguno 
de los brotes), y bajo la suposición de una temperatura electrónica homogénea en la zona de 
baja excitación, se ha considerado la aproximación Te([Nii]) ~ Te([Sii]) ~ Te([Oii]) explicada 
anteriormente en el análisis clásico. 

(i) Los cocientes de abundancia iónica del oxígeno, O^^/H"*" y 0^"'"/H"'", se pueden determinar 
a partir de las intensidades de las líneas [Olí] AA 3727,3729 Á y [Oiii] AA 4959,5007 Á, 
respectivamente, usando para cada ion su correspondiente temperatura electrónica. Sin 
considerar la contribución de oxígeno neutro que, en proporción, es igual a la fracción 
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Tabla 4.9: Ecuaciones para el cálculo de las abundancias de las líneas prohibidas (Hágele et al., 2008). 

Ion* Ecuación 
"0+ log (^f^) + 5,992 + ^- 0,681 • logí^ + log(l + 2,3 • Ue) 
02+ log r»)+l(5007)\ + 6^i¡4 + M51 _ 0^55 . ^^^^^ 

g+ log N6n7)+/(673i) j ^ 5^423 + °f5 - 0,28 • logíe + 1,0 • n,) 
S2+ log / W+g9532)< ^ 5^g ^ o¿7i^_ 0^22 . logí, 

N+ log (^^^^flf^) + 6,273 + % - 0,592 • logí^ 

Ne2+ log ^^^j+ 6,486 +if^- 0,504- log íe 

Ar2+ log ^^|)j+ 6,157+ ^-0,508- log íe 

Ar3+ log ÍW\+ 5^705 + i¿^- 0,156 -logí, 

Fe2+ log (M) + 3,504 + if^ - 0,483 • logí^ 

* Como 12+log(X/H+) 

te es la temperatura electrónica de la línea, en unidades de 10*^ K 

Ug es la densidad electrónica, en unidades de 10^ cm~^ 

de hidrógeno neutro, gracias a la reacción de intercambio de carga (Osterbrock, 1989): 

0+ + ií° ^ 0° + H+ 

Dado que a la temperatura a la que se encuentran los objetos aquí observados, la mayor 
parte del oxígeno se encuentra en forma de O"^ y O'^^, la abundancia total de oxígeno 
es bien representada por la expresión: 

(ii) Las abundancias iónicas del azufre, S^/H+ y S^^/H"^, están bien representadas por las 
líneas [Sil] AA 6717, 6731 Á y las líneas [Siii] AA9069, 9532Á, respectivamente. Depen- 
diendo del grado de excitación de la región, se podría esperar una contribución relativa- 
mente importante de S^'^. Para esos estados de ionización inobservables, la abundancia 
total de azufre se puede calcular usando un ICF para S'*"+S^"'". Una buena aproxi- 
mación es la fórmula de Barker (1980), basada en los modelos de fotoionización de 
Stasinska (1978) y toma la forma: 



ICF{S+ + S 



2+^ 



O 



2+ 



0+ + 02H 



-l/o 



(4.22) 



con aw2.5, medida para una muestra de objetos con la línea de [Siv] a 10.5 /u por 
Pérez-Montero et al., (2006). Utilizando este ICE como una función del cociente 0^"'"/0 
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en lugar del cociente O^/O se reduce el error propagado por esta cantidad, ya que los 
errores asociados a 0^+ son considerablemente menores que O"*" (Hágele et al., 2008). 

(iii) La abundancia iónica del nitrógeno, N+/H^, se deriva a partir de la intensidad de las 
líneas [Nii] AA 6548, 6584 A. En caso de no disponer de alguna de ellas por su proximidad 
a la línea Ha, se recurre a la relación teórica entre ambas líneas del [Nii] antes citada. 
Para el cálculo de la abundancia total de nitrógeno se parte de la relación: 

N N+ , , 

o = OT ("-23) 

luego la abundancia total del nitrógeno, N/H, es calculada como: 

log{N/H) = log{N/0) + log{0/H) (4.24) 

(iv) La abundancia iónica del neón, Ne^"*", se deriva usando la línea [Neiii] A 3868 A. Dado 
que su estructura de ionización es bastante similar a la del oxígeno, clásicamente la 
abundancia total del neón ha sido calculada asumiendo: 

^e Ne^+ ,^ ^^, 

-o=ly^ (4-25) 

La temperatura que se considera para este ion es la representativa de la zona de alta 
ionización Te([Neiil]) ?aTe([Oiii]) (Pcimbert & Costero, 1969). Izotov et al. (2004) in- 
dicaron que el hacer esta suposición puede llevar a una sobreestimación de Ne/H en 
objetos con baja excitación, donde la transferencia de carga entre O^"*" y H'' se vuel- 
ve importante. Por este motivo aquí se calcula con el ICF dado por Pérez-Montero et 
al. (2007). Dada la alta excitación de los objetos observados, no hay diferencias signifi- 
cativas entre la abundancia total derivada para el neón usando este ICF y la estimada 
con el modelo clásico (Hágele et al., 2008, 2011). 

(v) La abundancia iónica del argón, Ar^"*", se calcula usando la línea [Ariii] A 7136 A asumien- 
do que Te([Arlii]) « Te([Siii]) (Garnctt, 1992). También es posible medir, en algunos es- 
pectros, las líneas del [ArivJA 4713 y 4740 A. En ese caso la abundancia iónica del Ar^^, 
se calcula asumiendo que Te([Ariv]) «Te([Oiii]). No obstante, la primera de las líneas 
del [Ariv] suele aparecer fusionada con otra línea de HeIA4711 A, difícil de corregir, por 
lo que es mejor usar la segunda, más intensa, para calcular la abundancia de Ar'^''". La 
abundancia total del argón es calculada usando el ICF(Ar^+) y el ICF(Ar^+ + Ar'^"*') 
derivados a partir de modelos de fotoionización por Pérez-Montero et al. (2007). 

(vi) La abundancia iónica del hierro se calcula usando la línea [Feíii] A 4658 A y la tempe- 
ratura asumida es la del [Oiii]. La abundancia total es posible calcularla usando el 
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ICF(Fe^"'") de Rodríguez &: Rubín (2004). En este trabajo únicamente fue posible medir 
la línea [Feíii] A 4658 A en el espectro de ranura larga del brote B. 

Las abundancias iónicas de los elementos más pesados que el helio, los factores de co- 
rrección de ionización ICFs, las abundancias totales con sus correspondientes errores son 
mostrados en la Tabla 4.10 para los datos de ranura larga, y 4.11 para los datos échelle. 

4.5. Discusión 

4.5.1. Datos de ranura larga 

Cuatro temperaturas electrónicas, Te([Oiii]), Te([Oii]), Te([Siii]) y Te([Sii]), han sido de- 
rivadas de mediciones directas en el brote B. La buena calidad de los datos permitió alcanzar 
precisiones del orden de 2 %, 3 %, 10 % y 14 %, respectivamente para cada una de las tempera- 
turas. La obtención de las cuatro temperaturas fue gracias a la buena relación S/N del espectro 
que permitió la detección de las líneas débiles aurórales, como [Oiii] A 4363 A, [Sil] A 4068 A, 
[Slii] A 6312 Á y [Olí] AA 7319, 7330 Á. Dado que la línea auroral [Nii] A 5755 A no ha sido po- 
sible medirla en el espectro, se ha considerado la aproximación Te([Nii]) ~Te([Oii]) para la 
determinación de la temperatura Te([Nii]). 

Se ha calculado la densidad electrónica a partir del cociente de líneas del [Sil] , obteniendo 
un valor de 100 partículas por cm^, aunque es necesario destacar que el error de la misma es 
muy grande, por lo cual este valor sirve para determinar el orden de magnitud de la densidad. 

Las abundancias iónicas del O"*", 0^+, S+, S^^, N+, Ne^"*", Ar^"*" y Ar^"*" se han deducido 
utilizando las expresiones dadas en la Sección 4.4.2, las cuales dependen de la intensidad de 
las líneas intensas del ion cuya abundancia se quiere calcular, así como de la densidad (en 
algunos casos) y temperaturas electrónicas. Como en este brote fue posible medir la línea del 
hierro 4658A se obtuvo también Fe^"*" considerando la temperatura del [Oiii] para su cálculo. 
Además, fueron calculadas las abundancias totales del O, S, N, Ne, Ar y Fe, y los cocientes 
logarítmicos N/0, S/0, Ne/0 y Ar/0. Tanto las abundancias iónicas como las abundancias 
totales y los cocientes logarítmicos respecto al oxígeno, junto con sus respectivos errores, están 
tabulados en la Tabla 4.10. 

El brote B tiene una abundancia de oxígeno 0.22 veces el valor solar (Allende Prieto et al., 
2001, 12+log(0/H) = 8.69). Este valor de abundancia (8.04±0.03) es similar, dentro de los 
errores, al encontrado por López-Sánchez & Esteban (2009b) que es igual a 8.10zb0.06. La 
abundancia de azufre es 0.15 veces el solar (Grevesse & Sauval, 1998, 12+log(S/H) = 7.33), 
y la abundancia de nitrógeno es 0.09 veces el solar (Holweger, 2001, 12+log(N/H) = 7.93). 
Respecto a las abundancias relativas (ver Figura 4.4), se estimó un valor de log(N/0) = - 
1.16 el cual es 0.28 dex menor al valor solar (Asplund ct al., 2005, log(N/0) = -0.88) y 
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Tabla 4.10: Propiedades físicas, abundancias químicas iónicas y totales derivadas de las líneas 
prohibidas y de la línea de recombinación del helio para la medida global en los datos de 
ranura larga: brotes B y C 



n([Sll]) 


100: 


100: 


Te([Oni]) 


1.26±0.03 


1.01±0.02'^ 


Te([SlIl]) 


1.19±0.12 


0.88±0.13''í' 


T.([Oii]) 


1.20±0.04 


1.21±0.01" 


Te ([Su]) 


0.79±0.11 


1.21±0.01" 


T.([Nn]) 


1.20±0.04" 


0.93±0.01í' 


12+log(0+/íí+) 


7.34±0.05 


7.82±0.06" 


12+\og(0^+/H+) 


7.95±0.03 


7.84±0.03'^ 


12+log(0/H) 


8.04±0.03 


8.13±0.04* 


12+log{S+/H+) 


5.95±0.20 


6.16±0.07" 


12+log(52+//í+) 


6.17±0.23 




ICF(S'+ + 52+) 


1.41±0.03 




12+log(S/H) 


6.52±0.22 




log(S/0) 


-1.52±0.22 





12+\og{N+/H+) 

12+log(N/H) 

log(N/0) 



6.18±0.04™ 
6.88±0.20"' 
-1.16±0.07" 



6.70±0.07í' 
7.02±0.23í' 
-1.11±0.09P 



12+log(iVe^+/í/+) 


7.33±0.03 




ICF(Afe2+) 


1.08±0.01 




12+log(Ne/H) 


7.36±0.03 




log(Ne/0) 


-0.68±0.04 




12+log(Ar^+/_fJ'+) 


5.65±0.10 




l2+log{Ar^+/H+) 


4.62±0.08 




lCF{Ar^++Ar^+) 


1.03±0.01 




12+log(Ar/H) 


5.70±0.10 




log(Ar/0) 


-2.34±0.11 




12+log(Fe^+/_ff+) 


5.04±0.10 




ICF(Fe2+) 


5.21±0.62 




12+log(Fe/H) 


5.76±0.11 




He+/H+{\4'Í71) 


0.064±0.003 




He+/H+{X 5876) 


0.096±0.002 




He+/H+ (\6678) 


0.072±0.003 




He+/H+(\ 7065) 


0.100±0.008 




ífe+/ií+adoptcd 


0.086±0.015 




He'^+/H+{\A686) 


0.0013±0.0002 




(He/H) 


0.087±0.015 





densidades en cm '' y temperaturas en 10* K 

"* resultados usando las temperaturas predichas por los modelos, eomo se describe en el texto. '^ resultados usando las 
temperaturas derivadas con el método empírico D07. '' resultados usando las temperaturas estimadas con el método 
empírico encontrado por Hágele et al. (2006). ^ resultados usando las temperaturas derivadas con el método empírico 
P07. * resultados usando los datos necesarios para su cálculo que fueron derivados de modelos o relaciones empíricas. 
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0.44 dex mayor al valor típico para este tipo de objetos (el valor del plateau se encuentra en 
Iog(N/0) =-1.6, ver Amorín ct al., 2010; Pérez-Montero et al., 2011) y un valor de log(S/0) = 
-1.52±0.22 que está de acuerdo, dentro de los errores, con el valor solar (Grevesse &¿ Sauval, 
1998, log(S/O)0=-1.36). 

Ha sido posible estimar la abundancia total de He a partir de la obtención de las abundan- 
cias iónicas He"*" y He^"*" dado que fue posible medir las líneas, HelAA4471, 5876, 6678 y 
7065Á, y HeiiA4686Á respectivamente. El valor "adoptado" para He+/H+ es 0.086±0.015 y 
0.0013ib0.0002 para el He'^"^/H+ dando un valor de la abundancia total del He estimada en 
esta región de 0.087ib0.015, valores típicos para las galaxias H ii (Hágele et al., 2008). Los 
resultados pueden verse en la misma Tabla antes mencionada. 

En el brote C no fue posible obtener de forma directa las temperaturas electrónicas, dado 
que no se pudo medir las líneas aurórales correspondientes para sus cálculos. En este caso, se 
recurrió tanto a los modelos como a las relaciones empíricas para determinar las diferentes 
temperaturas. Con el método de P07 descripto en la Sección 4.4.1, se estimó la temperatura 
del [Nii], encontrando un valor para la misma de 0.93ib0.01. Y haciendo uso de la aproximación 
Te([Nii])ss500K + Te([Siii]) dada por D07 se estimó un valor para la Te([Siii]), con el que 
se derivó Te([Oiii]) a partir de la relación de H06. La temperatura electrónica Te ([Olí]) fue 
calculada a partir de Te([Olli]) utilizando los modelos de Pérez-Montero & Díaz (2003) y de 
la aproximación Te([Sii]) ~ Te([Oii]) se obtuvo la temperatura de [Sil]. También en este brote 
se estimó una densidad de 100 partículas por cm^ con errores muy grandes. Los resultados 
de las abundancias iónicas y totales usando las temperaturas estimadas con los modelos y los 
métodos empíricos, y los ICFs son expuestos en la Tabla 4.10. En este brote no se ha medido 
la línea [Feíii] A4658 A, por lo que no se pudo derivar una abundancia de este elemento. 
Tampoco, como se ha comentado en la Sección 4.3, se ha podido medir la línea [Nii] A 6548 A, 
con lo cual para el cálculo de las abundancias del nitrógeno, se ha recurrido a su relación 
teórica con la línea [Nil] A 6584 A. Para este brote se encontró un valor de la abundancia de 
oxígeno de 8.13±0.04, que es 0.28 veces el solar. 

Comparando los resultados de los datos de ambos brotes en los espectros de ranura larga, 
se puede ver que la abundancia de oxígeno es mayor en el brote C que en el B pero ambas 
por debajo del valor solar, con valores típicos en esta clase de objetos (ver por ejemplo, 
Terlevich et al., 1991; Hoyos & Díaz, 2006). Las temperaturas representativas de la zona de alta 
ionización son mayores en la región B. Las temperaturas correspondientes a la zona de baja 
ionización, en el caso de Te ([Olí]) aparentan ser similares en ambos brotes, mientras que para 
Te ([Sil]) y Te([Nii]) las temperaturas fueron derivadas utilizando distintas aproximaciones. 
En el brote B, la temperatura Te ([Sil]), obtenida de las medidas directas y es 4200 K mayor 
que en la del brote C, donde se derivó de suponer la igualdad con Te ([Olí]). Es de notar que en 
el brote B, Te ([Olí]) es 4100 K más alta que Te([Sii]). En cuanto a la temperatura de [Nii] en 
el brote B, el valor obtenido es mayor que la del brote C en 2700 K, pero nuevamente fueron 
estimados utilizando dos aproximaciones distintas. A modo de comparación se ha derivado la 
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temperatura del [Nii] para el brote B usando el método P07. El valor encontrado para B es 
Te([Nii]) = 10500±100 K, que comparado con el valor de Te([Nii]) encontrado para el brote 
C usando el mismo método P07 es ~1200 K mayor que éste. El valor estimado para N/0 es 
un poco mayor en el brote C que en el B, aunque igual dentro de los errores (ver Figura 4.4). 
El grado de excitación, log(0^"'"/0"'"), calculado para el brote B es 0.61 y para C es 0.02. Sin 
embargo, debe tenerse en cuenta que, para el brote C las abundancias iónicas se han deducido 
utilizando las temperaturas derivadas de los modelos y relaciones empíricas. Sin embargo, 
igualmente se puede considerar que es mayor el grado de excitación en el brote B que en el 
C, concordando este resultado con lo estimado por López-Sánchez & Esteban (2009a). La 
extinción observada en ambos brotes B y C es similar dentro de los errores y consistente con 
baja extinción. 

4.5.2. Datos échelle 

En esta subsección se analizarán los resultados alcanzados para las regiones estudiadas 
con los datos de los espectros échelle. Cada brote será analizado por separado, estudiando 
las temperaturas, densidades, abundancias iónicas y totales estimadas para la medida global 
y para cada componente cinemática separadamente. Luego se dará un análisis comparativo 
entre brotes. 

Haro 15 A 

El análisis que a continuación se detalla, se refiere tanto a la medida global como a las 
dos componentes angostas y la componente ancha del brote A, salvo que se especifique algún 
análisis particular realizado sobre algunas de las componentes o medida global. Las líneas 
aurórales [Olll] A 4363 Á, [Sil] A 4068 Á, [Siii] A 6312 Á, [Olí] AA 7319, 7330 Á y [Nii] A 5755 Á no 
se han podido detectar en el espectro de la región central de Haro 15, el brote A. En este caso, 
utilizando la relación entre temperaturas Te[Siii] vs. Te[Olll] dada por el ajuste encontrado 
en H06 se ha podido estimar la temperatura electrónica del [Oiii], previa determinación de 
Te([Slli]) con el método D07 usando el parámetro empírico S023. Luego, siguiendo el análisis 
clásico, se han podido obtener de los modelos, las temperaturas Te([Oii]), Te([Sii]) y Te([Nii]). 

Las densidades estimadas en los cuatro casos aquí analizados, están por debajo de la 
densidad crítica por desexcitación colisional. La medida global tiene una densidad de 300 
partículas por cm'^; la componente angosta 2, 200 partículas por cm^; y tanto la componente 
angosta 1 como la componente ancha, 100 partículas por cm'^. Aún teniendo en cuenta que 
los valores de la densidad son una estimación de su magnitud, dado que los errores son muy 
grandes, parecería que la componente angosta 2 tiene una mayor densidad que las otras dos 
componentes. Sin embargo, y dado que la medida global es la que da una densidad mayor, no 
podemos afirmar nada al respecto. 

En la Tabla 4.11 se muestran los resultados tanto de las propiedades físicas como de 
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las abundancias iónica y totales, y los ICFs. En este brote no se ha medido la línea del 
[Feíii] A4658 A, pero sí se tiene una medida de las líneas del He. De la medida global, se 
pudo determinar un valor "adoptado" de He^/H"'", derivado a partir de las líneas 4471, 5876, 
6678, 7065 A, siendo el mismo de 0.077±0.023. Para las componentes, angosta 1, angosta 
2 y ancha, se encontró un valor "adoptado" de He"'"/H+ usando sólo la medida de la línea 
5876Á, dando 0.051±0.005, O.lOliO.Ol, 0.111±0.024 respectivamente en cada caso. El valor 
estimado para la componente angosta 1 es bajo, si lo comparamos con los de otras galaxias 
de este tipo, mientras que para las otras dos componentes los valores son típicos de galaxias 
H II (ver Hágele et al., 2008). Se debe destacar que el valor "adoptado" para la medida global 
(0.077±0.023) es aproximadamente el promedio de los otros valores pesados por la luminosidad 
de las componentes, y está dentro de los valores típicos para esta clase de objetos. 

En lo que respecta a la abundancia total del oxígeno, la medida global presenta un valor de 
8.17ib0.05, 0.3 veces el valor solar, mientras que para las distintas componentes cinemáticas se 
obtuvieron 8.04±0.05 (angosta 1), 8.11±0.08 (angosta 2), y 8.31±0.13 (ancha), que son 0.22, 
0.26 y 0.42 veces la abundancia solar, respectivamente. Esto estaría indicando que, con errores 
de 1(7 (que son los dados en la tabla), la abundancia del oxígeno discriminada por componentes 
cinemáticas presenta valores diferentes entre la componente angosta 1 y la ancha, teniendo la 
componente ancha el valor más alto. La componente angosta 2 presenta el valor más próximo 
a la medida global, valor medio entre la ancha y la angosta 1. Sin embargo, si se consideran los 
errores de 2a las abundancias totales de oxígeno son las mismas para todas las componentes 
cinemáticas. La abundancia total del oxígeno en este brote derivada para cada una de las tres 
componentes y en la medida global, presenta, dentro de sus errores, los característicos valores 
bajos que se encuentran las galaxias Hii: 12+log(0/H) entre 7.94 y 8.19 (Terlevich et al., 
1991; Hoyos & Díaz, 2006). 

La abundancia total del azufre para la medida global es 0.20 veces el valor solar, mientras 
que las componentes angostas 1 y 2 dan 0.10 veces el solar y la ancha 0.32 veces el valor solar. 
En la abundancia total del nitrógeno los valores van de 0.12 a 0.19 veces el valor solar. 

Las diferentes abundancias iónicas del azufre, S"'"/H"'" y S^^/H"*", la abundancia iónica 
N"'"/H"'" y la abundancia iónica Ar^^/H"*" obtenidas para la componente ancha son más altas 
que las encontradas en las componentes angostas, entre 0.17 a 0.61 dex respecto a la angosta 
1, y entre 0.28 a 0.73 dex respecto a la angosta 2. Sólo se pudieron estimar las abundancias 
iónicas y total del [Neiii] con la medida global dado que la línea del neón en 3868 A no pudo 
ser resuelta por componentes cinemáticas. Los datos están expuestos en la Tabla 4.11. 

Respecto a las abundancias del azufre y el nitrógeno relativas al oxígeno, los valores S/0 
y N/0, la componente ancha tiene un valor del S/0 próximo al de la medida global, pero 
si se tienen en cuenta los errores, todos los valores son muy similares entre si (ver Figura 
4.4). El valor de N/0 en la componente angosta 2 es el más alto, pero su diferencia con 
las otras componentes está dentro de los errores como para considerarla significativa. Esto 
puede ser una evidencia que la evolución química de las diferentes componentes es muy similar 
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(ver discusión en Pérez-Montero et al., 2011) o que directamente son distintas fases de un 
mismo gas. Este brote presenta un exceso muy alto de N/0 (ver Figura 4.4) respecto a lo 
que generalmente se encuentra en las galaxias Hii (log(N/0) =-1.6; Pagel et al., 1979; AUoin 
et al., 1979), (ver discusión en Amorín ct al., 2010; Pcrcz-Montero et al., 2011). 

La extinición observada en el brote A es consistente con baja extinción. La componente 
angosta 2 es las que presenta mayor extinción y la componente angosta 1 es la que muestra 
menor extinción. La ancha tiene una extinción con un valor promedio de las otras dos com- 
ponentes cinemáticas, y la medida global muestra un valor similar, dentro de los errores, al 
de la componente ancha. 

Haro 15 B 

A diferencia del espectro obtenido en el modo de ranura larga del brote B, aquí sólo fue 
posible derivar dos temperaturas electrónicas con el método directo para la medida global y 
para las dos componentes cinemáticas angosta y ancha, Te([Oiii]) y Te([Siii]). Las precisiones 
alcanzadas, respectivamente para Te([Oiii]) y Te([Slli]) son: para la medida global, del orden 
de 5 % y 13 %; 15 % y 18 % en la componente angosta y 4 % y 14 % en la ancha. La temperatura 
Te([Oli]) se obtuvo de los modelos debido a la no detección de las líneas aurórales [Olí] AA 7319, 
7330 A. Y como no fue posible medir la línea débil auroral [Sil] A 4068 A, la temperatura 
Te([Sll]) se consideró igual a Te([Oii]). La línea auroral [Nii] A 5755 A tampoco en este caso 
se pudo medir en el espectro, y se ha hecho la consideración Te([Nii]) sáTe([Oii]). 

La densidad electrónica obtenida para la medida global es de 150 partículas por cm^ con 
una incerteza de 50 partículas por cm'^. Las densidades estimadas para las dos componentes 
cinemáticas están dentro del rango de bajas densidades con valores ~100 y 150 partículas por 
cm^ para la componente angosta y ancha, respectivamente. En estos dos últimos casos, los 
errores en las densidades son muy grandes, por lo cual estos valores dan una idea del orden 
de magnitud de la densidad. 

Los resultados de las diferentes abundancias iónicas y totales junto con los ICFs derivados, 
y los errores correspondientes, se muestran en la Tabla 4.11. Si se comparan las abundancias 
iónicas derivadas para cada componente cinemática se puede ver que los iones pertenecientes a 
la zona de menor excitación (O^, S^, N^) muestran una mayor abundancia en la componente 
angosta, mientras que los iones de la zona intermendia (S^"*", Ar^+) y los de la zona de alta 
excitación (O^^, Ne^"*") muestran el comportamiento opuesto, siendo mayores los valores de 
las abundancias determinados en la componente ancha. En este espectro no se pudo medir la 
línea del [Feíii] en 4658Á. 

Los valores de las abundancias totales son: la del oxígeno derivada para la medida global 
es 0.3 veces el valor solar, y la componente angosta y la ancha son 0.22 y 0.4 veces el valor 
solar, respectivamente. Estos tres valores son típicos de las galaxias H ii. Para las abundancias 
de azufre los resultados son: en el caso de la medida global 0.09 veces el solar, y 0.06 y 0.12 
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veces el valor solar en las respectivas componentes angosta y ancha. En el caso del nitrógeno 
la medida global es 0.08 veces el valor de la abundancia solar, y las componentes angosta y 
ancha son 0.06 y 0.11 veces el valor solar. La componente angosta parece tener un valor alto 
en abundancia de nitrógeno (ver valores para la muestra de galaxias Hii de Pérez-Montero 
& Díaz, 2003). Haciendo una consideración estadística con errores a 2a (los valores de los 
errores dados en la Tabla son a lo"), los valores de las abundancias totales son iguales para las 
distintas componentes y para todos los elementos, dentro de esos errores, y en algunos casos 
incluso a la. Respecto a las abundancias del azufre, nitrógeno, neón y el argón relativas al 
oxígeno, ambas componentes presentan, prácticamente valores iguales, incluso comparando 
estos valores con los estimados para la medida global. Esto podría estar indicando que ambas 
componentes presentan una misma evolución química y/o que son dos fases distintas de un 
mismo gas. 

Con la medida global, ha sido posible estimar la abundancia total de He a partir de la 
obtención de la abundancia iónica de He+ dada por la medición de las líneas Hel AA 4471, 5876, 
6678 y 7065Á. El valor "adoptado" He+/H+ que se ha obtenido es 0.080±0.017, similar al 
encontrado en la medida del mismo brote en el espectro de ranura larga y en otras galaxias H ii 
(Hágclc ct al., 2008). En el caso de la componente angosta, el valor "adoptado" se estimó con 
la medida de las líneas Hei AA 5876 y 6678 A y es de 0.043±0.006 (valor bajo para este tipo de 
objetos) y en la ancha 0.094ib0.014, valor típico para esta clase de galaxias. Este es el mismo 
comportamiento que se encontró en el brote A, entre la componente angosta 1 y las otras dos 
componentes. 

El brote B muestra baja extinción. La componente angosta presenta mayor extinción que 
la componente ancha y la medida global muestra una extinción cuyo valor es el promedio de 
los valores de las otras dos componentes. La extinción medida en el brote B observado en el 
modo de dispersión simple, considerando los errores, tiene un valor levemente menor al que 
se observa en la componente angosta observada en el modo échelle. 

HarolS C, Haro 15 E, HarolS F 

A continuación en esta subsección se analizarán conjuntamente los tres brotes restantes, 
C, E y F, analizando las medidas globales y las diferentes componentes cinemáticas, cuando se 
han podido ajustar. Los espectros de estas tres regiones tienen un bajo S/N lo cual dificultó la 
medición de no solo las líneas débiles aurórales sino de algunas líneas más intensas necesarias 
para el cálculo de temperaturas y densidades. 

En ninguno de los tres brotes, fue posible detectar ninguna de las líneas aurórales, además 
dada la baja señal de los espectros, en las medidas globales de los brotes C y F como en las 
componentes cinemáticas de los brotes C y E, no se pudieron medir las líneas del [Olí] AA3727, 
3729 A siendo imposible aplicar el método D07 para determinar la Te([Slii]) ni el método P07 
para determinar la temperatura Te([Nii]), dado que estos métodos requieren de esta líneas. 

84 



4.5. DISCUSIÓN 85 

Por tal motivo, en estos casos se consideró la temperatura electrónica de [Olii] como lO^K 
como fue comentado en la Sección 4.4.1. Con este valor se estimó la temperatura de [Olí] 
usando la relación entre ambas de Pérez-Montero & Díaz (2003), ya descripta en la misma 
Sección. Y las temperaturas Te ([Sil]) y Te([Nii]) se consideraron iguales a la del [Olí], y 
Te([Siii]) se determinó usando su relación con la temperatura de [Olii] derivada por H06. 

En cambio, en la medida global del brote E, si fue posible medir las líneas [Olí] AA3727, 
3729 A, con las cuales se pudo aplicar el método P07 para estimar la temperatura del [Nii] 
y luego con esta temperatura estimar la Te([Siii]), usando la diferencia sistemática de 500K 
entre temperaturas hallada por Díaz et al. (2007), y Te([Oiii]) del ajuste de H06. La tempera- 
tura Te ([Olí]) se obtuvo utilizando los mismos modelos y Te ([Sil]) considerando nuevamente 
su igualdad con la temperatura Te ([Olí]). 

Las densidades electrónicas estimadas de la medida global como de la componente angosta 
del brote C son de 400 partículas por cm^. En el caso de la componente ancha se encontró un 
valor mucho más alto de 1400 partículas por cm^. Estos datos dan una estimación del orden de 
magnitud de la densidad electrónica, dado que la incerteza en el cálculo es muy grande. Esto 
en parte puede deberse a la poca señal en el espectro y a que la línea [Sil] A 6731 Á debió ser 
corregida por emisión de líneas telúricas de cielo, dificultando los ajustes de la componente 
ancha (ver Tabla 4.5). La densidad estimada a partir de la medida global del brote E es de 500 
partículas por cm'^, mientras que la componente angosta 1 presenta una densidad electrónica 
de 900 partículas por cm^ y la componente angosta 2, 400 partículas por cm"^. En el caso de 
la componente angosta 1 del brote E, el alto valor de la densidad electrónica tal vez se deba 
al mismo problema presentado en el brote C. En la medida global del brote F, la densidad 
estimada es de 100 partículas por cm^. Debido a los errores muy altos, todos estos valores 
dan únicamente una idea del orden en magnitud de la densidad. 

Las abundancias iónicas en estas condiciones fueron difíciles de estimar. En las medidas 
globales de los brotes y en las dos componentes cinemáticas de los brotes C y E, se pudieron 
obtener las abundancias iónicas de O^"*", S"*" y N+. En la medida global del brote E, además 
fue posible obtener la abundancia iónica O"*" y las abundancias totales del oxígeno y azufre, 
y se pudo estimar el cociente N/0. Todos estos datos junto a las propiedades físicas de los 
brotes, están expuestos en la Tabla 4.11. En estos brotes no se midió la línea [Feíii] A4658A, 
por lo que no se pudo derivar una abundancia de este mismo elemento. 

La abundancia de O'^"'" en el brote C derivada a partir de la componente angosta tiene un 
valor más alto, ~0.3 dex, que la derivada a partir de la componente ancha, pero esto podría 
deberse a una diferencia en el grado de ionización de las componentes. Pero, para saber si 
hay una diferencia real entre las abundancias se debe comparar, obviamente, las abundancias 
totales y no las iónicas. En la medida global, el valor de 7.83ib0.03 dex es aproximadamente el 
promedio de dichas abundancias iónicas derivadas para cada componente (7.85, teniendo en 
cuenta que la escala es logarítmica). Esto se debe a que las intensidades de las líneas de emisión 
corregidas por enrojecimiento son prácticamente iguales. Las abundancias iónicas del S^ son 
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iguales, dentro de los errores, en las dos componentes y en la medida global. Y la abundancia 
de N+ es mayor en la angosta, 0.40 dex respecto de la ancha, pero esto nuevamente puede 
deberse a una diferencia en el grado de ionización de las componentes. Tanto en la abundancia 
del O^"*" como en la de N"*", la medida global es aproximadamente el valor promedio de las 
estimaciones para las dos componentes cinemáticas. 

Las abundancias del O^^, S"*" y N"*" en el brote E presentan valores más altos, derivados 
a partir de la componente angosta 2 que para la derivada a partir de la angosta 1. Para la 
medida global nuevamente se ha encontrado, como era esperable, que es aproximadamente 
el promedio pesado por la luminosidad de las dos componentes cinemáticas. La abundancia 
total del oxígeno para la medida global es 0.31 veces el valor solar. En este brote se pudo 
estimar un valor de la abundancia de helio (0.124ib0.040), que es un poco elevado pero aún 
en el rango de las galaxias Hii (Hágele et al., 2008). El mismo fue obtenido por la medida de 
la línea 5876A y fue tomado como el valor "adoptado" . 

La extinción observada en los brotes C y E es consistente con baja extinción. El brote 
C presenta valores de extinción con errores muy grandes, y solo se puede decir que las com- 
ponentes cinemáticas como la medida global presentan valores de extinción similares. En el 
brote E, la componente angosta 1 tiene un valor de extinción mayor a los valores de extinción 
de la componente ancha y de la medida global, ambos, como se ha comentado más arriba, con 
valores de c(H/3) compatibles con cero dentro de los errores observacionales, y por esta razón 
fueron igualados a cero. En el brote F las intensidades de las líneas no fueron corregidas por 
enrojecimiento. . 

Comparando la abundancia total de oxígeno derivada en el brote A con el valor obtenido 
por López-Sánchez & Esteban (2009b), 8.37±0.10, y teniendo en cuenta los errores, se puede 
ver que la abundancia obtenida para la medida global es menor a este valor. Esto puede deberse 
a las diferentes hipótesis consideradas para determinar la temperatura electrónica y por ende 
las abundancias. En el brote B, el valor obtenido en la medida global se encuentra en muy buen 
acuerdo con el valor de abundancia total de oxígeno obtenido por López-Sánchez &: Esteban, 
en el mismo brote, 8.10±0.06. Debe destacarse que para este último brote tanto el trabajo de 
López-Sánchez &; Esteban como el nuestro se realiza una derivación directa de la temperatura 
de [Oiii], y que están en excelente acuerdo (12000±600 K en este trabajo y 12900±700 K 
en el trabajo de López-Sánchez & Esteban), y una estimación de Te([Oii]) a partir de la 
de [Oiii], en el caso de López-Sánchez &: Esteban, utilizando la relación de Garnett (1992). 
Estas temperaturas también están en muy buen acuerdo, (11000ib200 y 12000±500 K, en este 
trabajo y en el de López-Sánchez & Esteban, respectivamente). 

Estudios recientes realizados en esta galaxia indican que el cúmulo estelar ionizante en el 
brote B sería el más joven de la galaxia dado que este brote presenta un color azul y una 
alta emisión FUV, indicativo de una reciente actividad de formación estelar en curso, un 
resultado apoyado por la detección de rasgos Wolf-Rayet en el espectro (López-Sánchez & 
Esteban, 2010). Esto está de acuerdo con los anchos equivalentes (EW) de la línea H/3 (ver 
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Tablas 4.2, 4.3 y 4.4) derivados de los datos de este trabajo, que reflejan que el brote B es 
más joven al tener un ancho equivalente mayor al resto de los brotes (Tcrlevich ct al., 2004). 
Por otro lado, del análisis realizado en este trabajo, se encuentra que la abundancia total de 
oxígeno derivada en los brotes A y B es muy similar. 

Si se tienen en cuenta las abundancias iónicas del oxígeno derivadas para el brote B 
utilizando los datos de dispersión simple, puede notarse que la abundancia para el O"*" es 
menor (0.16 dex) que el valor derivado de la medida global del échelle. Ambos valores son 
compatibles si consideramos los errores a 2a. Por otro lado, la abundancia de oxígeno O^^ 
derivada del espectro de ranura larga y de la medida global del échelle son compatibles dentro 
de los errores. La abundancia total de oxígeno de los datos de dispersión simple es menor 
(0.12 dex) si se la compara con la estimación global del échelle, y mucho menor (0.25 dex) si 
se la compara con la estimada en la componente ancha, pero es igual a la derivada a partir de 
la componente angosta. Sin embargo, si se tienen en cuenta los errores a 2a, el valor derivado 
de los datos de ranura larga y el global del échelle son compatibles. También se observa que 
el valor "adoptado" He^/H"'" de la medida global en el brote B del échelle, es similar al 
encontrado en el espectro de ranura larga y es compatible con los valores típicos para las 
galaxias Hii. 

Se ha comentado sobre el exceso de N/0 en el brote A, el cual es muy alto respecto a 
los valores encontrados en las galaxias Hii, aunque ya se han observado este tipo de excesos 
en otros objetos (ver discusión en Pérez-Montero et al., 2011). El brote B muestra también 
excesos en N/0 pero el mismo no es tan marcado como en el caso del brote A. Los valores 
del S/0, para los brotes A y B, aunque presentan una gran dispersión, están en el rango 
de este tipo de objetos (Hágele et al., 2008) como se puede apreciar en la Figura 4.4. El 
brote A muestra valores mayores de este cociente, que son compatibles con el valor solar 
dentro de los errores, al menos para la componente ancha. El valor del Ne/0 en el brote 
A (medida global) y en el brote B, dentro de los errores, son consistentes con el valor solar 
(log(Ne/0) = 0.61 dex; Grevesse & Sauval, 1998) y se encuentran en excelente acuerdo con 
los valores de la literatura (Hágele et al., 2006, 2008, 2011). El valor del cociente Ar/0 para el 
brote A presenta una dispersión mayor respecto de la que presenta el brote B, esto se puede 
ver en la Figura 4.4. Los valores del Ar/0 en el brote A son 0.08 y 0.39 dex mayores al valor 
solar (log(Ar/0) = -2.29 dex; Grevesse & Sauval, 1998) con la excepción de la componente 
angosta 2 que es 0.14 dex menor al valor solar. En cambio los valores obtenidos en el brote B 
son entre 0.1 y 0.2 dex menores al solar. Los valores de las abundadas relativas N/0, Ne/0 
y Ar/0 en el brote B observados en ranura larga, son consistentes con los observados en el 
modo échelle. En cambio, el valor de S/0 en el modo de dispersión simple se aparta de los 
valores del échelle, siendo mucho mayor. 

En el caso de las abundancias totales de He, se ha encontrado un comportamiento similar 
para la componente angosta 1 del brote A y la angosta del brote B respecto al resto de las 
componentes, y es que el valor "adoptado" de He^/H"*" es bajo para este tipo de objetos, pero 
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en cambio el valor obtenido para las componentes angosta 2 y ancha del brote A y para la 
componente ancha del brote B son valores típicos para las galaxias H ii. 

Respecto de los brotes C, E y F como es posible ver, el análisis comparativo sobre las 
abundancias iónicas y totales entre brotes es aún incompleto, dado que los espectros no han 
tenido la suficiente señal para realizar estas estimaciones y será necesario obtener nuevos 
datos para un análisis más profundo. Sólo se puede decir que la abundancia total de oxígeno 
del brote E para la medida global es muy similar a las abundancias obtenidas en los brotes 
AyB. 

La extinción medida en los brotes A, B, C y E, considerando los errores, presentan valores 
similares y todos compatibles con baja extinción, aunque se puede decir que el brote B tiene 
un valor levemente menor al del resto de los brotes. 



4.5.3. Estructura de ionización 

Desde hace varios años, se vienen utilizado las líneas más intensas del espectro para deducir 
metalicidades a través del uso de parámetros empíricos y modelos. Encontrar una cantidad 
basada en la intensidad de las líneas de emisión y que se correlacione con la metalicidad no 
es una tarea sencilla, dado que la mayoría de las líneas de emisión dependen de la geometría 
de la nebulosa a través del parámetro de ionización, o la "dureza" de la radiación ionizante. 
La intensidad de las líneas de emisión en el oxígeno no tiene un comportamiento lineal con 
su abundancia debido a la importancia de dicho elemento en el proceso de enfriamiento de la 
nebulosa. De este modo, al aumentar la metalicidad, el enfriamiento se vuelve más eficaz, la 
temperatura electrónica disminuye y con ella la intensidad de las líneas. Por ello las líneas más 
intensas tienen un comportamiento dual, de modo tal que a metalicidad baja son directamente 
proporcionales a la abundancia y a metalicidades altas ocurre justo lo contrario (Osterbrock, 
1989). 

El parámetro de ionización, definido como el cociente entre el número de fotones ionizantes 
y la densidad de átomos de hidrógeno, normalizado por un factor, da una idea del grado de 
ionización. Físicamente se interpreta como la velocidad máxima que puede alcanzar el frente 
de ionización dentro de la nebulosa. 

^_ Q{H) 



áncnnr'^ 



donde Q{H) es la tasa de producción de fotones de Lyman, r la distancia a la estrella o cúmulo 
central, uh, la densidad de átomos de hidrógeno neutros o ionizados y c es la velocidad de la 
luz. Su valor efectivo es, esencialmente, lo que vale en el radio de Strómgren, r^: 

U = -^^ = {Q{H)ennf" 



A-KCnHT 
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Figura 4.4: En los cuatro paneles se muestran los cocientes relativos al oxígeno: S/O, N/O, No/O y Ar/O en función 
de la abundancia de oxígeno [12 + log (O/H)]. Con diamantes están representados los objetos observados: en amarillo, 
brote B, en violeta brote A. Los objetos observados en Hágele et al., (2006; 2008; 2011) están representados por círculos 
turquesa, y las galaxias Hll de la literatura recopilados por Hágele et al., (2008) para el S/O y N/O, y por Pérez- 
Montero et al., (2007) para el Ne/O y el Ar/O con cuadrados grises. Con su habitual símbolo se representa los valores 
de las abundancias solares: del oxígeno de AUende-Prieto et al. (2001), del nitrógeno de Holweger (2001) y Grevesse & 
Sauval (1998) para el azufre, el neón y el argón. 
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donde e es el llamado factor de llenado que indica las irregularidades en la densidad de la 
nebulosa e influye en la emisividad y en la profundidad óptica. Si e << 1, como generalmente 
ocurre, la nebulosa es ópticamente delgada (Pérez-Montero, 2003). 

Para la determinación del parámetro de ionización, es bastante común usar cocientes 
de abundancias de especies en estados consecutivos de ionización. Sin embargo, a causa de 
la dificultad que tiene su determinación es igual de útil, y más fácil de medir, el cociente 
de líneas de un mismo elemento en dos estados diferentes de ionización como es el caso 
de [Oii]/[Oiii] o de [Sii]/[Siii]. Ambos tienen un comportamiento que también depende de 
la temperatura efectiva del campo de radiación (Teff) de la estrella o cúmulo ionizante. El 
cociente [Oii]/[Oiii] disminuye al aumentar Teff, al ir convirtiéndose más O"*" en O^^, y se 
considera válido en nebulosas altamente ionizadas. El cociente [Sii]/[Siii], por el contrario, 
aumenta al aumentar Teff debido a la ionización de S'^^ en S^+ y el efecto es más pronunciado 
cuando el grado de ionización es mayor. Por tanto, se utiliza más en nebulosas de baja 
ionización (log U< 2). Otros indicadores del parámetro de ionización son los cocientes entre 
líneas colisionalmente excitadas y de recombinación, como es el caso de [Sil] A 6717,6731 A / Ha 
o de [Oii]/H/3 que son mucho más independientes de la temperatura efectiva y están menos 
afectados por enrojecimiento pero que, en cambio, dependen fuertemente de la metalicidad 
(Pérez-Montero, 2003). 

Vílchez & Pagel (1988) demostraron que el cociente de las cantidades 0"'"/0^+ y S+/S^+, 
llamado "softness parameter" (cuyo significado en inglés es "parámetro de suavidad" ) y deno- 
tado con la letra griega rj, está intrínsecamente relacionado con la forma del continuo ionizante 
y depende ligeramente de la geometría. Suponiendo una geometría esférica simple y un factor 
de llenado constante, el efecto geométrico puede ser representado por el parámetro de ioni- 
zación que, a su vez, puede estimarse a partir del cociente [Oii]/[Oiii]. Y, como se vio más 
arriba, [Oii]/[Oiii] depende de Teff que, a su vez, depende de la metalicidad. El siguiente 
análisis sólo es posible realizarlo en los brotes con alto S/N en sus espectros, como lo es en 
los brotes A y B (ver Tabla 4.12). 

La contraparte puramente observacional del parámetro t] está dada por el parámetro r/' 
definido también por Vílchez & Pagel (1988) como el cociente entre las cantidades [Oii]/[Oiii] 
y [Sil] /[Sin]. La relación entre ambos parámetros estará dada por el siguiente ajuste: 

log 7]' = ¿0577 — 0,14íe — 0,16 

donde te representa la temperatura electrónica en unidades de lO^K. La relación entre r/ y ry' 
está dada a través de la temperatura electrónica, pero muy levemente, por lo que un cambio 
en la temperatura de 7000 a 14000 K implica un cambio en el logaritmo de un 0.1 dex, dentro 
de los errores de observación. Siempre el logr/' es inferior a logij. 

En el panel superior de la Figura 4.5 se muestra la relación entre log(0+/0^^) y log(S^/S^^) 
las cuales fueron derivadas usando métodos directos, modelos de fotoionización y relaciones 
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Tabla 4.12: Cocientes log(S+/S^+) y log{0+/0^+) para obtener logí?, y cocientes log ([Sn]/[Slll]) y log ([Oll]/[Oni]) 
para obtener logr;' , para las distintas componentes y medidas globales de los brotes A y B. Se agrega los correspondientes 
valores para el brote B de dispersión simple (DS). 



brote 



log(S+/S2+) log(0+/02+) log([Sii]/[Siii]) log([Oii]/[Oiii]) 



H15 A global 


-0.28±0.18 


0.12±0.07 


0.12±0.05 


0.07±0.03 


angosta 1 


-0.21±0.18 


0.32±0.06 


0.17±0.07 


0.26±0.03 


angosta 2 


-0.49±0.18 


0.12±0.12 


-0.17±0.08 


0.02±0.05 


ancha 


-0.10±0.30 


-0.01±0.19 


0.37±0.14 


-0.01±0.09 


H15 B global 


-0.63±0.13 


-0.90±0.11 


-0.45±0.04 


-l.OiO.Ol 


angosta 


-0.11±0.18 


-0.23± 0.28 


-0.07±0.05 


-0.34±0.02 


ancha 


-1.02±0.15 


-1.17±0.09 


-0.78±0.06 


-1.25±0.02 


DS 


-0.23±0.30 


-0.60±0.06 


-0.30±0.05 


-0.78±0.01 



empíricas, según corresponda el caso, para las medidas globales y las diferentes componentes 
cinemáticas de los brotes A y B (representados con diamantes violetas y amarillos, respectiva- 
mente). Por otra parte, en el panel inferior se representa el cociente entre el log([Oii]/[Oiii]) 
vs. log ([Sli]/[Slii]), el cual no requiere el explícito conocimiento de la temperaturas de línea 
involucrada en la derivación de los cocientes iónicos, y por lo tanto no depende del método 
para estimar estas temperaturas. Los objetos estudiados en Hágele et al., (2006; 2008; 2011) 
son representados con círculos turquesas, y las galaxias Hii, cuyos datos fueron obtenidos de 
la literatura con cuadrados (ver descripción y referencias en Hágele et al., 2008). Las galaxias 
Hll se localizan en la región donde el logry está entre los valores -0.35 y 0.2, correspondiente 
a valores altos de la temperatura efectiva del campo de radiación ionizante. En la Figura 4.5 
las líneas diagonales, tanto en el diagrama superior como en el inferior, representan valores 
constantes de logrj y de logr/' , respectivamente. En la Tabla 4.12 se muestran los valores de 
los cocientes log(0+/02+) y log(S+/S2+), y los cocientes log([On]/[Olii]) y log([Sii]/[Siii]) 
para las distintas componentes y medidas globales de los brotes A y B. Tanto en la Tabla 
como en la Figura 4.5, y a modo de comparación, se agregaron los valores correspondientes a 
la medida global del brote B observado con el modo de ranura larga. Todos los objetos están 
señalados con sus respectivos nombres y modos de observación. 

De la Figura 4.5 se desprende que el brote B tiene una estructura de ionización (o valor de 
r]) muy similar comparando entre si la medida global y las distintas componentes cinemáticas, 
e incluso si estos valores son comparados con los valores del espectro de dispersión simple del 
mismo brote, y siempre considerando los errores, lo que implica que la temperatura del campo 
de radiación ionizante en el brote es muy parecida para las distintas componentes cinemáticas. 
Esto podría deberse a que el cúmulo estelar responsable de la ionización de cada componente 
cinemática sería el mismo. En el brote A sucede algo similar, dado que tanto las componentes 
cinemáticas como la medida global, presentan una estructura de ionización muy parecida 
y la posición en ambos diagramas es compatible, aunque los valores en el diagrama de r] 
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Figura 4.5: Panel superior: log(0"'"/0^^") vs. log(S^/S^''') para las medidas globales y las diferentes componentes 
cinemáticas de los brotes A y B, representados con diamantes violetas y amarillos, respectivamente. Para el brote B, 
se agrega el dato de dispersión simple. Las diagonales en este diagrama corresponden a valores constantes de log 77. En 
el Panel inferior: log ([Oll]/[Oni]) vs. log ([Sll]/[Sin]), también para las mismas medidas. Las diagonales en este otro 
diagrama corresponden a valores constantes de log r¡' . En ambos diagramas, los valores están representados junto a 
objetos estudiados por Hágele et al., (2006; 2008; 2011), representados con círculos turquesas, y a los datos de las 
galaxias Hll, representados con cuadrados, y cuyos valores fueron obtenidos de la literatura (ver descripción y referencias 
en Hágele et al., 2008). El subíndice g se refiere a la medida global, el n se refiere a la componente angosta (1 o 2, para 
el brote A), el subíndice b se refiere a la componente ancha, y el subíndice sd se refiere al dato de dispersión simple 



92 



4.5. DISCUSIÓN 93 

están desplazados ligeramente hacia temperaturas efectivas menores del campo de radiación 
ionizante. La única componente que se distancia del resto del grupo, estando en la zona de alta 
ionización, es la componente ancha, pero es similar si tenemos en cuenta los errores grandes 
en la determinación de los cocientes. Esta diferencia respecto al brote B podría deberse a un 
estado evolutivo diferente dado que todas las componentes del brote A se ubican hacia una 
región con temperaturas efectivas del campo de radiación ionizante menores respecto a las 
componentes del brote B, en acuerdo con la presencia de un cúmulo estelar ionizante más viejo, 
más evolucionado, en el brote A, como fue sugerido por López-Sánchez & Esteban (2010). En 
el brote A los valores de ry y r/' son mayores, en casi todos los casos, a los típicos valores de 
las galaxias Hii (Hágele et al., 2008), solamente la componente ancha muestra valores de este 
tipo, sobre todo en el diagrama r/' . En cambio para el brote B, los valores de r/ y ?/' están en 
el rango encontrado para las galaxias H II. 

4.5.4. Abundancias químicas usando parámetros empíricos 

Las líneas de recombinación emitidas por los metales proporcionan las abundancias más 
precisas gracias a su débil dependencia con la temperatura electrónica en el interior de la 
nebulosa. Desafortunadamente, la gran mayoría de las líneas intensas de emisión que se ob- 
servan en el espectro óptico y/o que se pueden medir con suficiente precisión son de naturaleza 
colisional y sus intensidades dependen exponencialmente de la temperatura electrónica. En 
principio, esta temperatura puede ser determinada a partir de los apropiados cocientes de 
líneas. Sin embargo, estos cocientes necesitan de la detección y la medida de las líneas auróra- 
les, que, como ya se ha citado anteriormente, son intrínsecamente débiles y no son detectables 
en algunas situaciones. Es este el mismo caso que sucede en las regiones con alto contenido 
metálico, donde la eficiencia del enfriamiento ejercido por los metales provoca que dichas 
líneas sensibles a la temperatura electrónica sean indetectables, al igual que en regiones Hii 
en galaxias distantes y objetos con bajo brillo superficial donde no se las puede detectar 
debido al bajo S/N de las mismas. En estos objetos, cuya derivación directa de las tempe- 
raturas electrónicas no es posible, son utilizados diferentes calibradores empíricos basados en 
las intensidades de las líneas de emisión más conspicuas. 

Los métodos empíricos se basan en las propiedades de enfriamiento de las nebulosas io- 
nizadas, esto se traduce como una relación entre las intensidades de líneas de emisión y la 
abundancia del oxígeno. Los métodos empíricos de derivación de abundancias dependientes 
de las líneas intensas que han sido ampliamente estudiados en la literatura, se basan en la 
calibración directa de la intensidad relativa de algunas líneas brillantes de emisión versus la 
abundancia de algunos iones relevantes presentes en una nebulosa (ver por ejemplo, García- 
Lorenzo ct al., 2008; Cairos et al., 2009; García-Benito, 2009; García-Benito et al., 2010, y 
referencias en ellos). 

Entre los principales estimadores empíricos utilizados en la literatura, en este apartado se 
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analizan los siguientes parámetros empíricos usados cuando fue posible disponer de las líneas 
necesarias: O23, S23, O3N2, S3O3, N2 y Ar303. A continuación se detallan cada uno de ellos 
y luego se realiza una comparación entre los resultados obtenidos para los diferentes datos 
analizados en esta Tesis. 

Parámetro O23 

El parámetro 023, también conocido como R23, propuesto por Pagel et al. (1979) se define 

como: 

_ 1(3727, 29Á) + /(4959Á) + /(5007Á) 

^^ ~ J{H^) 

El mayor inconveniente que existe en utilizar el O23 es el comportamiento no lineal con 
la metalicidad, este parámetro tiene una naturaleza dual. El comportamiento dual de la 
calibración llevó a Skillman (1989) a hablar en términos de las ramas superior e inferior de 
la calibración y a utilizar el valor del cociente de las líneas de [Oiii] A5007A/[Nii] A 6584 A 
como un criterio para decidir a qué rama pertenece un objeto y usar la curva de calibración 
respectiva. Entonces es posible distinguir tres regiones diferentes: una rama más baja donde 
O23 aumenta con una abundancia cada vez mayor; una rama superior en la que ocurre lo 
contrario, el parámetro disminuye con una abundancia cada vez mayor y una región de cambio 
de orientación, que se produce en torno a log O23 ~ 0.8 y 12+log(0/H) ss 8.0. La dependencia 
del O23 con el parámetro de ionización da lugar a incertidumbres para la metalicidad en 
aquellos objetos que se encuentran en la rama inferior, aunque la dispersión nunca es tan 
elevada como en la zona de cambio de orientación. Para la rama inferior tanto Skillman como 
McGaugh (1991) tuvieron en cuenta la dependencia de O23 con el parámetro de ionización. 
Pero es McGaugh quien encuentra una mejor calibración basada en los modelos teóricos 
en los que la ionización es proporcionada por cúmulos estelares de diferentes metalicidades, 
incluyendo una corrección para la variación del parámetro de ionización. En este trabajo se 
calculan la abundancia de oxígeno 12+log(0/H) utilizando las calibraciones de Kobulnicky 
et al. (1999) basadas en los modelos de McGaugh. Y se utilizó el valor del cociente entre las 
líneas [Nii]A6584Á y [Oil]A3727Á y el parámetro N2 ( = [Nii]A6584Á / H/3, abajo descripto), 
como un criterio para decidir a qué rama pertenece un objeto y así poder usar la curva de 
calibración respectiva (Kcwley & Ellison, 2008). Si el logaritmo del primer cociente es menor 
a -1.2 entonces pertenece a la rama inferior, y si el logaritmo del parámetro N2 es menor a 
-1.3 (entre -1.1 y -1.3 no es posible decir nada) pertenece a la rama superior. La incerteza 
del ajuste da, para la rama inferior un error cuadrático medio de 0.13 dex y para la rama 
superior 0.19 dex. 

12 + log{0 / H)ioy,er = 12 - 4,944 + 0,767x + 0,602^2 - y(0,29 + 0,332x - 0,331x2) 
12 + log{0/H)upper = 12 - 2,939 - 0,2x - 0,237^^ - 0,305x^ - 0,0283x^- 
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y(0,0047 - 0,0221x - 0,102x^ - 0,0817x3 - 0,00717x^) 



donde 



X = log{023) 

, ,^ , , /[Oiii]4959 + 5007\ 
y = log{0,,) = log [^ ^^^^3^^^^^^ J 



Parámetro S23 

El parámetro S23 fue definido inicialmente por Vílchez & Esteban (1996) como: 

/(6717Á) + /(6731Á) + /(9069Á) + /(9532Á) 



S, 



23 



I[HP) 



y usado por primera vez como calibrador empírico de la abundancia de azufre por Díaz & 
Pérez-Montero (2000). Una de las ventajas de este parámetro es que las líneas de azufre son 
intensas, incluso para los objetos de mayor metalicidad. Además, el máximo de emisividad 
del azufre ocurre a temperaturas menores que en el caso del oxígeno, con lo cual un cambio 
de orientación en la curva del S23 frente a la abundancia del oxígeno se produce a metalicidad 
solar y no a 12+log(0/H) = 8.0 como se observa en el parámetro 023- Como las líneas de [Siii] 
cumplen con la relación teórica entre ambas de: I(9532Á) ~2.44- I(9069A), sólo se necesita 
disponer de una de ellas. 

Para calcular la abundancia de oxígeno (12+log(0/H)) se utilizó la calibración del paráme- 
tro empírico S23 dada por Pérez-Montero &: Díaz (2005): 

12 + log{0/H) = 8,15 + 1,85^05523 + 0,58(^0^523)^ 

con una dispersión de 0.2 dex en todo el rango de abundancias, disminuyendo a 0.10 dex para 
las galaxias HII. 

También es posible calcular con este parámetro, la abundancia de azufre 12+log(S/H) 
utilizando la calibración del parámetro empírico S23 dada por Pérez-Montero et al. (2006): 

12 + log{S/H) = 6,622 + l,860/oc/523 + 0,382(/oc/523)^ 

con una incerteza de 0.19 dex, pero aquí no fue usado con este fin. 

Parámetro SO23 

El parámetro SO23 es la combinación entre los parámetros S23 y O23 denominado y definido 
por Díaz & Pérez-Montero (2000). El mismo aumenta monótonamente con la abundancia del 
oxígeno hasta el régimen de abundancia sobre-solar. La calibración encontrada por Pérez- 
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Montero et al. (2005) está dada por la relación: 

12 + log{0/H) = 9,09 + l,03x- 0,23x2 
donde 






con una dispersión de 0.27 dex. 



Parámetro O3N2 

El parámentro O3N2 está definido por la relación: 

03ÍV2 = log nmii^mm 

dada por Alloin et al. (1979). Este calibrador es independiente del enrojecimiento y de la 
calibración en flujo. El uso del cociente I{Ha)/I{Hf3) en la definición indica que las líneas 
fueron medidas respecto a la línea de recombinación más próxima para minimizar así la 
dependencia con el enrojecimiento o la incerteza en la calibración en flujo. El ajuste para 
el cálculo de la abundancia de oxígeno, 12+log(0/H), que más se usa es el de Pettini & 
Pagel (2004) válida solo para 03N2<2, y con un grado de incerteza de 0.25 dex. 

12 + log{0/H) = 8,73 - 0,3203^^^2 



Parámetro S3O3 

El parámetro S3O3 fue descripto por Stasinska (2006) y está dado por: 

/(9069Á) 

^'^' - 7(^5571) 

Su calibración está dada por: 

12 + log{0/H) = 8,70 + 0,28x + 0,03x2 + 0,lx^ 

donde 

X = logSsOs 

con una desviación estándar de 0.25 dex. 
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Parámetro N2 

El parámetro N2 se define como: 

'I(6584Á)\ 



N2 = log 



I{Ha) J 



definido por Storchi-Bergmann et al. (1994). Este estimador, en contraste al O23, tiene un 
carácter univaluado en relación a la metalicidad. Además, es prácticamente independiente de 
la correción por enrojecimiento y de la calibración en fiujo de los espectros, dada la proximidad 
en longitud de onda de las dos líneas involucradas. Sin embargo, este parámetro tiene una 
alta dispersión asociada a los parámetros de la nebulosa, como por ejemplo, el parámetro 
de ionización y la temperatura de la radiación ionizante. También existe un alto grado de 
incertidumbre debido a la variación de la relación de abundancias entre nitrógeno y oxígeno, 
N/0, que a bajas metalicidades es constante al ser el nitrógeno de origen primario y no 
depender de la metalicidad. Para metalicidades mayores, en torno 12 + log(0/H) = 8.0, N/0 
aumenta con la metalicidad debido a la aparición del nitrógeno secundario. La calibración más 
usada para calcular la abundancia de oxígeno es la de Denicoló, Terlevich &: Terlevich (2002) 
la cual tiene asociada una incerteza de 0.23 dex. 



12 + log{0/H) = (9,12 ±0,05) + (0,73 ±0,10)iV2 



Parámetro ArsOg 



Este parámetro también fue definido por Stasinska (2006), en busca de las mismas buenas 
propiedades que el calibrador S3O3, y se lo define como: 

^ ^ /(7136Á) 

ArsOs = — ^ 

^ ^ /(5007Á) 

y su calibración está dada por : 

12 + log{0/H) = 8,91 + 0,34x + 0,27j;2 + 0,20j;^ 

donde 

X = logArsOs 

con una desviación estándar de 0.23 dex. 

A continuación se analizarán los resultados obtenidos para la determinación de la abun- 
dancia de oxígeno con los diferentes parámetros empíricos en los diferentes objetos estudiados. 
Se comparan estas abundancias con las derivadas previamente usando las temperaturas esti- 
madas. En la Figura 4.6 se puede ver las abundancias de oxígeno y sus incertezas para cada 
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brote observado en Haro 15 calculadas usando los diferentes calibradores empíricos para las 
diferentes componentes y la medida global. De arriba a abajo: brotes A, B, C, E y F. De 
izquierda a derecha los diferentes parámetros empíricos: S23, O23, SO23, O3N2, S3O3, N2 y 
Ar303. Con estrella azul se muestra la medida global en el échelle; con barra de error en negro 
los datos de dispersión simple; el triángulo turquesa hacia arriba se refiere a la componente 
angosta del échelle (para el brote A: el triángulo lleno representa la componente angosta 1 y 
el cuadrado verde lleno representa la componente angosta 2); el triángulo rosa hacia abajo es 
la componente ancha. Las líneas horizontales continuas representan el rango de la abundancia 
(teniendo en cuenta su error) derivada con la temperatura obtenida por el método directo 
(solo para el brote B en el modo de dispersión simple); la línea horizontal a trazos representa 
el rango de la abundancia (teniendo en cuenta su error) calculada usando la temperatura 
derivada de los modelos o de las relaciones empíricas entre temperaturas (los colores de es- 
tas líneas horizontales son los correspondientes a cada medida o componente, ya sea para 
los datos de dispersión simple (en color negro) o para las diferentes componentes (turquesa, 
verde o magenta) y para la medida global (azul) de los datos échelle). Por simplicidad, no 
se pondrá el rango de la abundancia correspondiente a la medida global en los brotes A y B 
ya que su valor es el promedio pesado por luminosidad de las abundancias derivadas para las 
distintas componentes cinemáticas. Luego, sólo se gráfica este valor para el brote E dado que 
es la única estimación previa de su abundancia que se pudo realizar. En el caso del brote C 
es posible obtener y mostrar la cantidad derivada de los datos de dispersión simple, y en el 
caso del brote F es el único caso en que no se tiene una estimación previa de la abundancia 
total de oxígeno. 

En la siguiente discusión de los resultados por brotes, siempre se tendrán en cuenta los 
errores. 

En el brote A existe una buena concordancia entre la abundancia total de oxígeno de la 
medida global determinada usando los diferentes parámentro empíricos, con la respecti- 
va abundancia derivada usando las temperaturas obtenidas de los modelos y relaciones 
empíricas (aunque no se muestra el valor con líneas horizontales en la Figura). La única 
que escapa a este resultado es la determinación de la abundancia total de oxígeno con el 
parámetro empírico N2, que presenta una diferencia de 0.14 dex respecto a la determi- 
nada con la temperatura obtenida con los modelos y relaciones empíricas si tenemos en 
cuenta la barra de error. En la componente angosta 1 existe una buena concordancia en 
la determinación de la abundancia usando los parámetros empíricos: S23, O23 y SO23, 
con la abundancia determinada usando las temperaturas obtenidas de los modelos y 
relaciones empíricas (línea horizontal discontinua turquesa). En cambio, la abundancia 
determinada para la misma componente angosta 1 usando los parámetros empítricos: 
O3N2, S3O3 y N2 y Ar303 es mayor, aunque solamente para N2 se encuentra una gran 
diferencia (^■0.25 dex) y en los otros tres la diferencia entre las barras de error es del 
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orden de 0.1 dex. La abundancia determinada para la componente angosta 2 usando los 
parámentros empíricos presenta buena concordancia con la abundancia derivada previa- 
mente usando las temperaturas obtenidas de los modelos y relaciones empíricas, excepto 
en el parámetro empírico N2, que sobreestima el valor de la abundancia de oxígeno en 
~0.15 dex. La abundancia determinada para la componente ancha usando los diferentes 
parámetros empíricos, a excepción del S23 que es levemente mayor si se tiene en cuenta 
la barra de error, son muy similares a la abundancia determinada usando las tempe- 
raturas obtenidas de los modelos y relaciones empíricas (línea horizontal discontinua 
magenta) . 

En el brote B, en el dato échelle, la abundancia de oxígeno de la medida global determinada a 
partir de los parámetros empíricos O23, O3N2 y N2 muestran un valor similar al derivado 
usando las temperaturas (aunque no se muestra este valor con líneas horizontales en la 
Figura). En cambio, la abundancia determinada a partir de los parámetros empíricos 
S23) SO23, S3O3 y Ar303 presenta valores menores a la abundancia determinada usando 
las temperaturas, siendo el S23 y el SO23 (relacionados entre sí) los que muestran las 
mayores discrepancias. La componente angosta presenta valores similares en la deter- 
minación de la abundancia de oxígeno usando los diferentes parámetros empíricos y la 
determinada derivada usando las temperaturas (línea horizontal discontinua turquesa, 
que en este caso es discontinua porque Te([Oii]) fue determinada a través de modelos). 
La abundancia de oxígeno obtenida para la componente ancha y determinada a partir de 
los parámetros empíricos es siempre menor (excepto para el O23 que da 8.17±0.13) a la 
determinada usando las temperaturas obtenidas (línea horizontal discontinua magenta) . 
Esta diferencia es mucho más marcada para el parámetro SO23, de aproximadamente 
0.75 dex (teniendo en cuenta las barras de error grandes, sobre todo en el parámetro 
empírico), aunque se observan diferencias muy grandes también para el S23 y el ArsOs 
(0.4 dex considerando las barras de error), y una diferencia de ~0.2 dex para los tres 
restantes (O3N2, S3O3 y N2). En el caso del dato de dispersión simple, los valores de 
la abundancia total de oxígeno derivada con los parámetros empíricos están en muy 
buen acuerdo con la abundancia derivada usando la temperatura obtenida por el méto- 
do directo (línea horizontal continua negra), excepto para los parámetros empíricos 
S23 y SO23 que dan valores de la abundancia de oxígeno menores, pero marginalmente 
coincidentes en el primer caso, y aproximadamente 0.1 dex inferior para el segundo, 
considerando la barra de error. 

En el brote C, los valores de las abundancias determinadas para las componentes ancha y 
angosta, y para la medida global del échelle utilizando los parámetros empíricos O3N2 y 
N2 , están en buen acuerdo con el valor de la abundancia total del oxígeno estimada para 
el mismo brote con los datos de dispersión simple. Además, concuerdan con el valor de 
la abundancia de oxígeno determinada usando las temperaturas obtenidas de modelos 
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y relaciones empíricas para los datos de dispersión simple (línea horizontal discontinua 
negra). Sólo para los datos de ranura larga, el valor de la abundancia determinada 
usando el parámetro empírico O23 presenta un valor ligeramente menor al derivado 
usando las temperaturas. 

En el brote E existe una buena concordancia entre las abundancias determinadas con los 
parámetros empíricos O3N2 y N2 para las dos componentes angostas y la medida global, 
y coinciden con la abundancia determinada con las temperaturas estimadas de los mo- 
delos y relaciones empíricas (línea horizontal discontinua azul) , excepto para la angosta 
2 con el parámetro N2. La abundancia determinada con el parámetro O23 en la medida 
global está un poco por debajo de los demás, pero en buen acuerdo, también con la 
abundancia total de oxígeno determinada a partir de las temperaturas. 

En el brote F, la estimación de la abundancia para la medida global con los parámetros 
empíricos O3N2 y N2 es similar, mostrando valores muy parecidos a los encontrados 
para los demás brotes de formación estelar. 

Se debe tener en cuenta que las relaciones empíricas usadas aquí como para los parámetros 
1] y 1]^ fueron derivadas a partir de medidas para regiones completas y que nunca se ha 
hecho un análisis y/o determinación por componentes. Este es un trabajo que merece una 
discusión y un tratamiento más completo que está más allá de lo que se puede realizar con los 
presentes datos. Es necesario realizar nuevas observaciones y ampliar la muestra de regiones 
con determinaciones de abundancias por componentes y también, dentro de lo posible, con 
una determinación directa de la abundancia iónia de O"*", lo que daría una medida directa de 
O para las distintas componentes. 
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P IffUircl 4t.01 Abundancias de oxígeno y sus incertezas para cada brote observado en Haro 15 calculadas usando los diferentes parámetros 
empíricos para las diferentes componentes y la medida global. De arriba a abajo: brotes A, B, C, E y F. De izquierda a derecha los diferentes 
parámetros empíricos: S231 O23, SO23, O3N2, S3O3, N2 y Ar303. Con estrella azul se muestra la medida global en el échelle; con barra 
de error en negro los datos de dispersión simple; el triángulo turquesa hacia arriba se refiere a la componente angosta del échelle (para el 
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4.6. Resumen y conclusiones 

En este Capítulo se ha realizado un análisis de las características físicas del gas ionizado 
discriminado para cada componente y para el flujo global de las líneas de emisión del gas en los 
brotes más brillantes de formación estelar de la galaxia BCD Haro 15. Para esto se estimaron 
temperaturas y densidades electrónicas, y se obtuvieron abundancias iónicas y totales de las 
diferentes especies. 

Como en las relaciones entre temperaturas electrónicas no fue posible obtener una estima- 
ción realista de las incertidumbres asociadas a los modelos dado que los errores de los modelos 
no se pueden cuantificar, los errores que se han dado en las abundancias son cotas inferiores. 

De los espectros de los brotes B y C obtenidos en modo dispersión simple tomados con la 
cámara WFCCD del telescopio de 100" du Pont de LCO, Chile, se puede resumir los siguiente 
resultados: 



en el brote B, cuatro temperaturas electrónicas, Te([Oiii]), Te([Oii]), Te([Siii]) y Te([Sii]) 
fueron derivadas de mediciones directas con precisiones del orden de 2%, 3%, 10% y 
14% respectivamente, y se estimó la Te([Nii]) considerado la aproximación Te([Nii]) ~ 
Te ([Olí]). Mientras que en el brote C, para determinar las temperaturas electrónicas 
Te([Oiii]), Te ([Olí]), Te([Siii]), Te([Sii]) y Te([Nii]) se recurrió tanto a los modelos co- 
mo a las relaciones empíricas. Es de notar que en el brote B, Te([Oii]) es 4100 K más 
alta que Te ([Sil]). En cuanto a la temperatura de [Nii] en el brote B, el valor obtenido 
es mayor que la del brote C, aún estimándola usando el método P07. 

la densidad electrónica calculada a partir del cociente de líneas del [Sil] para ambos 
brotes está por debajo de la densidad electrónica crítica por desexcitación colisional, tal 
como sucede en los procesos de formación estelar en las galaxias H ii. 

para el brote B, se dedujeron las abundancias iónicas del O"*", O^"*", S"*", S^"*", N"*", Ne^+, 
Ar^"*", Ar^"*", Fe^"*", He"*" y He^"*". También en este brote se calcularon las abundancias 
totales del O, S, N, Ne, Ar, Fe y He, y los cocientes logarítmicos N/0, S/0, Ne/0 y 
Ar/0. En el brote C, se estimaron las abundancias iónicas O^, O^^, S^ y N^, y se 
calcularon las abundancias totales del O y N, y la abundancia relativa N/0. El brote B 
tiene una abundancia de oxígeno 12 + log (0/H) = 8.04ib0.03 (0.22 veces el valor solar), 
similar, dentro de los errores, al encontrado por López-Sánchez & Esteban (2009b) en 
el mismo brote. La abundancia de azufre en este brote es 6.52ib0.22 (0.15 veces el valor 
solar) y la abundancia de nitrógeno es 6.88ib0.20 (0.09 veces el valor solar). El valor 
estimado para N/0 es un poco mayor en el brote C que en el B, aunque igual dentro de 
los errores, mostrando en ambos casos una sobre abundancia respecto a lo encontrado, 
generalmente, para las galaxias H ii. 
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■ se ha encontrado que es mayor el grado de excitación en el brote B que en el C. 

■ la abundancia total del He estimada en el brote B tiene un valor típico, como los que 
se encuentran para las galaxias H ii. 

De los espectros de alta resolución tomados con el espectrógrafo échelle del telescopio de 
100" du Pont de LCO, Chile, de los brotes A, B, C, E y F se puede resumir los siguiente 
resultados: 



■ en el brote A, para la medida global y para las dos componentes angostas y la com- 
ponente ancha, se utilizó la relación encontrada por H06 entre temperaturas Te[Siii] 
vs. Te[Olll] para estimar la temperatura electrónica del [Olll], previa determinación de 
Te([Siii]) con el método D07. Luego, siguiendo el análisis clásico, se obtuvieron de los 
modelos las temperaturas Te ([Olí]), Te ([Sil]) y Te([Nii]). 

En el brote B, se derivaron Te([Oiii]) y Te([Siii]) con el método directo para la medida 
global y para las dos componentes cinemáticas, con precisiones del orden de 5 % y 
13% para la medida global, 15% y 18% en la componente angosta y 4% y 14% en la 
ancha. La temperatura Te ([Olí]) se obtuvo de los modelos y las temperaturas Te ([Sil]) 
y Te([Nii]) se consideraron iguales a Te([Oii]). 

En las medidas globales de los brotes C y F como en las componentes cinemáticas de los 
brotes C y E, se consideró la temperatura electrónica de [Oiii] como lO^K. Luego, con 
este valor se estimó la temperatura de [Olí] usando la relación entre ambas de Pérez- 
Montero & Díaz (2003), las temperaturas Te ([Sil]) y Te([Nii]) se consideraron iguales 
a la temperatura de [Olí], y finalmente Te([Siii]) se determinó usando su relación con 
la temperatura de [Oiii] derivada por H06. En la medida global del brote E, se pudo 
aplicar el método P07 para estimar la temperatura del [Nii], y con esta temperatura se 
estimó Te ([Sin]) haciendo uso de la diferencia sistemática de 500K entre temperaturas 
hallada por D07. Luego, Te([Olll]) se estimó usando el ajuste de H06. La temperatura 
Te ([Olí]) se obtuvo utilizando los modelos y se consideró que Te ([Sil]) es igual a la 
temperatura de [Olli]. 

■ las densidades electrónicas estimadas tanto para las medidas globales de los cinco brotes 
como para las diferentes componentes cinemáticas de los brotes A, B, C y E fueron 
calculadas a partir del cociente de líneas del [Sil] y están por debajo de la densidad 
electrónica crítica por desexcitación colisional. Estos valores dan una estimación del 
orden de magnitud de la densidad, dado que la incerteza en el cálculo es muy grande, 
salvo en el caso del brote B, donde la densidad electrónica obtenida para la medida 
global es de 150 partículas por cm^ con una incerteza de 50 partículas por cm'^. 
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■ diferentes abundancias iónicas y totales junto con los ICFs fueron derivados tanto para 
las medidas globales de los cinco brotes como para las diferentes componentes cinemáti- 
cas de los brotes A, B, C y E. En el brote B, los iones pertenecientes a la zona de menor 
excitación (O^, S^, N"*") muestran una mayor abundancia en la componente angosta, 
mientras que los iones de la zona intermendia (S^^, Ar^"*") y los de la zona de alta 
excitación (O^"*", Ne^"*") muestran el comportamiento opuesto. 

Las abundancias totales de O, S, N, Ne y Ar, donde fue posible obtenerlas, se encuen- 
tran en el mismo rango de las abundancias totales medidas para las galaxias H il. La 
abundancia total de oxígeno derivada en el brote B está en muy buen acuerdo con el 
valor obtenido por López-Sánchez & Esteban (2009b), pero en el brote A las abundan- 
cias aquí derivadas son menores (~0.2 dex) al valor encontrado por López-Sánchez & 
Esteban, aunque si se tienen en cuenta las barras de error la diferencia se reduce a ~0.05 
dex. Del análisis realizado se encuentra que, dentro de los errores, la abundancia total 
de oxígeno derivada en los brotes A, B y la medida global del brote E son muy similares. 

Si se comparan los valores derivados en el brote B para los datos échelle con los estimados 
a partir de los datos de ranura larga se puede apreciar que la abundancia total de 
oxígeno es menor para los datos de dispersión simple comparándola con la estimación 
global (aunque muy parecido se considera los errores), y más aún si se la compara con 
la componente ancha, pero es igual al derivado para la angosta. 

Se encontró un exceso de N/0 en el brote A y B, siendo en A muy alto respecto a los 
valores encontrados en las galaxias Hii, aunque este tipo de excesos ya fue observado en 
otros casos (ver discusión en Pérez-Montero et al., 2011). Los valores del S/0, para los 
brotes A y B, aunque presentan una gran dispersión, están en el rango en el que suelen 
encontrarse este tipo de objetos. Y los valores del Ne/0 en el brote A (medida global) 
y en el brote B, dentro de los errores, son consistentes con los valores de la literatura. 
Finalmente, el valor del cociente Ar/0 para el brote A presenta una dispersión mayor 
respecto de la que presenta el brote B, y en ambos casos dentro de los valores típicos para 
esta clase de objetos. La similitud, dentro de los errores, en las abundancias relativas 
S/0 y N/0 en las tres componentes del brote A podría ser una evidencia que la evolución 
química de las diferentes componentes es muy similar o que directamente son distintas 
fases de un mismo gas. Una situación similar ocurre con las abundancias relativas S/0, 
N/0, Ne/0 y Ar/0 de las dos componentes del brote B. Los valores de las abundadas 
relativas N/0, Ne/0 y Ar/0 en el brote B observados en ranura larga, son consistentes 
con los observados en el modo échelle, salvo el valor S/0 que es mayor en dispersión 
simple que en échelle. 

■ se encontró que las abundancias totales de He tienen un comportamiento similar para 
la componente angosta 1 del brote A y la angosta del brote B respecto al resto de las 
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componentes, y es que el valor "adoptado" de He"'"/H+ es bajo para este tipo de objetos, 
pero en cambio el valor obtenido para las componentes angosta 2 y ancha del brote A, y 
para la componente ancha del brote B son valores típicos de los encontrados para galaxias 
H II. Sin embargo, se necesitan datos mejores para estudiar este comportamiento y tratar 
de reducir los errores de estas estimaciones. 

■ la extinción medida en los brotes A, B, C y E, considerando los errores, presentan valores 
similares y todos compatibles con baja extinción, aunque se puede decir que el brote B 
tiene un valor levemente menor al del resto de los brotes. 

■ será necesario obtener nuevos datos de los brotes C, E y F para un análisis más profundo 
de las abundancias iónicas y totales. 

■ analizando el brote A y el brote B por separado, y considerando los errores, las diferentes 
componentes cinemáticas de cada brote presentan una estructura de ionización similar 
entre ellas, lo cual implicaría que la temperatura del campo de radiación ionizante 
dentro de cada brote es muy parecida para las distintas componentes cinemáticas. Esto 
podría deberse a que el cúmulo estelar responsable de la ionización de cada una de las 
componentes sería el mismo. Para el brote B, este resultado es válido tanto para los 
datos échelle como para los datos de dispersión simple. 

■ las abundancias totales de oxígeno estimadas con los diferentes parámetros empíricos, 
en general, muestran valores consistentes con las abundancias derivadas usando las 
temperaturas electrónicas obtenidas ya sea de los modelos y relaciones empíricas o por 
métodos directos, excepto en determinados casos donde los valores de las abundancias 
determinadas con algunos parámetros dan valores que sobreestiman o subestiman la 
abundancia total de oxígeno obtenida previamente. 

Cabe destacar que las relaciones empíricas usadas aquí como las usadas para los paráme- 
tros r/ y r/' fueron derivadas a partir de medidas para regiones completas y que nunca 
se ha hecho un análisis y/o determinación por componentes. Este es un trabajo que 
merece una discusión y un tratamiento más completo que está más allá de lo que se 
puede realizar con los presentes datos. 
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Conclusiones y Trabajo a futuro 



En esta tesis he identificado y anaUzado las características físicas y cinemática de una 
muestra de Regiones Hii Gigantes en galaxias observables con telescopios del Hemisferio Sur. 
También he definido la naturaleza y las propiedades físicas de las componentes gaseosas en 
algunos brotes de formación estelar masiva de la galaxia BCD Haro 15, permitiendo analizar 
las abundancias químicas iónicas y totales de las mismas. Para llevar a cabo este trabajo se 
han utilizado técnicas de espectroscopia de alta resolución como lo es la espectroscopia échelle, 
y también he utilizado, para una galaxia, espectros de dispersión simple de baja resolución. 

Analizando 11 regiones Hii luminosas en galaxias espirales y en una galaxia BCD, he 
obtenido los siguientes resultados generales importantes: 



v^ he encontrado comportamientos cinemáticos complejos en seis regiones de formación 
estelar de las galaxias espirales NGC 6070, NGC 7479 y en cinco brotes de formación 
estelar violenta en la galaxia compacta azul Haro 15. En cuatro de ellas se hallaron 
indicios de la presencia de componentes con velocidades radiales distintas, prácticamente 
imposibles de detectar espacialmente. 

^ he realizado un estudio detallado de las propiedades físicas de los principales brotes de 
formación estelar de la galaxia BCD Haro 15 utilizando espectroscopia de ranura larga 
y échelle obtenida en LCO. 

Dentro de esos resultados puedo decir que: 



Fue posible determinar la dispersión de velocidades dada por el ensanchamiento del 
perfil de líneas de emisión nebular para cada región. La alta relación S/N, junto a 
la resolución del espectro échelle permitió resolver el perfil de las líneas de emisión y 
calcular la dispersión de velocidades del gas ionizado. Este análisis fue hecho midiendo 
el ancho observado del perfil de las líneas de recombinación de hidrógeno después de 
corregir el mismo por el perfil instrumental y por la contribución térmica. Fue posible 
confirmar la naturaleza gigante de las regiones candidatas. 
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De la muestra de seis GHllR en las galaxias NGC 6070 y NGC 7479, se ha encontrado 
que todas ellas muestran evidencia de un residuo presente en los perfiles de las líneas 
de emisión. En estas regiones se ha intentado ajustar una componente ancha o dos 
componentes simétricas separadas en velocidad respecto a la componente principal. 
Todas las regiones estudiadas presentan componentes cinemáticas supersónicas. 

Los cinco brotes de la galaxia Haro 15, a excepción del brote F, muestran múltiples 
componentes supersónicas. En el brote F el ensanchamiento del perfil integrado sola- 
mente pudo ser detectado en la línea Ha. En esta línea de emisión se pudo ajustar una 
componente ancha supersónica de 22 kms~^ y una angosta subsónica de 8 kms~^, 
dispersión de velocidades típica de regiones Hii clásicas, aunque la presencia de una 
componente ancha sea casi exclusiva de las Regiones H ii Gigantes. 

Si bien las características de estas Regiones Hii Gigantes presentan una correlación 
en el plano log(L) — log(o"), la presencia de múltiples componentes cinemáticas infiuye 
fuertemente en la ubicación final en dicho plano. Esto representa un alerta importante 
ante el uso indiscriminado de dicha regresión cuando no se tienen datos completos 
sobre el comportamiento dinámico de los brotes de formación estelar. Las componentes 
angostas identificadas en los perfiles complejos siguen la relación de sistemas virializados, 
en cambio los perfiles ajustados con una componente Gaussiana simple muestran una 
pendiente más plana. 

El análisis de las observaciones realizadas por Hágele y colaboradores (2006; 2008; 2011) 
y la metodología definida por ellos para obtener abundancias elementales precisas de O, 
S, N, Ne, Ar y Fe en el gas ionizado, me ha permitido estimar en los brotes de la galaxia 
Haro 15 por lo menos cuatro temperaturas de línea: Te([Oiii]), Te([Siii]), Te([Oii]) y 
Te ([Sil]), y una densidad electrónica, Ne([Sli]) para los datos de ranura larga. Para 
los datos échelle, cada brote fue analizado por separado, estudiando las temperaturas, 
densidades, abundancias iónicas y totales estimadas para la medida global y para cada 
componente cinemática separadamente. En el brote B fue el único brote para el que se 
derivaron Te([Oiii]) y Te([Siii]) con el método directo para la medida global y para las 
dos componentes cinemáticas, con precisiones del orden de 5% y 13% para la medida 
global, 15 % y 18 % en la componente angosta y 4 % y 14 % en la ancha. En el resto de 
los brotes (A, C, E y F) las temperaturas electrónicas fueron estimadas usando, según el 
caso, los modelos, relaciones entre temperaturas, métodos empíricos o incluso en algunos 
casos cuando no fue posible derivar ninguna temperatura, se consideró la temperatura 
electrónica de [Oiii] como lO^K y de ella, usando relaciones entre temperaturas, se 
obtuvo el resto. Las densidades electrónicas Ne([Sii]) fueron calculadas a partir del 
cociente de líneas del [Sil] y todas están por debajo de la densidad electrónica crítica 
por desexcitación colisional. 
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Con estas medidas y un tratamiento cuidadoso y realista de los errores observacionales, 
se obtuvieron estimaciones de las abundancias químicas iónicas y totales del O, S, N, Ne, 
Ar, Fe y He, según fuere el caso. Este análisis detallado de las abundancias fue realizado 
en cada brote de formación estelar con más de una componente y en cada una de las 
componentes individuales, ya sea en los datos de dispersión simple como en los datos 
espectroscópicos échelle. Las abundancias totales de O, S, N, Ne y Ar, donde fue posible 
obtenerlas, se encuentran en el mismo rango de las abundancias totales medidas para 
las galaxias H ii. La abundancia total de oxígeno derivada en el brote B está en muy 
buen acuerdo con el valor obtenido por López-Sánchez & Esteban (2009b), pero en el 
brote A las abundancias aquí derivadas son menores (~0.2 dex) al valor encontrado por 
López-Sánchez & Esteban, aunque si se tienen en cuenta las barras de error la diferencia 
se reduce a ~0.05 dex. Del análisis realizado se encuentra que, dentro de los errores, la 
abundancia total de oxígeno derivada en lo brotes A, B y la medida global del brote E 
son muy similares. Como en las relaciones entre temperaturas electrónicas no fue posible 
obtener una estimación realista de las incertidumbres asociadas a los modelos dado que 
los errores de los modelos no se pueden cuantificar, los errores que se han dado en las 
abundancias son cotas inferiores. 

Se encontró un exceso de N/0 en el brote A y B, siendo en A muy alto respecto a los 
valores encontrados en las galaxias H ii, aunque este tipo de excesos ya fue observado en 
otros casos (ver discusión en Pérez-Montero et al., 2011). Los valores del S/0, para los 
brotes A y B, aunque presentan una gran dispersión, están en el rango en el que suelen 
encontrarse este tipo de objetos. Y los valores del Ne/0 en el brote A (medida global) 
y en el brote B, dentro de los errores, son consistentes con los valores de la literatura. 
Finalmente, el valor del cociente Ar/0 para el brote A presenta una dispersión mayor 
respecto que la del brote B, y en ambos casos el valor del cociente está dentro de los 
valores típicos para esta clase de objetos. La similitud, dentro de los errores, en las 
abundancias relativas S/0 y N/0 en las tres componentes del brote A podría ser una 
evidencia que la evolución química de las diferentes componentes es muy similar o que 
directamente son distintas fases de un mismo gas. Una situación similar ocurre con las 
abundancias relativas S/0, N/0, Ne/0 y Ar/0 de las dos componentes del brote B. Los 
valores de las abundadas relativas N/0, Ne/0 y Ar/0 en el brote B observados en los 
datos espectroscópicos obtenidos con ranura larga, son consistentes con los observados 
en el modo échelle, salvo el valor S/0 que es mayor en dispersión simple que en modo 
échelle. 

Se encontró que las abundancias totales de He tienen un comportamiento similar para 
la componente angosta 1 del brote A y la angosta del brote B respecto al resto de las 
componentes. El valor "adoptado" de He"'"/H"'" que se obtuvo es bajo para este tipo 
de objetos. En cambio, el valor obtenido para las componentes angosta 2 y ancha del 
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brote A, y para la componente ancha del brote B son valores típicos a los encontrados 
para las galaxias H ii. Es necesario obtener nuevos y mejores datos para estudiar este 
comportamiento, además de tratar de reducir los errores de las estimaciones de las 
abundancias iónicas y totales de helio. 

■ La extinción medida en los brotes A, B, C y E, considerando los errores, presentan 
valores similares y todos compatibles con baja extinción, aunque se puede decir que el 
brote B tiene un valor levemente menor al del resto de los brotes. 

■ En los brotes A y B de Haro 15, la estructura de ionización y, por ende la temperatura 
del campo de radiación ionizante, es muy similar para cada componente del respectivo 
brote estudiado. Esto implicaría que el cúmulo ionizante de cada una de las componentes 
del brote en cuestión sería el mismo. La excepción estaría en la componente ancha del 
brote A que difiere levemente de su grupo pero es similar si consideramos los errores. 
El brote B tiene una ionización más alta que el brote A. 

■ Para cada brote de Haro 15, y cuando fue posible disponer de las líneas más intensas, 
se determinó la abundancia de oxígeno a partir del análisis de los diferentes parámetros 
empíricos: O23, S23, SO23, O3N2, S3O3, N2 y Ar303. Las abundancias totales de oxígeno 
estimadas con los diferentes parámetros empíricos, en general, mostraron valores con- 
sistentes con las abundancias derivadas usando las temperaturas obtenidas ya sea de los 
modelos y relaciones empíricas o por métodos directos, excepto en determinados casos 
donde los valores de las abundancias determinadas con algunos parámetros dan valores 
que sobreestiman o subestiman la abundancia total de oxígeno obtenida previamente. 
Cabe destacar que las relaciones empíricas usadas tanto para los parámetros empíricos 
como para los parámetros r¡ y i]^ fueron derivadas a partir de medidas para regiones 
completas y que nunca se ha hecho un análisis y/o determinación por componentes. Es- 
te es un trabajo que merece una discusión y un tratamiento más completo que está más 
allá de lo que se puede realizar con los presentes datos. 

Los resultados de esta Tesis amplían de distintas formas el estudio de las Regiones Hii 
Gigantes. Si bien no es la primera vez que se combina la información cinemática - que per- 
mite separar componentes indistinguibles en espectroscopia de baja resolución - junto con 
un análisis detallado de las propiedades físicas de cada componente, este trabajo presenta 
una mayor completitud frente a otros trabajos publicados dada la gran cantidad de líneas 
analizadas, entre las que se encuentran las débiles líneas aurórales. Además, a esto se suma la 
posibilidad de descomponer sus perfiles en más de una componente angosta y una componente 
ancha. Para ello, he debido avanzar sobre las técnicas convencionales de análisis espectral, 
desarrollando una metodología de trabajo original, necesaria para llevar esta investigación a 
un resultado convincente. 
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Trabajo futuro 

La influencia de las estreUas masivas deja su impronta en la fase gaseosa, principalmente 
excitando, ionizando y modificando estructuralmente el ISM a partir de su radiación UV y 
vientos. La intensidad de las líneas de emisión del gas ionizado permite analizar las propieda- 
des físicas de dichas regiones, su dinámica interna, y caracterizar, utilizándolas como sondas, 
la cinemática de las galaxias que las albergan. Siguiendo esta línea y teniendo en cuenta los 
progresos más recientes en este campo, especialmente la nueva información que ofrecen los 
grandes telescopios, el plan de trabajo a futuro se estructura en dos objetivos principales: 



(•■ Continuar la identificación y consiguiente análisis de Regiones Hii Gigantes en galaxias 
cercanas. La identificación es posible mediante espectroscopia de alta resolución que 
permite verificar la naturaleza supersónica de la cinemática del gas ionizado. Por otro 
lado, estudiar la estructura cinemática de dichas regiones con espectroscopia 3D que 
permita combinar información cinemática y espacial, con una unívoca determinación de 
los nodos individuales que contribuyen al perfil global observado. 

(•" Definir la naturaleza y propiedades físicas de los brotes de formación estelar masiva, 
como así también entender los mecanismos que permiten su retroalimentación y la 
formación de nuevas generaciones de estrellas que producen a la vez, el enriquecimiento 
químico del Universo Local. El análisis de los parámetros físicos tales como densidad 
electrónica, temperatura electrónica, y composición química permitirá caracterizar los 
brotes de formación estelar más extremos analizando su entorno. 

Se dispone de espectros de alta resolución obtenidos con los telescopios de 2.5 y 6.5m del 
Observatorio de Las Campanas (LCO) para la detección de Regiones Hll Gigantes en más 
galaxias espirales y BCDs. Recientemente, se re-observaron dos de los brotes de la galaxia 
Haro 15 (A y B) con un mejor cociente S/N, con lo que se espera poder detectar las débi- 
les líneas aurórales y poder realizar ajustes en estas líneas que sean independientes de las 
soluciones de las líneas más intensas. 

Para estudiar comportamientos cinemáticos complejos en las regiones de formación estelar 
haré uso de las facilidades instrumentales de espectroscopia 3D de los Telescopios Gemini, 
tanto en el rango óptico (GMOS), como en el infrarrojo cercano (NIFS). Cuento con obser- 
vaciones asignadas con alta prioridad (Banda 1) para obtener espectroscopia 3D de campo 
integrado en una región de formación estelar que presenta más de una componente en el perfil 
de sus líneas de emisión. La propuesta aceptada plantea hacer uso de la alta calidad de imáge- 
nes de los telescopios Gemini en su modo IFU (sigla en inglés de Integral Field Unit) para 
mapear el campo de velocidades. De esta manera espero dilucidar si el perfil, aparentemente 
supersónico de algunas regiones HIT gigantes puede resolverse en varias componentes, cada 
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una de ellas con distintas velocidades. También, junto a Guillermo Bosch y Guillermo Hágele, 
obtuve tiempo de observación durante el semestre 2009A, con el cual planteamos mejorar la 
resolución estudiando regiones visibles con Gemini Norte utilizando la resolución espacial al 
límite de difracción en el infrarrojo cercano, utilizando óptica Adaptable con estrella guía 
láser en el instrumento GNIRS. El uso de óptica Adaptable permite obtener imágenes de la 
misma calidad que la del Telescopio Espacial (FWHM «0.08 segundos de arco) en un te- 
lescopio de 8 metros de diámetro. Esto nos permitirá estudiar la estructura cinemática de 
dichas regiones con una resolución sin precedentes necesaria para distinguir la presencia de 
distintas componentes, con velocidades diferentes, que puedan contribuir al ancho observado 
en el perfil de la emisión global. 

Se pretende ampliar la muestra tanto de Regiones Hii Gigantes como de galaxias Hii para 
las que se pueda obtener una abundancia precisa utilizando lo que se conoce como el método 
directo. Usaremos los datos estudiados en Firpo et al. (2005, 2010) y nuevas observaciones de 
regiones Hil en otras galaxias BCDs. Esto permitirá producir y calibrar relaciones empíricas 
entre las diferentes temperaturas de línea que reemplazarían a las que se utilizan actualmente 
y que están basadas en secuencias de modelos de fotoionización muy simples, y pobremente 
contrastados, debido a la falta de datos idóneos. Para ello, se requiere el uso de espectrógrafos 
de amplia cobertura espectral y resolución espectral intermedia. La muestra debe ser diseñada 
para abarcar objetos en los extremos bajo y alto de temperaturas electrónicas. Esta cobertura 
total en temperaturas electrónicas proveerá una limitación muy necesaria de la estructura en 
temperaturas de las galaxias Hii. 
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